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CHAPTER I 


Electric Multipole Polarization in Atoms and Ions: 
Theoretical Methods and Ex''^erimental Evidences for 
Sternheimer Polarization Effects 
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I.l Introduction 

The study of polarization of, electron-distribution in 

atoms in presence of a perturbing electric or magnetic field 

still continues to be one of tlie most important research areas, 

although tbe early experimental measurements of electric dipole 

polarizabilities ' " in terms of dielectric constants for the 

noble gas systems were carried out in the twenties. The lowest 

nonlinear contribution to the induced electric dipole moment 

in helium has novr been measured using laser sources to produce 

strong fields in optical Kerr-effect studies. The theoretical 
4 

studies on one-electron systems in presence of external fields 
were undertaken immediately after the formulation of new quantum 
mechanics. For the treatment of multielectron atomic systems - 
where exact results cannot be obtained due to limitations in 
theory - several approximate methods^ with varying degrees of 
accuracy have been put forward. 

The perturbing field can be internal, e.g. originating 
from the valence electrons, the electric quadrupole and hexa- 
decapole moment of the nucleus and/or external electric field. 

This is the basis of Sternheimer ‘ s pioneering work on the 

6 7 8 

magnetic dipole, nuclear electric quadrupole and hexadecapole 

polarizations commonly known as Sternheimer shielding- anti- 
shielding factors. Their knowledge becomes essential in 
interpreting the experimental data on hyperfine interactions 
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which have largely become available during recent years via 
Y -resonance, magnetic -resonance and perturbed angular-correlation 
studies. In fact, the shielding-antishielding factors represent 
our lac]c of Imov/ledge of accurate electronic wave fionctions. 

In this thesis vre shall be concerned with the Sternheimer 
shielding-antishielding factors in a variety of atoms and ions 
that are of importance both from the theoretical and experimental 
points of vie’w. 

The different approximate methods used to calculate the 
shielding-antishielding factors are siniilar to those used for 
the calculation of polarizabilities and in this chapter we 
shall review these methods together, after presenting the 
general definition of static multipole-polarizabilities and 
shielding factors. As always, the accuracy of a given theore- 
tical method will, be judged from its ability to interpret the 
experimental observations. The dipole shielding factor for 
an N-electron atomic system is a Icnoxm quantity given by^ N/Z, 
vjhere Z is the atomic number and this provides an independent 
means of assessing the accuracy of the- arproxinate method Tised. 

1.2 Multipole Polarizabilities and Shielding Factors 

Consider the physical situation in which a distant point 
charge q is located at R a.u. from the nucleus of a non- 

degenerate N-electron atom, along the Z axis. The pertxirbation 
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of interest to us here can be expressed, as a symmetric sum of 
one electron operators. The perturbing potential (in a.u. ) can 
be expressed in terms of the electron coordinates with respect 
to the nucleus as: 


V, 





(I.l) 


v;here L varies from 1 to “ and x = The total wave 

function for the atom under the influence of can be written 


to the first order in X ; 


$ = (1.2!) 

The 2^-pole polarizability is defined as (-l)^L.' times the 
ratio of the induced 2^-pole moment to the Lth-order gradient 
of and is given by 

a = -2 < I rr p (cos 0 . ) I $ > (1.3) 

J-j ’ . - X Ju X ’ 

1=1 

The 2^-pole shielding factor is given by the ratio of the 
change in Lth-order electric field gradient at the nucleus due 
to the electron distribution to the Lth-order gradient due to 
external charge alone. This i-s given by 

= 2 < I I (Pj^(cos 8 j_) ($° > (1.4) 

i=l 

The polarizability largely depends on the outer parts of the 
\vave function, v/hereas the shielding-anti shielding factors 
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largely depend on the description of the wave function near 
the nucleus. The 2^-pole polarization of the spherical electron 
distribution can be given a physical interpretation in terms of 
angular and radial distortions. The angular distortion leads 
to a shielding effect (positive) and significantly contributes 
to the polarizability. The radial distortion produces 
antishielding (negative) effect and contributes to polarizability 
much less significantly. Thus the total shielding can be either 
positive (shielding) or negative (anti shielding) depending upon 
the relative magnitudes of the angular and radial contributions. 

The static dipole, quadrupole and hexadecapole polariza- 
bilities vjill be denoted hereafter as and respectively 

(vide eqn. 1.3) and the associated Sternheimer shielding- 
antishielding factors due to distant point charges be denoted 
as 6^/ Yoo and respectively (vide eqn. 1.4). The shielding- 
antishielding factors corresxxjnding to the valence shells at 

7 3 . 8c3. 

the nuclear site for quadrupole and hexadecapole cases will 

be denoted by R and respectively. The "shielding factor 

corresponding to the quadrupole term in the crystalline electric 

field^*^"^^ at tlie valence electron site and the atomic 

13 

quadrupole moment will be respectively given by and R . 

The various interactions of the nucleus with the surround- 
ing charges are always measured in terms of the product of a 
nuclear and electronic property. In the present thesis we 
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shall be concerned primarily with the nuclear quadrupole and 

hexadecapole interactions. The nuclear quadrupole interaction 

2 

is generally measured in terras of the coupling constant e qQ, 
vjhere eq is the total electric field gradient at the nuclear 
site due to the siarrounding charges and eQ is the nuclear 
quadrupole moment. 

Using a predetermined value of eo by means of some incle- 

2 

pendent experiment eq can be extracted from the e qQ data on a 
given system and the available Imov/ledge of the electronic 
structure of the system can he put to test. For this purpose 
the total field gradient is generally expressed as 


q = 


(1 - .) 


^-tlat 


+ (1 - R) 


'^val 


(1.5) 


where q 2 _^^ ^td give the field 'gradients due to the distant 

charges and valence electrons respectively and - 

-Ra represent the field gradients induced in the core elec- 
^val 

trons due to q. , and g , respectively. For the accurate 

"-L3.*C “V0-1 

calculation of q in solids therefore it is essential to have 
reliable theoretical values of Yoo and R,' Values of ^nd q^^2 

can be evaluated from the point-multi pole model and 'tine valence- 
electron wave functions. The quadrupolar polarizability ot^ is 
useful in the accurate calculations of based on point- 


ionic model, where the lattice potential V at the ionic site 
A is modified due to the induced quadrupole moment =ag 


V 



6 


The more fmdamental quantity Q, the nuclear quadrupole 
moment, is most conveniently derived from the quadrupole 
couplinq constants measured from the hrs studies on isolated 
atoms. The calculated field gradient q^^j_ in this case is 
corrected for the sternheimer effect according to 

q - (1-R) 

Thus the uncertainties in eQ and eq obtained from experiments 
to a large e:ctent depend on the reliability of (1 - ^oo ) and 
(1-R) values. 

The fourth-order gradient electrostatic poten- 

tial at the nucleus can be expressed as, 

n >16 = <1 >">16 1 *+ ™16 val 

where - n m, , . and -R,m, , represent the Sternheimer polarlza- 

“ 16 ion n lo vai 

tion effects in the core electrons arising from the distant 

ions and the valence electrons respectively. It is important 

3 

to note, here that values of (l - no, nre of the order of 10 

and can drastically change the fourth-order gradient, m^g 

due to distant ions. The m arises from the d and f valence 

16 val 

electrons. Using a predetermined value of H, the nuclear hexa- 
decapole moment, and the calculated value of (1 - ion' 

is possible to extract m. - , from the measured hexadecapole 

coupling constant e^'m^gll. It is then possible to obtain 
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experimentally information regarding d and f electron distri- 
butions in terras of an extremely sensitive fourth derivative 
terra 

Similar to the quadrupole moment case, the nuclear hexa- 
decapole moments derived from the experimental data on isolated 
atoms have to be corrected for the Sternheirner polarization 
effect (1 - • Thus, accurate calculations of (1 - ) 

and (1 - are highly desirable since experiraental data 

on hexadecapole coupling constants are now becoming available. 

The quadrupolar part of the static crystalline electric 
field at the valence site is shielded due to the Sternheirner 
polarization of the core electrons according to 

I 

A® = (1 - 03) (1.8) 

and it is necessary to interpret the experimental measurements 

T 

involving A° , e.g, in the study of tbe temperature dependence of 
Mossbauer quadrupole splittings# nuclear alignment experi- 
ments etc. The shielding factors, ^^1' external 

crystalline field at the core electron sites can lead to detec- 
table core level splittings in photoelectron spectra of heavier 
atoms. It is therefore essential to calculate and to 

interpret such data. 

The atomic quadrupole moment 0 ^ due to the aspherical 
np-electron density distribution can itself act as the source 



8 


of per tui .oation and induce a quadrupole raoment X'jithin the 
electronic shells in the atom, v/hich can change the value of 
0^ according to 

® - P-0 ) Op (1.9) 

Depending upon the atomic system under consideration can 
^ which case 'che ohserved atomic qua.drupole moment 
will have a sign opposite to that predicted Ijy the shielding 
uncorrected value of 0^. Therefore, it is of interest to 
calculate Rg for various atomdc systems. 

Although we shall ue dealing mainly with the quantities 
^03, R, ^ 2* ^co > 3.nd Rg in this thesis, the methods for calcula- 
tions of these and other polarizabilities and shieldings have 
much in comiaon and hence we shall outline the general theore- 
tical procedures that have been employed so far in the calcula- 
tion of polarizabilities and shielding— anti shielding factors. 

In the next section (l, 3 ) we give a review of various approximate 
methods available for the calculation of polarizabilities and 
shielding-anti shieldin'] factors. 

1,3 Approiciraa ce Methods for the Calculation of PolarizaJoilities 
and Shielding-Anti shielding Factors 

The various approximate methods available for the calcu- 
lation of pol ar izaoili ties and shielding— anti shielding factors 
can be broadly classified into the following categories: 
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(a) Statistical Methods [Gomhas (1944), Sternheimer (1950)]. 

(B) Coupled Ilartree-Foclc (CHF) Methods for Closed Shells 
[Dalgarno (1959), Kanelco (1959), Watson and Freeman 
(1963), Cohen and Roothaan (1964), and Lahiri and. 

Mulcher j i (1966)3, 

(C) Uncoupled Hartree(UCH) and Uncoupled Hartree-Foclc (UCHF ) 
Methods for Closed Shells [Sternheimer (l954). Das and 
Dersohn (1956), Dalgarno (1959), Kariplus and WoUcer 
(1963), Yosli-iraine and Hurst (1964), Langhoff, Karplus, 
and Hurst (196'G), Duff and Das (I960), Sadlej (1971 ) 
and others ]. 

(d) Perturljed Uncoupled Rartree-Fock (PUCHF) Methods 

Hirschf elder and co-v/orkers (1966), Kusher (1967), 
Epstein and Johnson (1967), Sadlej (1971, 1973), SteT-jart 
(1975) . 

(E) CHF Methods for Open Shell Systems [ Mukherj i (1975), 
Stev/art (1975)]. 

(F) Methods Including Electron Correlation Effects Explicitly 
[ Donath (1961), Kelly (1964), Das (1968), Nesbet (1970), 
Krauss (1972), Sinns (1973), Lindgren (1975)]. 

I. 3. A Statistical Methods 

Gombas^'' has developed a statistical procedure for 

calculating and has obtained satisfactory results for 
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heavy atoms, Sternheirner"^ has employed the Thornas-Fermi 
scheme to evaluate the electric field at the nucleus. Due to 
incorrect behaviour at large r this model is inadequate parti- 
cularly for polarizability calculations. In the case of the Thcmas- 
Fermi-Dirac model, Sternheirner found that the induced, field is 
far too large for licliter atoms. This method has been exten- 
sively used to ob'tain theoretical estimates of the total angular 
contributions to the multipole shielding-anti shielding factors. 


I.3.B Coupled Ilartree-Pocb. liethod for Closed- Shells 

We sliall first consider. the ’closed- shell’ IT-electron 
atomic systems v/hich can be satisfactorily represented by a 
single determinant (SD) vrave function v/ithin the Hartree-Fodc 
(HP) model- The coupled Hartree-Focl: (CHF) method is actually 
a generalization of the IIP idea which seems variationally the SD 

solution of lowest energy in a self-consistent manner either by 
using a perturbation expansion and obtaining the first order 
cor.rections to the unperturbed HF orbitals, or by directly 
solving the problem in presence of a finite perturbation without 
dsing the expansion. In the former method, Icnown as coupled 
Hartree-Poc]: perturbation (CHPP) method, the second order pro- 
perties can be calculated in two ways. One can solve the 
appropriate in'. omogeneo us differential equations or alterna- 
tively extremize the corresponding e no rgg^ functional. This 
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latter method will be called the coupled Hartree-Foclc perturba- 
tion variation (CHFPV) method. 

1 ^ 17 

Peng"^"^ and Alien‘d have discussed earlier the HF model 

for atomic system.s in presence of a static external field. 

1 8 

Langhoff/ Karjplus and Hurst (Ll'Ctl) have presented an instruc- 
tive analysis of the varioixs coupled and uncoupled methods in 
use. lie shall present here and in the next section the methods 
denoted as a, b, c, and d in the paper of LKH and identify 
various approximate methods in use to calculate multipole 
polarizabilities and shielding-anti shielding factors. In what 
follows, vre shall present a general outline of the CHP proce- 

5 

dure. 

Denoting the one-electron Foci; operator *and the external 
perturbatj.on as h(l') and h^(l) respectively, the problem reduces 
to solving the following coupled equations: 

{h(l) + X (l) “ J X^{l) =0, i = 1, 2, 3, . . .W (I.IO) 

where spin orbitals, X , obey the orthonormaliti^ criterion 


and h ( 1 ) 




H 


+ I < 

3=1 




.M 


P ) 
12 


X. > 

J 


Here P ^2 perrautation operator for electron 1 and 

Expanding x^ as 


( 1 , 11 ) 

( 1 . 12 ) 

2 . 
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(1.13) 

(1.14) 


and substj.tuting in eq (I.IO), the differential equations up to 


first order in X are obtained as. 


{ h°(l) -- } X° (1) = 0 

{ h°(l) - e°} xj (l) t-fh^d) ~ e^} (l) 



where h°(l) is 


the seroth-order Fock operator: 


(1.15) 


(1.16) 




(1-P12) 


X° > (1.17) 


Eq. (1.16) defines the coupled Hartree-Pock perturbation proce- 

190. ^1 

dure 'a' of LKK paper and has been used by Dulgarno, ' 

1 9*^0 *^0 18 

Kanelco, “ and Langlioff , Karpius and Hursr to calcvtlate 

3-nd Yoo for two-, three- and four-electron sequences 

and a^.-^ for he. Mg, and , Dalgarno and Kaneko have solved 

R 

the differential equation (I.ll) while Langhoff et al. have 
used the variational method based on product approximation for 
the perturbed function, viz., = f X°. 

Tlatson and Freeman^' have used a finite perturbation 
method within the unrestricted Hartree-Poclc (UliF) scheme and 



13 


23 


solved the equations corresponding to (1,5) to calcxilate 
quadrupole shielding-antishielding factors for heavier ions 
such as Cl”, Br”, l”. These authors have shown that within the 
UKF approximation, releocation of Kesbets sirametry and equiva- 
lence restriction^^ yields the angular and radial contributions 
to Yco (and R) respectively. Relaxing the sitnmetry condition, 
however, poses computational problems particularly for large 
systems and in the W method the angular estimates have to be 
based on less accurate uncoupled methods. Cohen and Roothaan 
have solved eg. (1.5) within the nartree-Foclc-Roothaan model 
v/ith the finite perturbation and have obtained i-or the 
members of two-, four- and ten-electron sequences. One of the 
drawbacks of the finite perturbation methods is that the HF 
calculations have to be carried through to a high degree of _ 
self-consistency in order to efficiently retain the effects of 
a small part of the Hamiltonian. This inherent error becomes 
magnified as the perturbation becoraes smaller and smaller. 

Lahiri and Mulcherji^^" have proposed a fully-coupled 
restricted Kartree..-Fock.-Roothaan perturbation-variation method 
and have calculated and values for He, Li, be. 

He and A isoelactronic sequences.. This is the most economic 
CHF method which gives Bo values close to the theoretical 
limit of K/Z. This indicates that the variational form used 
for the first-order wave function is sufxiciently rlexible. 



14 


Before closing this section on CKP methods the folloving 
reraarks on some of the procedures described here are in order. 

In the method 'a' of LKH the use of the product-approximation 
forces the perturbed wave functions to have the same nodes and 
a similar exponential decay as the unperturbed wave functions. 
Kanelco ' s procedure, due to its numerical character, accords a 
full variation to the perturbed wave function but it demands 
heavy computational effort. In Lahiri and Mul:herji's method the 
choice of an accurate variational function as a sum of Slater 
type orbitals (SID's) becomes complicated for systems beyond 
argon; hence calculations beyond argon series have not been 
attempted with this method. 

I.3.C Uncoupled Kartree (UCI-l) and Uncoupled Hartree-Fock (UCHF) 
Methods for Closed- shells 

The ’'various uncoupled metliods in use are the resifLt of compro- 
mises between accuracy and computational labour. With the assmp 
tion that the coupling terms in eq. (I.ll) do not contribute 
significantly -bo the total second order energy and therefore 
to second order properties such as polarizabilities and 
shielding factors, it has been the practice in uncoupled methods 
to completely neglect tliem. These methods can be best described 
as semi-empirical methods since they do not conform to any well 
defined physical model. 
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I . 3 . C (i ) Uncoupled Hartree Methods 

Our understanding of the polarization of closed shells 

is largely due to Sternheimer^~^' and the uncoupled 

perturbation-nuraerical method foinnulated by him. This method 

can be very briefly described as an uncoupled Hartree method 

(UCH) which uses Hartree or Hartree-Fock wave functions for 

the unperturbed state vritli-a local approxuifiation ■' to the 

effective potential in the first-order perturbed inhomogeneous 

Schrodinger equation. We shall discuss this method in greater 

25 

detail in Chapter II. Khi-ibchandani, Sharma and Das have 
shown that the inclusion of exchange among the perturbed vrave 
functions in the Sternheimer method significantly increases 
(in magnitude) for and Cl"” but values are not affected as 
much. This method rectifies the inconsistencies arising due to 
the' use of HF v/ave functions in the UCH method. Ray, Lee, Das 
and Sternheimer^ have extended the perturbation-numerical 
(differential equation) approach to incorporate intra- shell 
and inter- shell coupling effects which are ignored in the 
uncoupled founalism. 

27 

Das and Bersohn (DB) have formulated the variational 

193 . 

equivalent of Sternheimer ’ s method, Dalgarno has shown 
that for the lighter atoms and ions the small discrepancies 
betvreen the perturbation-n'umerical and variational methods can 
be ascribed to the different representations used by the two 
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authors for the unperturbed v/ave functions. Wilcner and 
29 

Burns have adopted the variational method of DB to obtain 
the multipole polarizabilities and shielding factors for 
heavier ions- We shall show in the subsequent chapters that 
for heavier ions the variational method underestimates the 
antishielding effect. The variational procedure of DE requires 
that the trial first-order perturbed wave functions are ortho- 
gonal to the associated unperturbed wave functions. The 

apparent difficulties concerning such an orthogonality require- 

30 

ment have been discussed by Ingalls. 

I.3.C(ii) Uncoupled Hartree-Pock Methods 
5 31 

Dalgarno ' has proposed the -uncoupled Hartree-Fock 

(UCHF) method in which both the unperturbed and perturbed v;ave- 

functions satisfy the Pauli principle. Dalgarno and McNamee^^^ 

have solved the appropriate differential equations -to calculate 

3 2 

o^/3ra / a and for Be. Yoshimine and Hurst, and Langhoff 
3 3 

and Hurst have used a variational version of Dalgarno ' s 

UCHF method^^ to calculate mul-tipole polarizabilities and 

shielding-anti shielding factors for a large number of atoms 

and ions. In this metlsod the energy functional is extremized 

without suEsjecting the trial perturbed wave functions to any 

34 

orthogonality constraints. Karplus and Folker have proposed 
a similar metliod in connection with the calculation of electric 
and magnetic suscep-tibilities of small molecules. There is. 
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however, an important difference between the methods of 

3 A 31 

Karplus and Follcer ^ and those of Dalgarno which has been 

X 3 

first iDointed out by LKH, In the following vie shall present 
LKH methods b, c and d as obtained from method a (see sec. 
I.2.B(ii)) and bring out this difference. 


Cancelling the self -potential term (j=i) in h°(l) in 
ecf. (1.17) v/ith (j=i) term of the explicit sum in eq. (1,16) 
and neglecting all the terms involving first order functions 
Xj, other than the one being considered, equations are obtained 
for LKFI^® UCHF method 'b'; 


{h°(l) - £°} xjci) + {h^(l') - X°(l) = 0 (I.18-) 

h°(l) is now a modified Foclc operator in tine sense that under 
independent variation of spin orbitals only x° is an eigen- 
function of h?. It should be noted that in the closed-shell 
version of eq. (1.18), x^ is not an eigenfunction of the 
eeroth-order Foclc type operator v/hich replaces h*^. Owing to 
the cancellation of the self-interaction term, the perturbation 
in method 'b' is obtained in presence of N-1 electrons only. 

neglecting completely the simmation terms in eq. {T.16} 
the equations corresponding to Lr-IH method 'c’ are obtained as, 

{h° (1) - } xj(l) +{h^(l) -e^}X°(l) =0 (1.19) 

where h'^ is the zeroth-order Foclc operator. Eq. (1.19) , is 



UCHF method. It can be seen that 


5 31 

equivalent to Dalgarno ' s ' 
due to the self-interaction term in h'^ the perturbed vzave 
functions are obtained in presence of an inappropriate core 
potential due to H electrons. This method provides an upper 
bound to the CHP second order energy. 


Substituting for h°(l), the Ternlcin'^ operator 
eq, (I. IS) or (1.19) LIT! have defined their UCHF method d as. 


{g°(l) - e°> x5;(l) + th^(l) - e^} xJCl) 


where g° ( 1 ) 
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The varietionel method of Karplus and Kbllcer is similar to 
LKH method ’d’. Therefore in the Karplus-I^lher method the 
perturbation is obtained in presence of (N-l) electrons as 
against the core potential due to N electrons in Dalgarno ’ s 
UCHF method. Duff and Das have formulated the comterparts 
of LKH methods 'a' and 'b' wherein the orthogonality restriction 
on perturbed one-electron wave function is relaxed. The counter- 
part of method b has a computational advantage over the original 
method due to simplifications’ achieved in the ortho normalization 
process. Sadlej has modified Karplus-'holker method by 
including intra-orbital contribution of the non-local part of 
the HP potential and self-consistencv; respectively. His 
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calculations of two-, four- and ten-electron systems 

are in good agreement with the CHF results. In a subsequent 
37 b 

paper Sadlej ' has omitted the nonorthogonality contributions 

and has obtained improved results for a-, and for ten-electron 

d 

atomic systems. 

The more approximate uncoupled Kartree perturbation 

3 8 

variation method of Pople and Schofield has been applied to 

3 9 

calculate, for all neutral atoms. Adelman and Szabo have 

described the Coulomb-approximation method to calculate a ^ 

and 6^ and have obtained good agreement v/ith the experimental 

values for monovalent ground and excited S-state atoms. 

The Bcj values differ significantly from N/Z. Ahlberg and 

Goscinslci^*^ have calculated a ^ and using a local linear- 

response method for singlet ground state atoms up to Z = 3S 

and obtained excellent agreement with the experimental and 

theoretical 3co values. These- authors have observed that the 

41 

full Slater exchange rather than the optimized counterpart 
yields more accurate results. 

I.3.D Perturbed Uncoupled Hartree-Foclc (PUCHF) Methods for 
Closed- shells 

Tuan, Epstein and Hirschfelder"^^ (TEH) have shown that 
starting with Dalgarno UCHF wave functions at the zeroth-order, 
the first order correlation correction to the UCHF second order 
energy'' can be obtained without solving any new equations. 
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According to these autliors/ the first order correction 
to the second order property can be obtained from 
K N 
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and the superscripts o and 1 refer to the unperturbed and the 
first order perturbed ^'j-ave functions respectively, 

43 

Musher has given a more rigorous treatment for such 

corrections and has shoi^rn that the UCHF result is equivalent to 

the first term in the double-perturbation theory expansion and 

the CHF result differs from first tv;o terms of the general 

expansion only by higher order terms. He has also shown that 

when per tixrbat ion theory is applied to improve upon CHF method 

the resulting solution offers no additional advantage. 

44 

Epstein and Johnson have carried out such perturbed uncoupled 
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Hartree-Foc]c (PUCHF) calculations of a-, and ct for H, Li, Be 
iso electronic series and have found that except for the near- 
neutral atomic systems these results are in satisfactory agree- 
ment with the CKF results. In a diagrammatic double-perturbation 

45 

analysis Caves and Karplus have shown that the first-order 
correlation correction to the second order UCHF energy coincides 
v/ith the first iterative solution in the CHF theory. Schulman 
and Musher"^® improved the situation by proposing a geometrical 
approximation on empirical grounds. These authors have shown 
that -the hydrogen atom dipole polarizability can be obtained 
from the double perturbation expansion based on a Hartree-Foclc 
potential, FT, by approximating the perturbation series by a 
geometrical series according ' to. 


geom _ o , , 

“cl " ad ^ ^ " 

vjhere corresponds to H 
d o 

tion correction to a°. 


(aj/a°)} ^ (1.24) 

and is the first-order correla- 
cl 


>mos, Tuan, Sadlej and Jaszunski^^ and Jamieson and 
5 0 

Gafarian have independently given theoretical justifications 

51 

for the GA method. Tuan and Davidz have performed PUCHF 
calculations of and for two-, three- and four-electron 
systems v/ith varying amount of self- interaction terms in the 
effective Hamiltonian and have shown that an optimum self- 
interaction term along vrith the GA can be a good substitute 

52 

for the CHF calculations of the second-order properties. Sadlej’ 
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has generalized the PUCHF method of TEH for the approximate 

53 

SCF unperturbed orbitals. Stewart has applied GA to Dalgarno 
UCHF by using the uncoupled differential equation (rather than 
using the appropriate functional method as in Saldej 's GA to 
Karplus-Itolker method) and obtained a-, and a in excellent 

O. 

agreement with CHF results. In the PUCHF calculations of 

polarizability t.he computer time required is one-fifth of . that 

required for CHF calculations. The first-order correlation 

corrections to the shielding factors which are bilinear in the 

perturbation pararaeter require twice the computer effort as 

54 

compared to the polarizability calculation. Litt has 
introduced an effective- field method which is equivalent to GA 
for the polarizability calculations if the zeroth-order appro- 
ximation is chosen to be Sternheimer approximation. The GA 
methods of Sadlej and Stevens seem to be worth pursuing for 
heavier systems. 

I,3.E Coupled HP Methods for Open- Shells 

55 

Mukher j i and. coworkers have extended the CHF perturbation- 

variation method to calculate multipole polarizabilities and 

shielding factors for open-shell atoms. This method neglects 

the core rearrangement and core x^olarization effects but 

includes all the first-order correlation and parts of higher 

order correlation effects. The core rearrangement effects 

55b 56 

have been recently included by these workers. Stewart 
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has solved the numerical version of CHFPV method of Mukherji 
and coworhers and has obtained good agreement with the latter 

4, o 

calculations and the experimental result for of A (“^p). 

I.3.F Methods Including Electron Correlation Effects 
Explicitly 

For accurate calculations of second-order properties it 
becomes necessary to go beyond Sternheimer polarization and 
include two-particle interaction in an explicit manner. 

57 

Donath has adopted a variational approach including 

configuration interaction to calculate for Me. Similar 
5 8 

calculations of other first row atoms have been performed. 

In the case of larger number of valence electrons a balanced 
treatment of the unperturbed and perturbed wave fimctions 
cannot be guaranteed since the multiconf igurational structure 
has to be ascribed to the perturbed v/ave function in the lowest 
approximation. Indeed such methods have produced low values 
for the dipole polarizabilities of G, N and Ne. For the 
cases with smaller number of valence electrons the variational 
configuration interaction method is certainly the most powerful 
one. 

59 

Sitter and Hurst have used CHF method using an extended 
polarization basis to compute of Ne. Stevens, Billingsley 
and Krauss^*^ used a coupled multiconfiguration self-consistent 
field (QiGSCF) method to calc u). ate a ^ for the first row atoms 
to within +5% in acouiracy. 
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Meyer and coworliers^ have performed the most accurate 

calculations of for the first row atoms ('v2% Accuracy) using 

pair natural-orbital configuration- inter action (PNO-Cl) and 

coupled electron pair approximation with pair natural-orbitals 

62 

(CEPA-Pt-TO ) . Ahlrichs and cov/orlcers have reported a computa- 
tionally improved version of these methods to calculate 
correlation energy in Me and smaller molecular systems. These 
methods have not been applied so far to the calculation of 
shielding factors. 

63 

Kelly first applied the linlced-cluster many -body pertur- 
bation theory (LQ-lBPT,) expansion method of Brueclcner and 
Goldstone^^ to atomic systems and obtained Be » “g and 

for Be including correlation contributions up to second order. 

65 

Similar MBPT calculations of Li, C, 0, Me have also been 

6 6 

carried out. Ebran calculated multipole polarizabilities 

67 

and shielding factors for Me, A, Kr and Xe. Ray, Lee and Das 

have calculated for Fe^"^, Rb"*" and Cs"^ using LQIBPT method. 

Similar calculation of R has been first reported by Lyons, Pu 

and Das for Li(^P). These authors have rxjted that the effect 

of second-order correlation correction to R is a small positive 

contribution with respect to the first-order correlation correc- 

69 

tion. Accurate calculations of Mesbet using variational 
Bethe-Go Id stone approach have suggested that the two-particle 
correlation corrections give rise to near cancellation of the 
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7 O 

single-particle correlation effect. Hameed and Foley have 
applied the perturbation method to second order within RHF 
scheme and obtained an agreement v/ith Nesbets' calculations 
to within 2 %, This is probably due to the limited nOTber of 
configurations v/hich Hameed and Foley used to construct the 
first-order wave functions. Garpman, Lindgren, Lindgren and 

7 1 

Morrison have formulated an elegant effective-operator form 
of LCMBPT method and applied it to calculate ^ init io, the 
hyperfine interaction in L.i(np) and Na(np) sequences- Instead 
of generating a complete set of excited states as in Kelly's 
approach the first- and second -order corrections to the HF 
value are obtained by solving inhomogeneous one- and two-particle 

9 

equations. Their value of R( P Li) is in excellent agreement 

7 2 

with Nesbets' calculation. Das and coworkers have recently 

shown that the discrepancy in their earlier calculation on 
2 

Li(‘"p) is a result of neglecting some diagrams corresponding 

to tvro-particle (ls-2p) correlation. The calculations of R 

73 

with similar accuracy have been reported for the first row 

2 + 7 2 

atoms, allcali metals and Fe ion. These calculations are 

in reasonably good agreement with the corresponding calcula- 

7 4 

tions based on Sternheimer method. One remarkable advantage of 
the Brueckner-Go Id stone approach lies in its transparency 
which makes it suitable to physically evaluate^^^ the numerical 
results at various levels of approximation. The convergence 
of the perturbation series seems to depend considerably on the 



26 


particTolar system aind this leads to practical difficulties in 
carrying out all the calculations with the same accuracy. The 
electron pair appro ache on the other hand, include 

three-particle correlation effects which become important 
beyond nitrogen. Such calculations are forraidable even with 
present day computers. 

It is evident from the survey of theoretical methods 
presented alsove that not many detailed calculations of "^ca , 

R/ 02 / ’^co ^ have been carried out using reasonably 

accurate uncoupled methods. There is, therefore, a need to 
undertalce a project involving extensive calculations of these 
Sternheimer shielding-anti shielding factors in a consistent 
manner. Such calculations will be equally Important for both 
experimentalists and theoreticians. An attempt has been made 
in this thesis -to malce some contribution in this direction. 

In the next section (1.4) we shall review the experi- 
mental evidences that support the viewpoint that the shielding- 
antishielding phenomena play an important role in the study 
of quadrupolar and hexadecapolar interactions. 

1.4 Experimental Evidences for Sternheimer Factors 

As mentioned earlier in this thesis, we are interested 
in the quantities Y„ , R, Ho, , and Rq . So we shall confine 

our attention in this section to the quadrupole and hexadecapole 


interactions. 
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In the earlier literature, the importance of Sternheiraer 

polarization effect corresponding to the electric quadrupole 

interaction has teen recognised mainly through nuclear magne- 

7 5 

tic resonance (lI'iR) and nuclear quadrupole resonance (NQR) . 
Later, other methods such as Mossbauer effect, time differen- 
tial perturbed angular correlation (TDPAC), low temperature 
colorimetry and atomic hyperfine structure studies have also 
emphasized the role of Sternheimer effects. In the following 
we shall present a survey of these experimental evidences for 
Sternheimer polarizations. 

I. 4. A Quadrupolar Interactions and (1 - Y„ ) 

Sternheimer”^^ estimated the effective antishielding 
factors for ®^Br, ^^Cl, ^^Cl and in the polar molecules 

of KBr, I'laCl, KCl, and KI, respectively from the available 
quadrupole coupling constants derived from molecular beam 
electric resonance studies'^ and the known‘s Q values and 
concluded for the first time that the experimental results give 
support to the calculated antishielding factors. 

I.4.A(i) Magnetic resonance studies 
7 8 

Pound has shown that the perturbing nuclear quadrupole 

interaction in nonccbic ionic solids can split and/or shift 

7q 7 23 27 

the high field ' MMR line in the cases of Li, Na, and Al. 

The spin lattice relaxation time (T^) of these lines were also 
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803, 

shown to be affected by such interactions. Bersohn and 

80]o c 

Volkoff ' have given perturbation expressions to extract the 

quadrupole coupling tensor components from the NlMR spectrum of 

noncijbic single crystals,. Several theoretical models such as 

81 

acoustic ionic -model of Van Kranendonlc, charge- transfer 

82 a 

covalency model of Yoshida and Moriya and overlap-covalency 
model of Kbndo and Yamashita have been proposed to explain 
the quadrupolar relaxation time T^. 

Q 1 

Van Kranendonk has estimated the Combined effect of 

o 1 27 

CO valency and (1 - y„) to be >10 from I NMR relaxation time 

83 

measurements in Id. Kawamura, OtsuJca and Ishiwatari have 

observed striking reduction in intensity of the satellite lines 
23 

of Na in llaCl crystal when small amounts of NaBr are added. 

In such mixed-crystal experiments the field-gradient results 
from the size differences in solute and solvent ions. From 
the second-order shift of the central line at 10% of solute 
concentration, kawamura et al. have estimated (1 - Ya, ) 'v 9 for 

• 4 - 84 

Na . Otsuka and Kawaraura have sirailarly studied the second- 

79 81 

order shift of Br and Br resonance lines in KBr-NaBr system 

— otr 

and concluded that (l v 40 for Br . Bersohn has inter- 

preted the quadrupole coupling constants available from NMR 
data on I'Ta, 'Al, '^Cu and Cu nuclei in ionic solids using 
the point-charge model and has found that the abandonment of 
the factor (l - ) representing , the polarization effect removes 

all agreement with experiment. 
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Wikner and Das have shown that the ratio of field gra- 
dients obtained by Kiddle and Proctor^^ for ^'^Rb and ^^^Cs NMR 
in tutton salts can be explained by considering the theoretical 

I , O'! 

ratio of (1 -Y ) for Rb and Cs"^ ions. Valiev has obtained 

CO 

2 

reasonable results for (1 - Y^ ) frora the 1-MR line-width studies 

3 "4“ 3 “f” 88 

on A1 and Ga in solution.- Wikner, Blumberg, and Hahn 

have included the effects of induced dipole moments due to 

81 

optical modes into Van Kranedonlc model and have assumed the 
theoretical value of (1 - ) for Na"*" to estimate (1 -Y^, } for 

Cl" as 20 from the ratio of spin-lattice relaxation times for 
the two nuclei. From the temperature dependence of the quadru- 
pole coupling for Al^'*' in cCnh^)^ Al (so^) 2 • and Al^"*" and 

Ga in isomorphous ionic solids. Burns has obtained good 

agreement with experiment on including theoretical (1 ) 

3 4- 3 + 29c 

for Al and Ga . Burns and Wilcner have concluded that the 

experimentally estimated values of (l --Yco) positive ions 

are in reasonably good agreement with the theoretical values 

but the empirical values for negative ions are 3-5 times 

90 

smaller than the calculated values. Das and Dick have criti- 
cized these estiraates for the negative ions on the basis that 
the empirical values for the negative ions have been derived 
in a crude fashion. Thus, the empirical estimate of Ya, for 
Br by Otsuka and Kawamura has been derived from a simple 
Compression theory which leads to erroneous estimates of the 
the displ aeon ents of the lattice points and also neglects the 
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polarization contribution of the nearest-neighbour ions which 

90 

are important. Das and Dic]c have included these contributions 
in interpreting the results on NaCl-Br, NaBr-Cl and KBr-Na and 
have observed that the NI-IR data can be satisfactorily explained 
on the basis of theoretically calculated values of Yoo for Na"^ 
and Br~ . 

91 

Wilson and Blume have measirred the quadrupolar shifts 
of NMR lines of the off-centre Li"*" and Cu"^ impurity ions in KCl 
and have noted that the ratios of first- and second-order effects 
obtained in the tv7o cases can be explained by assiwing the theo- 

QO 

retical values of (1 -Y„ ). Raymond has used the point- 
multipole model and obtained the best fit value of Yeo -4.9 

to explain 34 independent field gradient components extracted 

27 93 

from A1 RMR in three Al^SiO^ polymorphs. Sharma has used 

the overlap model along with the Sternheimer shielding-anti- 
shielding factors and obtained the nuclear quadrupole moment 

for ^'^Al from the NMR data”^® on Al^O^. His value of Q(^'^Al) = 

94 

0,148 b is in excellent agreement with the atomic beam result 
of 0.149b. 

There arc innumerable examples from pxare quadrupole 

95 

resonance experiments vj-hich have led to the determination 
of nuclear quadrupole interactions in solids. These have also 
provided evidences in favour of shielding-anti shielding effects. 
Hewitt®® has estimated (l«-Y^ ) = 16 for Nb®"^ from the NQR data 
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93 

on Nb in which compares well with the theoretical esti- 

97 

mate of 15. Barnes, Segal and Jones have shovm that using the 

empirically estimated values of for Cl”, Br” and l” as -10, 

-35 and -45 respectively, the NQR data on a series of solid layer- 

type metal halides can be e^qslained if a small increase in the 

value of s -hybridization due to covalency is assumed. Rao and 

983 . 

Narasimhamurthy using a point charge model with induced dipole- 

moments found thatYco = -7 5 for Cl” is necessary to fit the NQR 

99 

data on CUCI 2 . 2 H 2 O. Narath has used a point-charge model and 

observed that Ym = -27 for Cl” gives a reasonable agreement with 

the NQR data on C 0 CI 2 . 2 H 2 O. Tlie value of (1 - ^ca ^ estimated by 

Rao and Narasimhamurthy for Cl” appears to be too high. In fact, 

a later calculation on CUCI 2 . 2 H 2 O including charges and effective 

dipole moments by Raj and Amirthalingam^®^ favours a value of 

Yoo = -27.4 for Cl”. Liu^^*^ has reported the best fit value .of 

ionicity and Y^ for Cl” as 0.911 and -31 respectively to explain 

the NQR and NMR data on C 0 CI 2 . 2 H 2 O. implying a soft sphere model 

35 

to explain the experimental Cl NQR data on GdCl^ and a series 
of alkali-metal-hexachlorostanates. Current^ has estimated 
Y 'u -30 for Cl”. 

CO 

121 123 

The quadrupole coupling constant for Sb and Sb in 

102 

antimony metal have been measured by means of NQR and low 
temperature calorimet3:Y^*^^ methods. Krusius and Picket^ have 
obtained the extra-ionic contribution to the field gradient by 
subtracting the ionic field gradient corrected for anti shielding 
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and assuming^^^ Q(^^^Sb) = 0.5 b. Their value of 560 x 10^^ 

2 

V/cm compares well with the conduction electron contribution 
calculated by Hygh and Das.^^^ 

The quadrupolar relaxation time in liquid metals and 

1 07 

III-V alloys have been measured. From such experiments 

Haldar‘°^ has estimated values for Ga^"^ and In^''' close to 

-10.5 and -24.9 respectively. He has observed that the theore- 

2"h T 09 

tical value for Hg calculated by Rao and Murthy is boo low 

by a factor of 2 or more. Kerlin and Clarh^^*^ have found that 

in liquid GaSb the alloying can significantly enhance (1 - Yo, ) 

value for Sb over Ga due to an increase in negative charge 

on Sb. 

From electron nuclear double resonance (EtTDOR') studies on 
F-centres in KCl crystal, Feuchtwan^^^ has estimated that the 
theoretical (1 - Y^o ) for Cl is probably too large by a factor 
of 2 or more. Zimmerman and Valentin have made EMDOR studies 

9 4 - 

on uniaxial stress- induced quadrupole splitting for Eu doped 

in BaF^, SiF^, ^^^2 Using the theoretical values of 

2 4- 

Ye, for Eu these authors have obtained reasonably good agree- 
ment with the experiment. 


I.4.A(ii) Electric Field Gradients and Acoustic Effects 
in Magnetic Resonance 

Kaslter^^^ and Al'tshuler^^"^ have independently predicted 
the possibility of observing saturation effects in NMR by inducing 
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the fluctuating field gradient by means of ultrasonics. The 

coupling between small elastic strain and electric field gradient 

115 

can be represented by a fourth-order tensor and in Voigt 
notation two independent tensor components ^44 

sufficient for the calculation of quadrupolar transition proba- 
bilities for a cubic crystal. In one method^ the S components 

can be determined by studying the quadrupolar broadening of NMR 

115 

lines due to an applied static uniaxial stress or due to the 

presence of internal strains from a hnow'n concentration of 

dislocation^^^ or impurity atoms^^'^ in a crystal. In the second 
118 

method ultrasound of known distribution of elastic strain 
amplitudes at Larmor or twice the Larmor frequency is applied 
and the transition between the nuclear spin levels induced by 
the time-varying electric field gradients at nuclear sites are 
studied from the decrease in intensity (saturation) of NMR 
absorption in a pulse or cw NMR experiment. In the more accu- 
rate method of nuclear acoustic resonance (NAR) , the resonant 
acoustic absorption by a spin system is directly measured as an 
additional attenuation. In NAR the elastic strain magnitudes 
do not enter into the determination of S components for certain 
directions of acoustic wave propagation and polarization in a 
cubic crystal. The accuracy of these three methods is estimated 
to be +15%, a factor of 2, and + 6 % respectively. 

Otsuka^^^ has estimated (1 - ) 'V 50 for I~ from reduction 

1 27 

in intensity of I resonance due to dislocations in KI. His 
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estimate is however based on the choice of (l - ) -u 40 for 

- 117 

Br . Fukai has expressed the gradient elastic tensors for 

^"^Na, ^^Cl, ^^Br and in a series of allcali halide solid- 

solutions in terms of antishielding factor, nearest-neighbour 

overlap integrals and degree of co valency. He has found that 

the theoretical for the ions lead to 2/^44 Na'*‘, Br” 

and I in fair agreement with the experimental data. A better 

agreement results by assuming a reduction of 10%, 40%, 50% in the 

theoretical Ym values. Due to the lac]i of availability of good 

atomic wave functions the overlap integrals were calc-ulated in 

1 20 

an approximate manner. Marsh and Cassabella have included 
second nearest neighbour overlap effects and obtained satis- 
factory agreement with the measured S components in uniaxially 
strained sodium halides using free ionic Y^ for Na'*’ and 

contracted for Cl“ and Br~. The nonlocal terms in the 

117 

overlap effects were neglected by Fulcai and Marsh and 
1 20 121 

Casabella. Ilcenbery and Das have calculated the electric 

field gradient by including these terms in the overlap model 
and have observed that an increase in Y^ for positive ions 
and a substantial decrease in Y^ for negative ions would give 
a tolerable agreement with the experiniental results without 
invoking any charge- transfer covalency. 

The predictions of Kastler^^^ and A1 ' tshuler^^^ regarding 
spin-phonon coupling has been realized in a pulsed NMR ultrasonic- 

saturation experiment on NaClO^ crystal by Procter and cowor- 

122 123 

leers. Jennings, Tantilla and Krauss have compared the 
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relative saturation behaviour of I and Na resonances. 

Assuming the theoretical value of 7^ for Na'*' they have estimated 

1 1 R 

'[j 50 for I . Taylor and Bio cmbergen"^ have measured 

S components in Nad single crystal by ultrasonic- saturation 
method and estimated (1 ) 'v- 5 and 9 for Na"*” and Cl"’ respec- 

tively. This experimental value for (1 - Yc ') for Cl” appears 

12 4 

to be very lov^. Antotolslcii, Sarnatskii and Shutilov have 

7 

studied Li in LiF by ultrasonic-saturation and derived 

+ 125 

Yco = 0.5+0, 3 for Li . Anderson and Karra in a similar study 

have derived 4. 4 + 0 . 4 for Li’*'. The reason for this discre- 

pancy in the two measurements is not clear to us. Kanashiro 

has explained the experimental values of S components in alkali 

bromides by including overlap effects and assuming the free 

1 27 

ionic Y„ values. Kanashiro, Ohno and Satoh have studied 

23 

ferroelectric NaN02 by ultrasonic-saturation of Na resonance 
and found that a point multi pole model coupled with the theore- 
tical Y^ for Na"*" gives reasonable agreement with the experiments. 

Bolef and Menes have studied "^Br and I in KBr and 

KI single crystals respectively by NAR and estimated (1 - Y^ ) 

for Br and l” to bo 26 and 38 respectively. Sundfors has 

carried out sirailar measurements on III-V semiconductors and 

1 29 

S'undfors and Tsui have analysed their data on to derive 

(1 - Y„) for Al, Ga, In, As and Sb as 16, 45-59, 95-108, 79-82 
and 85-104 respectively. Sundfors, Wang, Bolef, Fedders and 
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and V7estlake^^° have studied NAR in a hydrogen-free pure 

Ta single crystal and obtained the effective (l - Yo, ) to bo 
10 times the theoretical value. 

I.4.A(iii) Electric I'ield Gradients from M‘dssbauer and TDPAC 
Studies 

The quadrupole splittings obtained from Mossbauer effect 
studies on ^'^’^Pe. in solids have been generally interpreted on 
the basis of eg. (l.S). The agreements obtained via eq. (1.5) 
with the experimental data can be regarded as indirect eviden- 
ces in favour of Sternheimer shielding-antishielding .effects. 

,131 

The Mdssbauer data on ferric compounds have been interpreted 

on the basis of ionic model. In the case of ferrous compounds,^' 

on the other hand, the valence contribution is dominating. 

Such interpretations have led to a range^^^ of values (0.1-0. 4b^ 

of nuclear quadrupole moment Q for Fe. Including the effects 

due to overlap and Sternheimer shielding-antishielding factors, 

Sharma^^^ has resolved this ferric-ferrous anomaly to obtain 

Q(57mpe) _ o,18 + 0.02 b and has satisfactorily explained the 

Mossbauer data on several rare earth garnets. 

1 R 3 "I" 3 4“ 

Bhide and Hegde have compared the efg at Sc and Fe 

44 

sites, as determined by Y-y directional correlation of Sc 

and ^'^^Ve Mossbauer effect studies in PbTiO^, respectively. 

Assuming the theoretical value of for Fe^'' these authors 

3 136 

have estimated Y^ = -19 for Sc Sharma 


has incorporated 
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the effects of ligand distortions due to internal crystalline 

electric field in the overlap model and analysed the time- 

differential perturbed angular-correlation (TDPAC) data on 

and ^^^Cd CdCl^. His analysis gives Q(^^^Cd) = 0.89+0, 05b 

which is very close to the value of 0. 9 b obtained from the 

138 111 

circular polarization measurements of 247-lceV Y-ray of Cd, 

2 

In their interpretation of the e qQ data on metals 

137 139 

obtained by TDPAC and other methods Das and coworhers 

have calculated the conduction electron and distant charges 

contribution to the field gradient in Be, Mg, Zn and Cd metals. 

» 

The nuclear quadrupole moments obtained from their analysis 
are consistent with the nuclear structure systematics and the 
calculated field gradients are in quantitative agreanent with 
the experimental values for Be and Mg. 

140 

Raghvan, Kaufmann and Raghvan have shown that the major 
portion of extra- ionic field gradient in metals and dilute ailloys 
is universally correlated to the antishielding-corrected' ionic 
counterpart and is of the sign opposite to it. This correlation 
is the strongest evidence in favour of the existence of the 
term (l - Y^^). 

A quantitative agreement between the estimates of (1 - Ye,) 
obtained on the basis of a point-charge model in solids and the 
theoretical values cannot always be expected since several other 
important electronic effects might have been ignored. The best 



38 


way is to use the calculated lattice contribution to the field 
gradient corrected for the theoretical (1 - Yoo ) and interpret 
the discrepancy in terms of charge- transfer and overlap- 
covalency effects in the case of molecular solids or conduction 
electron contribution in the case of metals and alloys. If an 
agreement is obtained with the point charge model in the case 
of molecular solids, it is possible that the electronic effects 
are not important. However, this should be substantiated by 
independent experimental observations. 


I.4.B Quadrupolar Interactions and (1-R) 


The analysis^*^^ of hyperfine structure measurements on 
isolated atoms has produced striking evidences in favour of 
the Sternheimer shielding- anti shielding factors. In order to 
evaluate Q from such data the field gradient, or equivalently 
the value of <r >^, is calc\xLated theoretically or deduced from, 
other experimental quantities. The standard procedure has been 

_3 

to estimate a value of <r > from the magnetic hyperfine struc- 

_3 

ture or from fine structure. Using this < r > value the uncor- 
rected experimental quadrupole moment is obtained. In the 
case of np excited alkali atoms, for example. 


1.06395 X 10 


-2 


^un corrected 


-3 

< r > 


^3/2 barns) (1,25) 


np 


wh 


ere the quadrupole interaction constant ^ 2/2 in MHz 


141 
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_3 

and <r in a.u. The Sternheimer correction is then applied 

to obtain the corrected spectroscopic Q as. 


Q ~ 0 / (l— R) 

^corrected ~ ^uncorrected ' ~ 


(I. 26) 


Use of appropriate Sternheimer correction factors have been 
helpful in explaining several anomalous ejq^crimental resiilts. 

It has been observed"^ from hfs studies on ^Cu in 

929 9 10 

3d 4s ( ^ 2^2 ^5/2^ excited states that the 

7 4a 

estimated Q values differ by more than 40%. Sternheimer has 

shown that this difference can be quantitatively explained if 

the Sternheimer correction corresponding to 3d and 4p valence 

143 

electrons are included in the interpretation. Murakawa has 

studied 5d and/or 6p states of ^^^La, ^"^^Lu and and 

estimated R_(5d) = -0.4 and R_(6p) = -0.1. Sternheimer"^ has 
Q Q 

calculated R^CSd) =-0.4 for La v/hich compares well with the 

1 44 

experimental estimate. Childs has measured the magnetic 

dipole and quadrupole coupling constants for 4f and 5d valence 

159 8 2 

electrons from the hyperfine structure of Tb in 4f 5d6s 

65 7 4b 

state. Similar to the case in Cu, Sternheimer has shown 
that the experimental ratio of the coupling constants can be 
quantitatively attributed to the differential Sternheimer 
shielding for 5d and 4f electrons. The experimental studies 

on several np excited state alkali metal atoms have also 

145 1451 

shown similar evidences. Belin and Svanberg have shown 
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that the quadrupole moments derived from the ejqjerimental 

2 

quadrupole coupling factors for np P 3/2 ~ 5,6,7, 8 ) 

87 

sequence of Rh shows a decreasing trend. The application 
of Sternheimer correction factor, R (np), in the sequence 
leads to a more consistent value of 0.13 b for all the cases 
studied. 


Svanbcrg^'^^^'' has similarly found that Q(^^K.) corrected 

2 

for Sternheimer shielding factor in np P 3/2 
sequence leads to a more consistent value of 0.059 b as 
against the uncorrected values of 0.82, 0.85 and 0.86 respec- 
tively. From a comparison ?oetween Q values derived from 

hyperfine structure^"^^ of muonic ^^Nb with the hfs measure- 

4 32 ,147 

ments on 4d 5s and. 4d 5 s states, Buttgenbach and Dicke 

have estimated -0.123. 


In the calculations of nuclear quadripole moments from 

the atomic hfs data according to the semiempirical method, 

(e.g., in ^"^Cu, ^"^Rb, and ^^Nb) no correction is applied 

71 

to the magnetic dipole interaction (a/p). Lindgren has 

shovm that such corrections are important for nd excited 

—3 

states of Rb. VJhen < r > is determined ^ i^.i.tio, the L=0 
terms should be included in the calculations of R. In the 
semi-empirical procedure, hov/ever, these contributions are 
already contained in a/y . Correlation effects should be also 
considered since they need not be the same for magnetic dipole 
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and electric quadrupole interactions. Thus/ the apparent 

agreement between the quadrupole moments of Li(2p) obtained 

from Sternheimer semi- empirical method and the calculation of 

Lindgren and coworlcers is due to the accidental cancellation 

of .polarization effect on the magnetic interaction with the 

23 3 

net correlation effect. In the case of Na( P) these 
secondary effects are-negligible and one expects the same to 
hold good for np states of heavier allcali metal atoms. 

The Mo'ssbauer experiments on the rare earth compounds 
have predicted R = 0.1-0. 3 for tripositive rare earth ions 
vrhich compares well with the theoretical values of ‘'^O.iS 
obtained using Sternheimer method. 

1 49 

Dunlap, Kalvius and Shenoy have presented a clear 

evidence of R for 5f electron from their Mossbauer quadrupole 

splitting studies on uranyl and neptunyl compounds. These 

authors have assumed that the difference in auadrupole split- 

—3 

ting in the tvso cases should correspond to (1 - <^~ > 5 f 

where < r is derived from tlie relativistic HF wave functions 

(3 "f" 

for Np . Their estimated R__,. value is +0.32. On the basis 

5f 

of the known values for R this value appears to be high. 

I.4.C Hexadecapole Interactions 

In a few recent experiments the electric hexadecapole ■ 

iSOs.— c 

interaction has been measured in solids and in the cases 
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of free atoms ' namely Ho and Es. Sternheimer has 

calculated and ri^ and found that the hexadecapole polariza- 
tion is more significant than the quadrupole case. The hexa-- 
decapole moment obtained for Ho has been corrected for the 
hexadecapole shielding effect (1 - R^) . Theoretical calcula- 
tions of too and R^ are however available in a very few cases. 


I.4.D Shielding Factors (1 - o ) and (l - ^ to the 
Crystalline Electric Field , ^ 


The evidence for the shielding factor (l - *^ 2 "^ 

static crystalline electric field at the valence electron site 

1 - Y 

148 

has been obtained as a ratio of for several effect 

^ " 2 

nuclear alignment and a combination"^ of electron paramag- 
netic resonance and optical spectroscopy experiments. These 

values are in reasonably good agreement with the theoretically 

153 

calculated ratio. Pelzl has analysed the available hyper- 
fine interaction data on rare earth metals and concluded that 
1 - 'i^eo / ‘^2 with respect to the free ionic values. 


Novalcov and Hollander^^"^ have measured the splittings of 

%/2 core levels caused by internal electric field 

gradients in the compounds of Th, U, Pu and Au by means of 

155 

photoelectron-spectroscopy. Gupta and Sen have shown that 
the measured splittings in “^^3/2 levels can be 

qualitatively explained as due to the differential Sternheimer 
shielding (1 - 4p and 5p sites respectively. Their 
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estimate of the ratio of splitting is 1.7 for Au (5d^^6s^) 
and 1,6 for Au (5d^*^6p^) which compares favourably with the 
meas-ured value of 'V' 2 in gold compounds. 

I.4.E Shielding Factor (l - R^) for Atomic Quadrupole 
Moments 

13 

Sternheimer has shown that the significantly smaller 

5 

values of atomic quadrupole moments measured on np (n+l)s 
configuration of rare gases by atomic beams involving 
magnetic and electric field gradients can be quantitatively 
explained as arising due to the shielding effects ) of 
(n+l)s-»-d polarization caused by the hole in np shell. 

1.5 Scope of the Present Work 

In overall, it appears that there is overwhelming experi- 
mental evidence in favour of the Sternheimer polarization 
effects especially those corresponding to the quadrupole 
interaction. The experiments on hexadecapole interactions 
are presently being carried out in various laboratories and 
it is expected that similar polarization effects will be 
exhilDited in these cases as well. An accurate knowledge of 
these shielding-anti shielding factors .is therefore ex trariely 
important in the interpretation of the experimental data. 
T^iccordingly, we have investigated in this thesis the sternheimer 
shielding-anti shielding factors of a variety of atoms and ions 
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using consistently one type of wave function, namely, the 

41 

Hartree-Fock- Slater wave function. The consistent use of 
this type of wave function enables us to compare all the 
res\ilts on various systems with some significance. 

The present thesis is divided into seven chapters: 


In Chapter II, we give the details of the method of calcu- 
lations adopted by us. The procedure of generating Hartree- 
Fock-Slater (HFs) v;ave functions for free ions, and ions in 
crystal using Watson model are given. The Sternheimer pertur- 
bation-nxmerical approach is disciassed and the eiq>ressions for 
various shielding-anti shielding factors as given by Sternheimer 
presented. 


In Chapter lil we present the results of our calculation 
and 


of «q/Y« 


for several free ions and ions in crystals. Estimates 
of relativistic and crystalline distortion effects on Y (and 


a ) have been made. We have also studied the effect of crystal- 
q 

line distortion on and R in the case of Ce^'*’. 


In Chapter IV we report our results of Yco and R for several 

ions of interest in Mossbauer spectroscopy. The effect of 

6 ■ ' 

crystal distortions on Y^ and R for 3d ions have been studied 

by correlating these values with the position of the outermost 

maximum in the total radial electron density distribution 

2 

function and the experimental data on e qQ for several iron 



45 


compoimds have been reanalysed using most reliable value of 
Yco and R and a nex-^ range of values for Q v ^Fe) has been 
derived. 

In Chapter V we present our calculations on the hexadeca— 
pole antishielding factor for several ions. Estimates of 
crystalline distortions on have been made in the cases of 
a few actinide ions. We also report our preliminary results 
on the second-order quadrupolar effect in the cases of a few 
actinide ions. A largo value of (1 - ilco) aided with the other 
favourable parameters such as nuclear quadrvpole moment, natural 
abundance, and NMR sensitivity has been considered as a cri- 
terion to identify those nuclei which would be possible 
candidates for studying nuclear hexadecapole interaction. 

In Chapter VI we present the results of our calcuilations 
on atomic shielding factor Rq for np(n+l)s excited state 
atoms of G, Si, Ge, Sn, and Pb. 

Chapter VII summarizes the various conclusions obtained 
in the present work. The scope of future work in this area 
i s al so o utl ined . 



CHAPTER II 


StGrnheimer Perturbation-Numerical Method and the 
Calculation of Shi elding- Anti shielding Factors 
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CHAPTER II 

Stemheimer Perturbation-Numerical Method and the 
Calculation of Shi elding- Anti shielding Factors 

II. 1 Introduction 

The polarizabilities and shielding-antishielding factors 
calculated from Hartree v/ave functions have been found^ to be 
significantly larger in magnitude than those calculated from 
Hartree-Fock (HF) wave functions. The rnarlced sensitivity of 
these one-electron properties to the treatment of exchange 
potential in the zeroth-order wave functions attracted us to 
test the merits and demerits of Slater exchange approxima- 
tion, according to which the nonlocal HF exchange potential 
is replaced by . 

V„ = - 6a { (3/8 tt) (ll.l) 

-<rv 

156 

One definite advantage in using Slater approximation 
is that it reduces the computational time relative to the HF 
case by an order of magnitude without significantly affecting 
the accuracy in the wave functions. In passing we may, mention 
that the Slater approximation provides the only feasible ab 
initio, method for studying one-electron properties of large' 
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molecules and metals. Saxena and Narasimhan^^^ have shovm 

that the diamagnetic susceptibility and nuclear magnetic 
shielding factors for the closed-shell atoms and ions calcu- 
lated from Hartree-Pock-Slater^^^ (HPS) \vave functions corres- 
ponding to a =1 in eq, (II. l) are in esccellent agreement with 

”L 6 O 

HF calculations.' One might therefore expect that the HFS 
wave functions wnuld be suitable for the study of polarizabi- 
lities. 

As regards the choice of the method used for calculating 

the first-order perturbations, we have found from the survey 

of reported results that using the zeroth-order wave functions 

of similar accuracy, the UCHP method due to Sternheimer gives 

the values of the quadrupolar polarization effects which 

compare most favourably with the CHF methods. For example, 

9 4r d 

Lahiri and Mulcherji ' have noticed that the CHF results of 
a and Y<a do not differ as much as a, and Bo, vrhen a comparison 

g 

is made with the other uncoupled calcinations. Sternheimer'^ 
has found that the CHF values of ct^ for up to Argon isoelec- 
tronic sequence differ from the results obtained from the 
perturbation-numerical calculations by 'u 5-15%. In recent 
LCMBPT and DE calculations ‘ on it has been concluded 
that the coupling effects amount to '''15% in most of the cases. 
The Sternheimer metlxfd therefore is an attractive substitute 
for the more difficult coupled calculations for heavier systems. 
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as far as quadrupolar polarization is concerned. Further, 
several reported calculations of a / , R and have been 

generally based on UCHP methods using Hartree or Hartree-Fock 
wave functions. A comparison of I-IF and HFS xvave f -unctions 
corresponding to -the zeroth-order is more justified, if a 
method of similar accuracy is used to obtain the first-order 
perturbed wave functions in both the cases. We have therefore 
chosen the HFS wave functions in our present study and have 
adopted the Sternheimer perturbation-nun eric al method to 
generate the appropriate quadrupolar and hexadecapolar pola- 
rizations in the case of a large number of free atoms and 
ions with the motivation of studying the Sternheimer shielding 
anti shielding factors in detail. We would like to make three 
comments which are related to the work embodied in this thesis 
(l) ' Since the HFS v/ave functions for all the neutral a-txsms 
are already available and since these can be generated 
with much less compuhational effort for most of the 
free ions of interes’t, it is worthwhile to compare the 
HFS polarizabili'ties and shielding results with the 
corresponding HP results and thus evaluate the merits 
of Slater exchange approximation in this context. 

The inter- shell and intra-shell consistency effects can 
be included by using the same moment, and external- 
charge perturbed wave f-unctions, if desired. 


( 2 ) 
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(3) Due to a consistent choice of the zeroth-order wave 

functions, which are already Icnovm to be good approxima- 
tions to the HF v;ave functions, coupled with the use of 
an accurate UCHF perturbation-numerical method, the 
general trends obtained from the study of a large number 
of atoms and ions are expected to be valid even in the 
case of more refined calculations. 

In the previous chapters we have already referred to 
the pioneering work of Sternheimer who made calculations of 
“q' ' cr 2 ' ^0 ^ variety of systems. 

Subsequently, calculations by others have also been reported. 

All these available values have been compared v/ith the present 
calculations, in the following chapters. A fairly recent 
compilation of and R factors for the ground state atoms 
and ions has been given by Lucken. 

The experimental estimates of (l - ) in highly ionic 

solids for the positive ions come out to be larger than the 

theoretical values V7hereas in the case of negative ions the 

estimated values are substantially lower than the theoretical 

counterparts. It is very likely that (l - ) values themselves: 

get substantially modified due to the changes in the outer 

orbitals of the ion in solids. To o'er Icnowledge, there have 

29c 163 

been only two calculations ' reported in the literature, 
both of which use the variation-perturbation method of Das and 
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Bersohn ' and attempt to estimate the effects of crystalline 
potential on polarizabilities and shielding-antishielding 
factors. We therefore thought it worthvzhile to calculate a » 
cr^ and R for several ions using the Sternheimer perturba- 
tion-numerical method and HPS wave functions with Watson 
model to simulate the ion within a crystal. 


II. 2 Calculation of Zeroth-order Wave- functions 

II. 2. A Free Ion 

The Schrodinger equation for a one-electron orbital 
in an atom can be written as (energy in Rydbergs) 


(-v£ + + vp 


Tpi O 

il. . Y , 

1 


(II. 2) 


v/here is the Coulomb potential energy, i.e.. , the potential 
energy of an electron in the field of the nucleus and of all 
electrons, including itself. is the exchange potential 

which can be thought of as a correction for the self- interaction 
energy of the electron. In the Hartree-Foclc approximation. 


V. 


XHF 


E, 


^C + 


(V2 


)§)/x; 


( 11 . 3 ) 


In the theoiry of f ree-electron gas obeying Fermi statistics, 
Slater^^®^ has shown that the exchange potential can be 


= -8 F (n) t (3/8»)} 


(11.4) 


written as 
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where p is the total electron density, and F(n) is given by 


F(n) 



In 


1 + 'L 
1 - 'n ' 


(11.5) 


v/here ri is the ratio of the momentum of the electron to its 
momentum at the Penni energy. In Slater's exchange approxima- 
tion to the atomic systems, F(n) is tahen as 3/4, which is 
the average value of F(n) over all electrons in a Fermi gas 
at the absolute zero of tanperature. The exchange potential 
then becomes, 

V^S = -6 {(3/’3it)p}^/^ (11.6) 

The x° ' s obtained from eq. (ll. 2) using eq. (II. 6) for are 
called Hartree-Fock-Slater (HFS) wave functions. 


We have used the computer prograra of Herman- Skillman 

-159 

(HS) " to generate the free-ion wave functions required in 
the present thesis over 441-point mesh. 


II. 2. B Crystal Ion 


The crystal ion wave functions have been generated by 

including an extra stabilizing potential V shown in Fig. II. 1 

164 

which is due to Watson. In Watson model, the ion with an 

electric charge n^^^, is sxirrounded by a hollow sphere of radius 

r. and carries a uniformly distributed charge of -n. o, units, 
ion ^ ^ ion ^ 

! 1 f K /' NPUR 

The Coulomb potential V_ in eq. (II.2) then becomes; ' 


AiX, Nt>. 


A 


5QS55 
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(II- 7) 


for 

r ^ 

r . and 

ion 



II 

V_ + 2a 

C r 

(II. 8) 


for r ^ ^ion' ^^spectively. 

T 65 

ThrouglTout this wrlc, r. has been taken from Pauling's 

ion 

jL 2 -f- 3 + 

table of ionic radii, except in, the cases of Fr , Ra , Ac , 

At and rare earths and actinides where Zachariasens values^^^ 
have been used. 

II. 3 External Charge-Perturbed Inhomogeneous Differential 
Equation; Quadrupolar Perturbation 

In the following paragraphs we shall present the method 
of computation adopted by us v/hich is essentially the same as 
that prescribed by Sternheimer”^^' in his original papers. 

The first-order perturbed Schrodinger equation can be 
written as 

( H - E ) x ^ = (E. - H ) x° (II. 9) 

O o 1 1 1 1 

where x? represents the first order perturbation of the unper- 
turbed. wave function X?; H and E correspond to the unperturbed 
Haiailtonian and eigenfunction for the unperturbed state n . The 
unperturbed wave function can be written as 
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u (nl) 

x° = ^ 

1 r 1 


(ll.lOO 


where v.^(nl) is r times the radial part of the lonperturbed 
function. '>J-^(nl).is norraalized according to 


wave 


/ ^i^(nl) dr = 1 

^ o 


{II. 11') 


The crystal potential at the ionic site can be 


written' 


V , , = I ^ r^- 
crystal ^ n n 

^ n,m 


(II. 12a) 


where A are the constants characteristic of the ionic distribu- 
n 

tion in the lattice and is the noiraal ized spherical harmonic. 
The . term corresponding to n=2, and m=0, is the quadrupolar poten- 
tial given by 


H = _A° r^ Y° 

1 '^2 ^ 2 


(II. 12b) 


Under the influence of this perturbing field the spherical 
shells get polarized according to, 

)P(1=0)-^ X?(1=0) + x^ (1* = 2) } 

)^(l=l) {x?(l=l) + ’>^■(1' = l) + (1' = 3) } 

^(1=1) {x°(l) + xV (1’ = 1) + xj; (1' = 1 + 2) , 

+ xj (1' = 1 - 2) } (11.13) 

The part of the polarization characterized by the same 1 value 
as the unperturbed v^ave function is called the radial perturbation 
and those with Al= +2 are called angular perturbations. 
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The perturbed wave function can be separated into radial 
angular variable as 


1 ^1 1 ' ) _ 

x: (1 = 1') = 

1 ■ r 1 


(IT. 14) 


Sifastituting eq- (II. 10) and (II. 14) in eq. (II. 9) and separa- 
ting the variables, u' can be obtained from the following 
inhomogeneous radial Schrodinger equation 


{ - + V - E } u' (nl-v 1 ' ) 

2 O O 1 

dr r 




(11.15) 


where -f- Vg and (nl 1 ‘ ) represents r times the 

perturbed radial v/ave function corresponding to the polarization 
(nl -♦-1'). According to Sternheimer local approximation. 


(V + V - E ) - -i— “ — -- ~ 

^ C + ‘XS o u ,2 2 

o dr r 


(11.16) 


We have solved the finite difference equivalent of eq. (II.15) 
by writing eq. (II. 16) as 

u’(r+6) - 2.u'(r) +u'(r-5) i n j-iT 
V (r) - E = -o . UA+IJ (11.17) 

o 

and substituting in eq. (II.15), 


(r+5) = u^ (r) {2. + 5 


2 l' (l'+l) - 1(1+1) 


u' (r+ 6) - 2u’ (r) + u ' (r-6 ) 

o _ o , . _t> , 

u'(r) 

. o, ■ ■ 

- (r-5 ) - u^ (r) ^ {r^-<r^>^^<S (II. 18) 
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We have used this form for the outward integration which 

requires the knowledge of the unperturbed wave functions and 

the perturbed v/ave functions at tvro initial points to start 

the direct solution. Tire two boundary conditions (O) and 

(»') = O define the solution uniquely. The starting slope of 

the solution is changed in an iterative way till the solution 

at large r decays exponentially. The behaviour of the solution 

for the first few points on the r mesh when fed in the computer 

results in a faster convergence. Three to four iterations are 

generally required to achieve convergence. For A1 = +2 the 

behaviour of (r ->-o) has been taken as r^"^^ whereas for 

1-1 

A 1 = -2 and it is taken as r . 


For large r the inward integration has been started with 
a multiplicative constant in u^ (r) according to the relation 


(r+A) = expE -{W (r)} 5] 


(11.19) 


where 

N(r) 


u’ (r+5 ) - 2u' (r) + u' (r-5-) 
o o o 



1' ( i; +1) -lii+i) 

I ‘ o 


ru' (r) 

p 

(r") 


(II. 20) 


In eg. (11.19) to start with we have taken N = -e . The multi- 

, o 

plicative constant is varied till the boundary condition at 
r=0 is bracketed within a desired limit. 



57 


We have employed both the inv/ard and outward integrations 
and matched them iteratively near the outer peak of the unper- 
turbed wave functions for all the excitations. 


II. 4 Moment -Perturbed Inhomogeneous Differential Equation: 
Quadrupolar Perturbation 

In an alternative approach the source of perturbation is 
the nuclear moment itself. 


H 


^ _0 (3 cos^e - 1) 


( 11 . 21 ) 


2r' 


The first-order perturbed differential equation takes the formj. 


. ^1' (I'-M) 

i • “ '2 + ■ - ■ ~— 

dr r 


+ V - S } = h { - <^> 


o ' -3 

r 


3 nl 


y(ll. 22) 


which has the equivalent difference equation as 


u. 


(r+5 ) = (r) { 2. + 


,2 l'(l'+l) - 1(1+1) 


+ 


u (r+6 ) - 2u (r) + u (r-6) 

__p__ p j 


- u, (r-6 ) - u (r) ( 6. { - - < 5^^) (11.23) 

J- o r j-3 hi 11 

The metlTod of solving the moment-perturbed equation is exactly 
the same as the external charge-perturbed case, except that the 
boundary condition near the nucleus is different. Thus, in the 
case of s -»■ d type excitation 
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u (O) = ^ 

1 o 


(11.24) 


where 


u (ns) 


{ 


} 

r-+-0 


(11.25) 


In the case ofrp-^f, (O) is zero and near the origin it 
behaves as 


where 


b r 

(r small) = { - } 


u (np ) 

b = { -P-™- } 

r ^ r.>.0 


(II. 26) 


(II. 27) 


For other excitations ^.^(O) = 0 and its behaviour near the 
nucleus has been similarly obtained from the appropriate loga- 
rithmic series. It should be noted that the constants a , b , 

o o 

are calculated in the Herman- Ski Urn an routine and then 
recalculated in the independent perturbation program to checlc 
that tbe wave functions have been correctly transferred. 


The invrard integration in the case of u^(nd-»-s) perturba- 
tion was very unstable and this is one of the reasons why we 
have chosen to match inward and outward integrations for the 
angular cases as well. 

II. 5 Moment-Perturbed Inhomogeneous Differential Equation: 
Hexadecapolar Perturbation 

In the case of hexadecapole moment-perturbed wave function - 

u „ the inhomogeneous eauation can be written as, 

1 ", 
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{ 


dr 


2 


1 ' (1 ’ + 1 ) 


^ +V - E>u = U{v -< i-> 5^^J(II.280 

2 o ol/H o5 

r r nl 


and only nd d and nf f perturbations have been considered. 

For nd d excitation, the boundary condition at r=0 is given by 

3 

c/6, vrhere c is the coefficient of r in the power series 


expansion of nd vrave function near the nucleus. 


II. 6 Method for Integral Evaluation 

The 'presentation mesh' of HS has been used as the 

'integration mesh' to perform the block-vrise integration through 

167 

adjacent inteirvals according to 



x=a 


F (X) 



6[_|{F(X.^l) 


+ F(X J}- 


12 


i (42 

j=i 



+ F(Xj)} 


+ 


11 
7 20 




F(X , ) + A^(X )} ] 

J '‘'-L J * 


(II. 29) 


In order to calculate the differences up to fourth it is 
necessary, that the two extra entries at both the ends of each 
block are included. This sequence of entries gives J=41 points 
in the first blocl: and J=44 in the rest. This is achieved 
before the integration is started and the final result is of 
course the sua over all the blocks. 
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Table II. 1; The various coefficients pertaining to the quadru- 
polar polarization, obtained as a result of 




integration over angular 
over magnetic substates 

variables and 

suffnmation 

G(nl->1, 

; 1 ■) 
e 

c(nl^l') 

C (L^) 

C(L2) 

C(L.3) . 

(ns-»d; 

p) 

3/5 

4/3(1) 



(np-»p; 

p) 

43/25 

4(0) 

4/25 (2) 


(np^f ; 

p) 

72/25 

36/25(2) 



(nd.vd; 

p) 

16/7 

4/3 (l) 

12/49 (3) 


(nd->-g; 

p) 

144/35 

72/49(3) 



(nf -»-h; 

p) 

16/3 

40/27 (4) 



(nf-»■f ; 

p) 

224/75 

24/25(2) 

8/27(4) 


(ns->d; 

d) 


4/5(2) 



(np->p; 

d) 


28/25 (l) 

36(175(3) 


(np-vf ; 

d) 


12/25 (1) 

144/175 (3) 


(nd^d; 

d) 


4(0) 

-12/49(2) 

16/49(4) 

( nd-)-g; 

d) 


144/245 (2) 

40/49(4) 


(nf-»-h; 

d) 


40/63 (3) 

80/99(5) 


(nf-»-f ; 

d) 


48/25 (1) 

-88/225(3) 

40/99(5) 

(ns-»-d; 

f) 


4/7(3) 



(np^p; 

f) 


103/175(2) 

4/21(4) 


(np->f ; 

f) 


72/175(2) 

4/7(4) 


(nd-»-g; 

f ) 


12/49(1) 

24/49(3) 

300/539(5) 

(nd^; 

f ) 


72/49(1) 

-44/147 (3) 

500/1617(5) 

(nf->-f; 

f > 


4(0) 

7 6/225 (2) 

-4/11(4'), 
500/1287 (6) 

(nf-»-h; 

f ) 


20/63(2) 

40/77 (4) 

7 00/1287 ( 6) 
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electric field at the valence electron site will be referred 


to as of„ and is given by. 


0 . ‘(nl-vl') = c(nl-»-l') ■ 


D 


(11.34) 


where f (r) = / 


^n 1 
e e 


^ J u u' r‘^dr’ + r^ /“u’u'r' ^dr' (11.35) 
^4 '' o 1 ■' O 1 

r o r 


and w(n^l ) denotes the unperturbed valence electron wave func- 

vj G 

tion. The corresponding exchange contribution v/ill be referred 
E* 

to as 02 is given by 

E 


a 2 '^(nl^l’;n^l^;L) = ~C(nl->l’;L) "^^2. 

^ ^n 1 

e e 

CO 


where 

K 

= / u ' 

with 

G 

1 



j.L +1 



L 

+ r 


e e 


,L 


(II. 36) 


e e 


(11.37) 


The coefficients C(nl l';n_l ;f') are given by, 

CO 


G(nl -d. ' ; ngl^;L ) 

4 T dm. I’m) 1 m^) C ' (1 ’m; 1 md 

m=-l - ® ® ® . 

d m ;l m^) 

® ® ^ ® (11.38) 

The integers L correspond to the expression in spherical 
harmonics of rj ^2 in the exchange interaction. Tlie coefficients 
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C have been compiled in Table II. 1. 

The total is given by 

o^(nl-^l') + l o ^ (nl-»-l’; n^l^,L)} (11.39) 

^ nl L ^ ® 


II. 7, B Calculation of and R from the Moment-Perturbed 
Equation 

From the quadrupole moment perturbed wave function 
■u.^(nl-^l) the Y„ can be expressed as. 


Y = T c(nl-+l‘) / u- (nl-^l') u (nl) r^dr 

00 « 


(11-40) 


nl 


Denoting the valence electron wave functions by w(nglg), the 
direct contribution to R, by (nl-»-l’) excitation of the core 
orbital is given by „ 


RD(nl-^l'; n^l^) = 


/ Y' (nl-4-1 ’ ; r)^ w (n 1_)^^ r~^dr 

0 0 


-3 

< r > 


(II. 41) 


n 1 
e e 


vj'here 

Y(nl->l’;r) 


with (nl"*"! ' ) 


= X t / Q.(nl-^l') dr’ + r^ / Q. (nl-»-l*)r'”^dr '} 

o ^ r ■ 

(11.42) 

= c(nl-*'l ' ) - Q u(nl) (nl-*-! ' ) r^ 


which gives the direct contribution to the density of the 
induced moment corresponding to (nl-*-!') excitation. The 
contribution to R, by the exchange term arising from (nl-»-l') 
is given by 
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C (nl-> 1 ’ ; n 1 ; Ll « 

Rp(nl->-l';n 1 ;L) = - ------ ® /u (nl)v7(n 1 )g dr 

<r“^> o ^ e e 

® ® , (II. 44) 


1 ■ ^ 

V 7 here gj^ (nl -^1 ’ ; n^l^) = -^-- 3 - / u (nl->-l') w (n^l ) r'^dr' 

CO 

+ /u (n> 1 ' ) w (n^l^) r^“^ dr‘ (11.45) 

„ X 0 0 


the net R is given by 


R = I {RQ(nl-»-l'; n 1 } +l Rp.(nl^l‘; n 1^,L) } (11.46) 

nl ® ® L ^ ® ® 

The constants C (nl-^-l ' ; n^l^;L) are given in Table II. 1. 

0 0 

II. 7. C Calculation of ’i„ from the Moment- Perturbed Equation 

The hexadecapole antishielding factor has been calcu- 
lated as. 



11 

8 

(nd->d) + I 

(nf-»-f) 

(11.47) 


nd 

nf 



where. 

(nd-»d) = 

80 r ' f i\ 

^ JUo(nd) 

u ' (nd-^d ) r'^ dr 

i , h 




0 

CO 


(11.48) 

and 

(nf-»- f ) = 

~ /u^(nf) 

O’- 

u ' (nf ->£ ) r^^dr 

X / A I 



The net angular contribution to hag been estimated by using 
the apprcximation 


, ang _ 9 _ 

~ 5 

, ang 


(11.49) 
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where are calculated from moment-perturbed equation 

( 11 . 22 ). 

II. 3 Typical Results of Calculations 

For illustration purposes we present here a few typical 
numerical and graphical results. Thus, our numerical values of 

. 3 

u^(5p-».p) for Trn ' have been given in Table II. 2. This function 

has been represented graphically along with 4 x u^(5p) in 

Fig. II. 2. The free ion HFS wave functions for -were 

generated in •v iJj min and single iteration in the case of u^CSp^-p) 

took '^2 sec in our IBM 7 044/1401 computer syston at the Indian 

Institute of Technology, Kanpur. We have given in Fig, II. 3 

the plots of u.(5d-»-d) along wdth 2500xu (5d) in the case of 
1 , B o 

Au"*". 

II. 9 Computer Programs 

The computer listings of the following three main programs 
used in the present work have been given in .Appendix I; (a) - the 
free and crystal ion HFS wave function program, (b) y„ and R 
routine, and (c) polynomial fitting program. 

The quantitative estimates of the effects due to crystal 

distortions on shielding-antishielding factors have been 

obtained by expressing these in terms of suitable polynomial 

—3 

series in r . <r > and p , the position of the outermost peak 
ion m 
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Table II. 2; The moment perturbed radial wave f unctions u^(5p-»-p) 
as a function of r (in a.u.} for Tm^'*’ ion. 


r (a. u. } 

u^ (5p-»-p} 

r (a.u.) 

u^ (5pH-i 

0. 0000 

0. 0000 

1.2563 

11.1452 

0. 0048 

1.5539 

1.5671 

20. 4203 

0. 0097 

2.6399 

1.8779 

21.7037 

0. 0016 

3.1168 

2.1887 

18.9976 

0.0019 

3.0298 

2.4-996 

15.0938 

0.0183 

3.0932 

2.5341 

14, 6514 

0.0280 

1.7986 

3.1553 

7.7172 

0.037 3 

-0.5717 

3.7774 

3.5415 

0.047 5 

-3.2644 

4.3991 

1 . 4893 

0.057 2 

-5.7451 

5.0208 

0.5955 

0.0583 

-5.9917 

5 . 0898 

0.5375 

0.0777 

-9.0211 

6.3331 

0.0724 

0. 097 1 

-9.1619 

7.5765 

0. 0087 

0.1166 

-6.8807 

8,8198 

0, 0009 

0.1360 

-3.1008 

10,0631 

0.0001 

0.1331 

-2.6295 



0.1770 

5.8894 



0.2153 

11,9292 



0,2547 

13.9431 



0.2936 

12. 2773 



0.2979 

11.9121 



0.3756 

2.1018 



0.4533 

-9.0335 



0.5310 

-17.0507 



0.6087 

-20.8103 



0.6173 

-20.9793 



0.77 27 

-18.1110 



0.9282 

-8.6681 



1.0836 

1.7156 



1.2390 

10.3431 













U n p© r t u r be cj function 
(BOdx Uo(5d)) 
Perturbed function 


'JY'! /- 








Wmmi 
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in the total radial charge density distribution of the atom 
ion. The various polynomial fittings have been carried out 
by Using standard method of minimizing the average percentage 

1 gp 

error via normal equations. ^ The listing of the computer 
program for this purpose is given in Appendix I under polynomial 
fitting prograra. The routines corresponding to (i) the inward 
and outv/ard integration, (ii) the integrals, and (iii) the 
polynomial fitting have been developed by the author. 



CHAPTER III 


Steriilieimer Quadripole Shielding-Antishi elding Factors 
for Several Closed Shell Ions, A Rare Earth Ion and 
Ions in Actinide Series: Estimates of Relativistic 
Effects and Crystalline Distortions 
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CHAPTER III 


Sternheimer Quadrupole Shi el ding -Anti shielding Factors 
for Several Closed Shell Ions, A Rare Earth Ion and 
Ions in Actinide Series: Estimates of Relativistic 
Effects and Crystalline Distortions* 


III.l Introduction 

From a survey of the reported theoretical calculations 

till 1971, it became clear to us that extensive calculations 

of and utilizing nonrelativistic HFS wave functions for 

free ions and ions in crystals represented by the Watson 
1. 6 4 

model, adopting the Sternheimer perturbation-numerical method 
will serve the following two pxorposes: 

(l) a comparison of our and values v/ith those of Feiock 
and Johnson^^^ (Fj), v/ho have used relativistic HFS wave 
functions and the uncoupled method similar to Sternheimer, 
v7ould give an estimate of net relativistic-effects on 
and a , and 

The work presented in this chapter has been published in 
the following papers: 

(l) K. D. Sen and P. T. Narasimhan, Adv. Nucl. Quadrupole Res. , 
1 , 277 (1974); (2) Phys. Rev. B (in press) ,(1976) ; (3) Phys.' 

Rev. A (in press) (1976). 
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(2) a comparison of our free ion and ion in crystal values 

woTold give an approxim3.te estimate of the changes brought 

about in and a due to the lattice potential. 

'Z 

The interpretation of nuclear quadri:pole coupling data 
is generally based on eq. (1.5). It is therefore necessary 
to make estimates of the net relativistic and lattice poten- 
tial effects on the shielding-anti shielding parameters. 

Ill, 2 Results and Discussion 

III. 2. A Y„ and a for Closed-Shell Free Ions and Ions in 
Crystals^ 

We have considered the following six sequences of closed- 
shell free ions isoelectronic with He, Ne, A, ICr, Xe, and Rn 
respectively. 


(i) 

Li"^, 

iD0 / ID ^ 

C^+ 


0®+; 


(ii) 

^ / 


. - 2 + 
Mg , 

Aa2+, 

Si^^ 



(iii) 

Cl", 


Ca^ h 

Sc^'^, 



Cr^+ 

(iv) 

Br", 

Rb"^, 

Sr2+, 

1 Yt^+, 

Zr^"^ 

■, Nb^ 

Mo 

(v) 


Cs^, 

Ba^’^, 

La^"^; 

and 



(vi) 

At“, 

Fr'^, 


, Ac^"^. 





With the use of nonrelativistic HFS wave functions for 
these ions we have generated the external charge-perturbed 
wave functions and calculated 7^ and ct according to the 
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method, described in Chapter II. The calculated values of ^ 

CO 

and a have been oresented in Table III.l and Table III, 2 
q 

respectively. For purposes of comparison, the values calculated 
by FJ have been given in column 3 of these tables. As expected, 

our (or a ) values in the case of lighter ions (Z<18) are 

q “ 

in excellent agreement v/ith FJ calculations. For the heavier 
ions FJ values are consistently larger than the presently 
calculated values. Comparing the negative halogen ions in the 
two cases v/e have obtained the upper limit of the net relati- 
vistic effect on Y in each series. It has been found that 

CO 

the relativistic effects increase the nonrelativistic Y<a value 
by 7% for the ions up to Kr isoelectronic series. For Xe 
series this effect becomes as high as'''l9%. 

The results of other available calculations of loa and a 

q 

have been compiled in column 4 of Tables III.l and III. 2, 
respectively.. Our values for the positive ions are in good 
agreement with those of Sternheimer and Lahiri and Miilcherji 
(see the references appearing at the end of Table III. 2), 
which are based on HF V7ave functions. For P”, Cl“ and Br” 
the present Yc values are -41, -82 and -195 respectively which 
are significantly larger than the antishielding effect obtained 
from HF v7ave functions. It would be useful to modify the Slater 
exchange potential in. the outer regions to make the wave func- 


tions more internal. 



Table III. 1*: Total Y„ values for the free ions and ions in crystals. The crystal ion 

values correspond to Pauling ionic radius. 
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In the case of multi negative ions corresponding to He, 

A, Kr, Xe, and Rn series it has not been possible for us to 

generate the self-consistent v/ave functions. Such a situation 
17 0 2— 

is V7ell knovm- in the case of 0 where the stable free ion 

does not exist and it is doubtful that a HP solution would 

converge to a state with all the 10 electrons bouud to it. On 

the other hand oxides, sulphides etc. exist as stable ionic 

solids and the approximation of ions with a certain radius 

is a useful physical concept which explains many properties of 

ionic solids. In order to obtain the wave functions for ions 

in crystals it is necessary to include the stabilizing po ten- 

171 

tial due to the lattice. Yamashita and coworkers have 

2 — 

obtained the wave functions of O in MgO crystal by a variational 

metliod. Watson has proposed a spherical potential surrounding 
2 - 

O ion according to the method described in Chapter II. As 

compared to the free atomic or "ionic cases these v/ave functions 

give better agreement with the experimental properties such as 

1 G 2 

dipole polarizabilities, diamagnetic susceptibilities, atomic 

171 

scattering factors, and Compton profiles. Pulcamachi and 

17 2ci 

coworkers have shown that the contracted wave functions of 

2 - ' 

Yamashita and cpX'^rorkers and Watson for O give more or less y 
similar results for atomic scattering factors and Compton 
profiles in Migi). Burns and Wikner^®*^ and Paschal is and Weiss^^ 
have earlier calculated 1^,, values in a few cases using the 
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164 27 

Watson model and the variational method of Das and Bersohn. 

We thought it worthv/hile to study these ions as well as other 

ions of interest to the experimentalist by means of a more 

accurate perturbation-numerical procedure and using HFS v/ave 

functions. 

With the use of the Watson model^^"^ to describe the 
lattice potential, V7e generated the SCF wave functions for all 
the ions considered atove and additional 10 negative ions given 
by 0^", P^“, Si‘'"; Se^", As^", Ge^"; Te^“, Sb^“, Sn'^". 

To our knowledge, the self-consistent numerical HFS wave f'onc- 
tions for these ten ions have been generated for the first 
time. We have used these wave functions to calculate Y«, and a 

q 

values by means of external charge -perturbed procedure and the 
results obtained have been given in column 5 of Table III.l and 
Table III. 2 respectively. For the purpose of comparison we 
have also included the y„ results obtained vrith the Watson 
model by Burns and Wikner^^'^ (BW) and Paschal is and Weiss^*^^ 

(PW) in colum.n 6 of Table III.l. 

For heavier positive ions the results based on the 

27 

variation-perturbation metliod of Das and Bersohn underestimates 
the antishielding effect. This has also been observed in the 

ease of free ion caloulations. In the cases of k 7 Ca^'*', Sc^’*’, 

4+ 163 . ■ 

and Ti , Paschalis and Weiss have quoted as respectively 

-22.83,-20.58, -20.34 and -25.51. This trend seems to be 
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4 -+ 

erroneous since one expects the 3p v;ave function in Ti to be 

3_l 

more internal than that in Sc . As compared to our free ion 
values the in crystal is''' 50% smaller for the negative ions. 

For the positive ions also a similar increase is obtained in 
a large number of cases. Thus the effect due to crystalline 
distortion is estimated to be *''3 times larger than the relati- 
vistic or coupling effects discussed earlier. 

Ejcperimental estimates of are available for a large 

number of ions in ionic solids (vide Chap. I, Sec. 4. A). The 

effective Y„ for ITa*^ has been estimated to be '''-9. The present 

92 

value of -7.7 compares favourably with this value. Raymond 

34 * 

has estimated the best fit value of 7^ for A1 in Al2Si0^ 

polymorphs as -4.9 v/hich is in satisfactory agreement V7ith our 

135 

crystal ion value of -5.7. Bhide and Hegde have combined 

the electric field gradient data on '^^sc and ^"^Pe to estimate 
3 + 

Vjo --19 for Sc which is to be compared with the present 

125 

value of -20. Recently .-jiderson and Karra have estimated 

Y^ = +.4.4 for Li'*’ in LiF single crystal via acoustic-saturation 

studies. The interpolated result of 0.34 due to Lahiri and 
2 ^ 

HrOcherji * and the present value of 0.28 are significantly 1 

155 

low’-er than the above experimental estimate. Gupta and Sen 

2 1 

have calculated Y„ for Li in Is 2p configuration and noticed 
that the theoretical value of -4.4 so obtained is in reasonably 
good agreement with Anderson and Karra's estSnatein LiP cirystals. 
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However, such an electronic configuration for the lithiim ion 

is highly iraprohable in solid LiF and the two values cannot be 

meaningfully compared. Vie have found that in an analogous 

experiment on LiP single crystal, Shutilov and cov7orkers^^^ have 

estimated to be 0.5+0, 3 v/hich is very close to the crystal 

ion value predicted by Lahiri and Mulcherj i, Paschalis and 
, , 163 

■'leiss and by us. It therefore appears that overall agree- 
ment between theoretical and ej^perimental estimates in the 
cases of positive ions in ionic solids is quite satisfactory. 

29c 

Burns and Wikner have estimated the values of for 
Cl , Br” and I~ to be centred around -10,. -35 and -45 respec- 
tively. As mentioned in Chapter I, this choice has been 

90 ■ 

criticized by Das and Dick on the grounds that the models 

used to interpret the experimental results in these cases are 

121 

rather crude. . However, Ikenbery and Das have subsequently 
found that a fairly contracted in the cases of halogen 
negative ions v/ould be required to interpret the quadrupole 
coupling data in al’cali halide solid-solutions. Since our free 
ion values for halogen negative ions significantly over- 
emphasize the antishi elding effect as compared to HF results 
we cannot directly compare our results with the experimental 
estimates. Beri et al. have recently included the consis- 
tency corrections to Yo, values based on HF wave functions for 
F~, Cl” and Br” ions. Assijming the net contraction effect to 
be similar in the cases of HF and HFS approximations, we have 
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estimated the contracted HP value of 7^ for F", Cl“ and Br” 
to be - 5 , -23 and -74, respectively. In view of the recent 
estimates of Yjo a, -30 for Cl” in CoCl^. 2K20^^*^ and alkali metal 
hexachlorostanates, GdCl^,^^^ the agreement can be considered 
as satisfactory. He thinlc that Burns and Wikner's estimate 
of contraction of Yea for halogen negative ions needs a careful 
rechecking. 

Our nonrelativistic values are in excellent agreement 
with the FJ values for the lighter ions. In the cases of heavy 
positive ions the present nonrelativistic values are consistent- 
ly larger than the corresponding relativistic values. This can 
be qualitatively explained as arising due to the relativistic 
contraction of the outer loops of the external p orbitals which 
leads to decrease in relativistic ct values. For the mono- 

q 

negative ions however, our values are slightly lower than the 
relativistic values of FJ. We do not have a satisfactory 
explanation for such a trend observed in the cases of negative 
ions. 

Due to the easy defo rm ability of the negative ions, it 
is natural to eiqDect a range of ionic radii for them. We have 
studied y_ and a values as a function of the radius of the 

■ q 

Watson sphere for F“, Cl”, Br”, l”, O^”, S^”, Se^”, and Te^ . 
These results have been shown in Table III, 3. It is evident 
from these results that y and a values critically depend upon 

CJ 
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the choice of the ionic radius. With increase in this radius, 

as expected, the and a values approach the free ionic 

result. We have obtained polynomial-fittings in each case for 

Y and a , as a function of the ionic radius r. . The poly- 
q ion ^ 

nomial 

2 or 3 

Y„(or a^^') = 

has been used and the various coefficients have been listed 

in Tabic III. 4. From such a correlation it is found that a 

typical reduction by'v 20;4 in the Pauling ionic radius results 

values of 

in 15-20>o and 25-30% decrease respectively in/mo no -negative 

and di-negative ions considered here. In the cases of negative 

ions due to their easy deformability, it is very likely that 

the size of the ion is different in different ionic solids. 

As estimated above, the values undergo significant changes 

as a result of the variation in r. . It would be highly 

desirable to have estimates of in a given systom by so-ai- 

empirical means so that more appropriate values of Y^ can bo 

2 

used in interpreting the experimental e qQ data in solids- 

'V',. ' / 

III.2.B ^ ^ in Crystal 

151-153 

The experimental estimates of the Sternheimer 

shielding-antishielding factors for rare earth ions in solids 

'''''''' X 

are generally higher than the corresponding free ion values. 
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We therefore thought it ^^^■orth■while to estimate the effects of 

crystalline potential on Y , a r and R for the rare earth 

“ q 2 

. ^ „ 3 + 

ion of Ce 

We have used the Watson model to calculate'^cs , a , 

2 

and R for Ce'^'^ to estimate the effects due to the lattice 

potential. The calculations ofY„ / a , and a were carried 

q i 

out using external charge-perturbed v/ave functions. In the 

calculation of R vre have used the moment-perturbed procedure. 

All the direct and exchange terms have been included in the 

calculation oi- o ^ and R. These results have been given in 

Table III. 5. Our free ion values of a , Y„ , <^2 ^ 

3-^ 

Ce ' are 1.26, -65, 0.71 and O. li respectively. These are in 

1 2b 

quantitative agreement v/ith the results of Gupta and Sen 
who have used the HPS wave functions and Sternheimer perturba- 
tion-numerical method to calculate Y„ for the free 

tripositive rare-earth ions. Our crystal ion values of 
Yo, <■ 02 ^ Ce^"^ are 3.52, -107, 0.53 and 0.2 respec- 

tively. In the cases of a , Y^^ , and R values we observe a 
significant increase over the free ion values. The decrease 
of crystal ion case is not clear to us. We note 

151-153 

here that the experimental estimates of <3^ are consis- 

tently 40-50% higher than the free ion values. We are 
studying in greater detail, the effects of lattice potential 


on 02 * 
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Table III. 5; The individual contributions to a ^ y , and R 

3 , q 'o 2 

in the case of Ce in crystal. Pauling ionic 

radius has been taken as 1.12 


Exct. 

“q 


^^2 

R •• • 

ls-+d 

0. 0000 

0.0116 

O. OOOO 

0. 0116 

2s->d 

0.0000 

0.0269 

0. 0000 

0.0248 

3s-k1 

0. 0000 

0. 0487 

0.0008 

0.0286 

4s-k3 

0. 0013 

0.0345 

0.017 8 

0.0150 

5s-k1 

0.5592 

0.2768 

0.2114 

0. 0039 

2p-*f 

0. OOOO 

0.0339 

0. 0000 

0. 0350 

3p-»-f 

0. 0000 

0.0808 

0.0006 

0.0472 

4p-»-f 

0. 0029, 

0.1710 

0.0226 

0.0265 

5p-^ 

2.7577 

0.6608 

0.4197 

0.0050 

3d-»s 

0.0000 

-0.0279 

-0. 0007 

-0.0159 

4d-»-s 

0.0011 

-0.0891 

-0.0199 

-0.0089 

3d-»g 

0.0000 

0.0915 

0.0005 

0. 0606 

4d.»g 

0.0058 

0.2301 

0.027 0 

0,0391 

Total (Ang. ) 

2.3317 

1.534 

0.3529 

0.2720 

2p.^D 

0. 0000 

-0.2462 

O, 0000 

-0.2084 

3p-».p: 

0.0000 

-1.4521 

0.0000 

-0.3646 

4p^p 

0.0006 

-8. 0337 

-0.0013 

0.1011 

5p^p 

0.63 84 

-95.805 

-0.1710 

0.4275 

3d-Hi 

0.0000 

-0.3244 

0.0000 

-0.0332 

4d-K3 

0.0019 

-2.5267 

-0.0046 

0.0016 

Total (Rad. ) 

0.6409 

-108.388 

0.17 69 

-0.07 60 

Grand Total 

3.5226 

-106.854 

0.5298 

0.1960 
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III.2.C for Actinide Free Ions and Ions in Crystals 

The quadrupole coupling rneasuranents on actinides have 

recently become available mainly through Mossbauer effect^ 

and atomic beam magnetic resonance^ methods. In this 

region the first calculation of was carried out by 
173 2 + 

Sternheimer on Am . He used the HFS wave functions corres- 
ponding to the neutral atom and iDredictcd that these functions 
overestimate the antishielding effect by a, 5%. Feioclc and 

1 r \ ■Nr -4- 

Johnson (Fj) have calculated Yco for Th and U using 
relativistic • HFS vrave functions. 

With the aim to estimate relativistic and lattice poten- 
tial effects we have calculated free ion values of for all 
tetrapositive actinide ions and crystal ion for Pa'^'*', 
and Lw" . We have used the moment-perturbed wave functions 
in both cases and the results so obtained have been presented 
along with the other available results in Table III. 6. In 
view of the fact that the actinide elements show variable 
valency we have also calculated Yc» for and 

The 5 f contribution has been included by considering the 

occupancy of 5f shell. Aiter the completion of our work, we 

155 

came across the paper by Gupta and Sen where these authors 
have calculated for Th'^f Pa'^'f Np'^'^ and Pu^"^ 

using external charge-perturbed wave fiinctio ns but without the 
inclusion of 5f contributions. A.1 though we used the moment- 
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Table 1x1.6. 


Tlie non-relativistic total values to zeroth order 
for the ions in actinide series (90 _<103). The 

electronic configuration considered in each case 


has been shovm in column 3. 


Ion 

Z 

Configuration 

Free 

(Presen't) 

-''°5ther-sT Ion 

Th"^'’' 

90 

(Rn)5f° 

-107.17 0 

-177.5^ 

-108.46^ 

Pa^'^ 

91 

(ii.n)5f'^ 

-105.325 

-104.9^ -145.23 

,j3 + 

92 

(Rn)5f^ 

-117.937 



92 

(Pji)5f^ 

-103. SOS 

-103.0^ 


92 

(Rn)5f^ 

- 95.530 


u6+ 

92 

(Pn)5f^ 

- 85.155 

-143.9^ 

- 88.3^ 


93 

(Rn)5f^' 

-102.170 


Pu^+ 

94 

(Rn)5f'^ 

-lOO. 970 

-100.4^ 


95 

(Rn)5f’^ 

-132.867 

-137.3^ 


95 

(Rn)5f^ 

- 99.786 



96 

(Rn)5f^ 

- 98.842 
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(Rn)5f'^ 

- 97.935 

-141.7 

Cf^'" 

93 

(Rn)5f''^ 

- 98. 890 


e4+ 

99 

(Rn)5f^ 

- 96. 47 



100 

, (Rn)5f^° ^ ^ ^ 

- 96.40 


Md^+ 

101 

(Rn)5f^^ 

- 95.417 


M 4+ 
No 

102 

(Rn)5f^^ 

-95.210 


T 4t 
Lw 

103 

(Rn) 5f 

- 94.211 

-139.7 

^ See 

reference 168; ^ see 

reference 173; 

See ref erence 146. 
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perturbed v/ave functions, the agreement in the case of 

is to within fv 3 i'o and this provides a good check on the accuracy 

of both these calculations. It is gratifying to note that the 

use of ionic v/ave functions actually reduces of by 

173 

as has been earlier predicted by Sternheimer. The 

for heavier tetraposi tive actinide ions in free state seem to 

be Constant at -95 and for actinide series of ions can be 

tak.en as -100(+5'). This is not xonexpected since like the 

tripositive rare earth ions, the tetrapositive actinide ions 

o 

are characterized by an ionic radii of 0.9 A. 

169 4+ 6“f- 

A comparison with FJ values for Th and U shows 

that in this region of the periodic table relativistic effects 

would increase Y by ''>65/0. Freeman and Watson^ have 

CO 

indicated that such an in'eroase might result because of a 
severe relativistic contraction of the inner loops of p shells, 
Hovever, one expects the contraction in outer loops of p shells 
duo to relativistic effects and a detailed analysis of the 
relativistic two-component unperturbed and perturbed functions 
is necessari'' before arriving at any physical reasoning for- 
such an effect. 

44 " 44 - 

In the cases of Pa Bk and Lw , it is found that 
the lattice effects would increase the free ion Y^ value by 
'u 45% which is comparable to the relativistic effect in this 
region, as mentioned earlier. 
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XII.2.D Some Comments on tho Watson Model for Ions in Solids 

The Watson model provides a moans to gonerato the zeroth 
order crystalline v/ave functions and makes it possible to 
estimate the effect of t'ne lattice potential on Sternheimer 
shielding-anti shielding factors. It preserves the clcctron- 
eutrality of the crystal and can be further extended to 
incorporate the proper symmetry of the ion within crystalline 
lattice. Such a model can also be used to estimate the 
changes in shielding-antishiclding factors due to partial 
compensation of charges in ionic solids. 

Tho choice of the radius of Watson sphere is arbitrary 

and it vjould be more appropriate to choose r. in some sani- 

■ ion 

empirical manner. The complete neglect of overlap, the point 
charge like potential outside the sphere and the use of similar 
potentials for toth positive and negative ions, however, remain 
the main drawbacks of the model. But keeping in mind tho fact 
that tho simple Watson model qualitatively reproduces the 
experimental trends in values it vjould be highly desirable 
to improve upon this model in order to obtain more accurate 
values of shiclding-antishielding factors for ions in solids. 



CHAPTER IV 


SternhGimer Shi el cling- Anti shielding Factors for Some Ions 
and Metals of Interest in Mossbauer and TDPAC Studies 
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CHAPTER IV 

Stemheimer Shielding-Antishielding Factors for Some Ions 
and Metals of Interest in Mo'ssbauer and TDPAC Studies* 

IV. 1 Introduction 

The experimental studies on the nuclear quadrupole coup- 
ling in solids via MJJsshauer effect and TDPAC studies have been 

largely concerned with the transition metal elements. In 

.. 57m 99 189 

Mossbauer effect studies the nuclei of Fe, Ru and Os 

occupy a central role and in order to interpret the experi- 
mental data according to eq. (1.17) it is necessari^ to have 
reliable estimates of and R. The calculations of Y„ and R 

for Fe^”^ and Fe^’*' have been carried out by several vrorkers 

29c 

using methods of varying accuracy. Burns and Wilmer have 

2+ 3 + 

reported Y as -6,31 and -9.47 for Mn and Fe . respectively 

CO 

27 

using the variational method of Das and Bersohn. Gupta and 
155 

Sen have calculated Y^ for most of the transition metal ions 

but their calculations do not include the important contribution 

174 

from the outermost d shell, Ray et al. have compared the 

The worlc presented in tliis chapter has been pxdDlished in the 
following papers? 

(i) K. D. Sen and P. T. Narasimhanh Phys, Rev. A (in press) 
(1976); (2) Phys. Rev. B (in press) (1976); (3) Phys. Rev. A, 

(in press') (1976) . 
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LQIBPT and DE calculations of Y for and found that the 

00 

two res^jlts including consistency terms agree within 12%. The 
net effect of intra- shell and inter-shell consistency terms 

175 

on is only 8 % of the first-order result due to Sternheimer. 

30 2 + 

Ingalls has calculated Y and R for Fe using the varia- 

'2n 2Hd 

tional metliod of Das and Bersohn. Preonan and Watson (FW) 

have calculated the radial contribution to R using the orbital" 

polarized unrestricted HP method. The generally employed 

less-accurate estimate of R = +0.32 for Fe“ is based on the 

30 

addition of FW estimate of radial terms to Ingalls’ calcula- 
tion of the angular terms. A more reliable value of +0.12 

7 4d 

has been given by Sternheimer. In the overlap model of 

Sharma^^®^ a major source of inaccuracy arises from the 

uncertainty in the shielding-anti shielding parameters (l-Yj^) 

and ( 1 -r) respectively. In the analysis of quadrupole splitting 

13 2 

obtained from Mossbauer effect studies on PeSiFg. 6 H 20 Ingalls 

_3 

has used a 10% covalency reduction in ^r > 3 ^^* Such changes 
in the electron density distribution of ions in solids are 
ejqpected to alter the values of (1 -7^ ) and (1-R) and it is 
of interest to make an estimate of these using some simple ^ 
assumptions. 

IV. 2 Results and Discussion 

5 

IV. 2. A Calculation of 7^ for nd Sequences 

24 ' 

Lahiri and MuJcherji have suggested that within a given 
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isoelectronic series the correlation of (or a ) values to 

q 

the position of the outermost pealc in the total radial elec- 
tronic charge density distribution, P , can be used to estimate 
the effect of distortion of the ion in going from the free 
state to V7ithin a crystal. A similar treatment of Yeo R 

in the cases of nd^ (n = 3,4,5) and 3d^ isoelectronic sequences 
has been carried out by us. Thus, in our calculations the 
following ions have been considered. 


n=3 : 

Cr'*', 

Mn^^, 

t:. 3 + 
Fe , 

Co^^ 

Mi^+, 

Cu^^: 

n=4 ; 

Mo^, 

Tc^^ 

Ru^^ 

Rh'^"^, 

Pd^+, 

Ag^+ 

n=^5 ; 


Re^'*', 

Os^"^, 

ir^ 


- 6-f 

Au 


The radius of maximura charge density has been calculated 
according to. 


m 


= I 1 X°*(j)x° (j) r? d £ 2 ^ 


(IV. 1) 


where i varies over all occupied core orbitals, including the 
outermost d orbitals, and x° gives the wave function of the 
j -th electron in the i-th orbital. 


5 

In the case of tlie Sd series v/e have not attempted such 
a correlation owing to a rather flat density distribution 
obtained in the outer region. Our results of y<a and (in a.u,) 
for these ions are given in Table IV. 1. For comparison purposes 
we have also included the earlier available calculations. Here 



102 


Table IV. i: Zeroth order ionic antishielding factor ("'^oo) 

for various ions. P gives the distance of the 

ifi 

farthest peak in the total electron density 
distribution (Eq, VI. l) 


lon 

- ( V J 

Prr, ( a* 
in 

Cr'*' 

16.67 

0.767 


(10.57)^ 



11.96 

(6.81-^,11.37^,9.19^) 

0.721 


9.64 

(6.17'^, 9. 47"*^-, 7, 97^) 

0.67 5 

Co^'^ 

8.13 

0.643 


(6.97)^ 


Ni""^ 

7.04 

0.600 

6 + 

Cu 

6.21 

0,571 

Mo'*’ 

39.36 

1.014 


(27.53)^ 


Tc^'*' 

30.62 

0.986 


(24.21)^ 


Ru^'^ 

25.71 

(21.40)^ 

0.938 

4-*- 
Rh ' 

22.27 

(19.08)^ 

0. 891 

. Pd^ + 

19.70 

0.865 

Ag*^+ 

17.68 

0.819 
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Table IV. 1 ( . . . Contd. ) 



Ion 


Pm (a.u.) - 


W‘' 77.31 

(56.85)^ 



61.67 

(e) 



52.66 

(e) 


■ _ 4 -;- 
xr 

46.10 

(39.78)^ 

(e) 


5j- 

Pt 

41.35 

(e) 


» 6+ 

Au 

37.49 

(e) 


^ See ref. 155; 

t) c 

See ref. 29c; See ref. 

175; 



e 5 

see ref. 17 4; Precise (3.^ values in the 5d sequence 

could not be obtained with HFS wave functions on account 

of their rather flat behaviour at larger r values. 
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32c 7-'- 3 + 

again, the variational calculations on Mn"" ‘ and Fe seem 
to underestimate the anti shielding effect. In these two cases 

our values are in excellent agreement with Stemheimer's 

i 7 5 -L 7 ^ 

calculations based on HF V7avs functions. Due to the 

increased binding effect as Z increases values decrease. 

The difference between the presently calculated values and 

those of Gupta and Sen is due to the contribution from the d 

shell Vwliich has not been included in the latter calculations. 

As is evident/ this contribution can by no means be neglected. 


In Table IV. 2 v;e have presented the results of the poly- 
nomial fittina according to 


= T a. p“ 
i^O ^ 


(IV. 2) 


5 5 * 

for 3d and 4d series- The average percentage error in the 
fitting is '^4%. From this correlation we have concluded that 


for a typical increase in P value by 10% the antishielding 

in 

factor is expected to increase by v 40 % for Fe^'*' and 60% for 

3+ / 

Ru respectively. 


IV. 2. B Calculation of R for 3d^ Sequence 

The calculations of R including all the direct and ^ 

exchange terms'^ have been carried out for the following 3d^ 

-f* S"!" 3 "I" 4 '4" 

series of ions: Cr, Mn , Fe , Cq , Ni . The resixlts of 

V _3 ' ... , 

individual contributions along with < r > and 3d eigen- 
value given in Table IV.3. 
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Table IV. 2; 

3 

Values of coefficients a. in |y 1 = T a. 

for Id"^ and 4d^ isoelectronic series. 

Series i 

I. 

3d^ 

0 

+ 9.94331934 


1 

+ 56.51866163 


2 

-247.487 45045 


3 

+241.13786680 

4d^ 

0 

+ 25.88997667 


1 

+254.74583419 ’ 


2 

-671.61179214 


3 

+426.21804032 
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Table IV. 3: 


The individual direct and exchange contributions and 

Rg to the Sternheimer valence shielding factor R for Cr, 

t 2+ 3+ 4->- 6 

Mn , Pe , Co and Hi ' in 3d configuration. Rg value 

have been given in parentheses in each case below Rg 

values. P and <r”^., , values are in Bohr units and the 
m 3d 

3d orbital eigenvalues are given in Rydberg units. 


o 


Perturba- 

tion 

Cr 

Mn*^ 

Fe’"^ 

Co^ + 

4+ 

Ni ' 

1 s-^dl 

0.0279 

0.0268 

0.0257 

0. 0247 

0.0237 


(-0.0008) 

(-0. 0008) 

(-0. 0008) 

(-0.0008) 

(-0.0008) 

2 S'+’d 

0.047 6 

0. 0444 

0.0416 

0.0389 

0.0367 


(-0.0127 ) 

(-0.0122) 

(-0.0118) 

(-0.0113) 

(-0.0109) 

3s.).d 

0.0236 

0. 0207 

0.0184 

0.0161 

0.0145 


(-0.0091) 

(-0.0079) 

(-0.0070) 

(-0.0060) 

(-0. 0053) 

2p-»-f 

0, 07 03 

0.0661 

0. 0624 

0.0590 

0. 0560 


(-0.0146) 

(-0. 0141) 

(-0.0137) 

(-0.0132) 

(-0.0128) 

3pL»-f 

0.037 6 

0.0335 

0.0300 

0.0271 

0.0248 


(-0.0137) 

(-0.0122) 

(-0, 0108) 

(-0.0097.) 

(-0.0037) 

3d*>-s 

-0.0090 

-0. 0035 

-0.0084 

-O. 0081 

-0.0078 

3d-»-g 

0.0167 

0.0171 

0.0169 

0.0164 

0.0158 

Rp(Ang'); 

0.2147 

0.2001 

0.1366 

0. 1741 

0.1637 

Rg(Angi : 

(-0.0509) 

(-0,047 2) 

(-0.0441) 

(-0.0410) 

(-0.0385) 

Tot. ( Ing') 

0.1630 

0.1529 

0,1425 

0,1331 

0.1252 

2p-vp 

-0.3626 

-0.3308 

-0. 3027 

-0. 2782 

-0.2569 


(0.2698) 

(0.2535) 

(0.2387) 

(0.2251) 

(0.2127) 
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Table IV. 3 (. 

. , Ccntd, ) 





Perturbat io n 

, Cr 

Mn ' 



.4+ 

Hi 

3p-».p 

0.1587 

0.1750 

0.1790 

0.1748 

0.1667 

( 

-0.1275) 

(-0.1581) 

(-0.1773) (-0.1848) 

(-0.1348) 

3d -^d 

0.1558 

0.0625 

0.0339 

0.0283 

0.0223 

(Rad) : 

-0.0431 

-0.0933 

-0.0848 

-0.0751 

-0. 067 9 

Rg (Rad ) s 

(0.1423) 

(0.0954) 

(0.0614) 

(0.0403) 

(0. 0279) 

Tot.. (Rad) 

0.0942 

-0.0021 

-0.023 4 

-0.03 48 

-0. 0400 

Net R 

0. 2530 

0.1550 

0.1191 

0.0983 

0. 085 2 


0.76733 

0.73273 

0.63735 

0.64335 

0.61227 

<r~\ - , 

3d 

2.967 471 

4.264823 

5.77 3 645 

7.49987 

9.458906 


- e. 


3 cl 


0.292534 


1.161839 2.418018 3.984175 5.32547 
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Tlie most significant contributions to R come from the 

radial excitation of np-»-p type. The exchange contribution R^, 

corresponding to these excitations is alvrays ''' 70% of the direct 

R^^ and is opposite in sign- In the case of 2p»-p excitations 

the direct contribution predominates over tlie exchange term 

and the net effect is anti shielding. As 2 increases within 

the isoelectronic series the valence electrons experience 

greater binding and a stronger overlap effect is obtained 

between 3d and p orbitals in the core. This is reflected in 

the increasing percentage of the exchange term with respect 

to the direct term corresponding to 2p-»p contributions to R 

presented in Table IV. 3. Thus, R^ (2p.^, 3d) is respectively 

74.4%, 76-6%, 73.8%, 80.9%, and 82.8% of R^ (2p-^p, 3d) for Or, 

Mn'*’, Fe^'*’, Co^"^ and Ni*^'''. Interestingly, an even larger 

effect corresponding to 3p.^ excitations can be noticed in 

Table IV. 3 , For Cr, Rg {3p.».p, 3d) is 80% of the direct term 

2 + 

and the net effect is shielding. Beyond the case of Fe ion 
in the series v^here the txro contributions almost cancel each 
other, the exchange term corresponding to 3d-»d excitation is 
zero and the direct term gives rise to shielding which decreases 
V7ith increasing Z in the series. As a result, except in the 
case of Cr atom, where due to a large 3d-»-d direct contribution 
to R one obtains a net shielding effect, the total radial 
contribution to R is anti shielding and increases in magnitude 
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with increase in Z. The total angular contribution to R is 
alv/ays shielding/ and predominates over the total radial 
contribution to R. With increasing Z the total angular 
contribution decreases and the net R is shielding and 
decreases with increasing Z. 

6 

letter the completion of this work on 3d isoelectronic 

155 

series a paper by Gupta and sen came to oxir notice. These 

2+ 3-f- 44" 

authors have calculated R for Hn , Fe , Co , and Hi as 
0.034/ 0.053, 0.062 and 0.065 respectively. We note here that 
these authors have not considered the contributions due to 
the 3d orbital and exchange tenns. Our calculations in Table 
IV. 3 show that the sum of these contributions amount to 0.121, 
0.065, 0.035/ and 0.020 for Mn', Fe^', Co^"^, and Ni^"*”, respec- 
tively. On subtracting these contributions from our total R 
given in Table IV. 2 we obtain the values 0.034, 0.054, 0,063, 

and 0.066 v/hich are in excellent agreanent with the results of 
155 

Gupta and Sen. We conclude from such a comparison that the 

exchange and 3d orbital contributions should be included in order 
to obtain reliable values of R. 

2 + ' 

Our calculated value of +0.12 for Pe based on HP S wave 
functions is in excellent agreement with Stemheimer 's"^ 

value of +0.12 obtained using analytic HP wave functions of 

177 ' ' 4' 

Clemonti. In order to estimate the effect of solid-state 

distortions of the ion on the f ield gradient due to valence 
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electrons we have expressed (1-R) < 

P according to 
m ^ 


^val 


i=l 


a. P 
1 m 


-1 



as a polynomial in 


(IV. 3 ) 


The values of a^ are presented in Table IV. 4. From this cor- 
relation we have concluded that an increase in p by 10% 

m 

brings down the valence contribution to the field gradient 
by '^47%. We note here that the experimental estimates of P^ 
can be obtained from single-crystal X-ray data and this can be 
used in correlations based on egs. (IV. 2) and (IV.3) 


IV. 2. C Reanalysis of Mdssbauer Quadrupole ^^litting Data on 
Ferric and Ferrous Compounds 

Vie have used the best value^"^^' of 1 -Y = 10.47 

CO 

3 “f" 2 + 

and 1-R = O. G8 for Fe and Fe respectively to recalculate 

5 Vn. 1 3 4r 

Q( ^Fe) following Sharma's analysis'^ of quadrupole coupling 

3+ 3 + 

data on Fe„0_, A1_0 - Fe and the T, and O, Fe sites in 
23 23 , d h 

yttrium iron garnet- These values come out to be 0.154 b, 

0.179 b, 0.146 b and 0.139 b, respectively, v/hich are to be 

Compared with the original values of O.IS b for Fe202 and 

0.204 b for Al^O^sFe^"^, We have also recalculated Q(Fe^'^^) 

from the FeSiFg.BH^O quadrupole splitting data by using 

132 -3 

1-R = 0.855 and = 12.07 corresponding to Ingalls' ^ ^3d 

value of 4.622 a.u. In order to evaluate the appropriate values 
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of Y and R v/e have exoressed these in terms of the following 

CO 

polynomials; 

4 n r" S — N 

R (orY„ ) I aj, <r~^>3^ (IV, 4) 

l'T=l 

The values of various coefficients along with the respective 
average percentage error in the various fittings have been 
given in Table IV. 4. In both the cases (R andico 5 - we have 
utilized the exact polimomial to evaluate the values of R 

_3 

and at < r > = 4.622 a. u. 

The new value comes out to be 0.156 b which is in closer 

agreanent with Q(Fe^'^'^) obtained from ferric compound data. 

134 

We note here- that in the interpretation ' of the latter no 

account has been talcen of the covalency effect although there 

is a good deal of evidence for such effects in the ferric 

57m 

systems. The present estimates of Q( ^Fe) may be compared 

177 

with the theoretical estimate of 0.16+0,02 b based on 
nuclear-model calculations. 

IV. 2. D The Results of R Factor Calculations for A Pew Other 
Ions of Interest in Mossbauer ^ectroscopy 

In Table IV. 5 we have presented the results of our 
calculations of R factors including the exchange contributions 
for Ru^'^(4d^), 0s^'''{5d^), Pr^'*’(4f^), 'IiTi^'^(4f^^) and Np^'^(5f^) 
ions, which are of interest in Mossbauer studies. The small 
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See equations ( IV. 2) and (IV.3) of text. 
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Table IV. 5: The separate direct and exchange contributions Rq 
and Rg to R for Ru^‘''(4d®), 0s^'’'(5d^), Pr^'‘'(4f^)> 
Tm^'*’(4f^^) and I'Tp^'^(5f^) ions. Rp, values appear in 
the parentheses belov; the R^ values 


Pertiarbation 

O 2+ 

Ru 

0s^+ 

■n 2 + 

Pr 

—I 3H” 

Tm 

T>T 8 + 

Np 

1s-k1 

0. 01 48 

0.0083 

0.0114 

0.0097 

0.0072 


(-0.0008) 

(-0.0005) 

(o.oooo) 

(0.0000) 

(-0. 0000) 

2s-k3. 

0.0180 

0. 0082 

0.0252 

0.0208 

0.0145 


(-0.0059) 

(-0.0027) 

(-0.0008) 

(-0.0008) 

(-0.0009) 

3 s-»d 

0.0096 

0.0044 

0.0328 

0.0249 

0.0139 


(-0.0005) 

(-0.0005) 

(-0. 0054) 

(-0.0044) 

(-0.0026) 

■ 4s4-d 

O. 0082 

0.0038 

0.0177 

0.0117 

0.0091 


(-0.0037) 

(-0.0002) 

(-0.0038) 

(-0.0023) 

(-0.0003) 

5s-*-d 


0.0042 

0.0036 

0.0028 

0,0187 



(-0, 0019) 

(-0. 0004) 

(-0.0003) 

(-0.0046) 

6s-»-d 





0. 0037 






(-0.0006) 

2p4-f 

0.0299 

0.0149 

0.0351 

0.0294 

0.0209 


(-0.0077) 

(-0.0041) 

(-0.0008) 

(-0.0008) 

(-0.0009) 

3p -*£ 

0.0164 

0.0077 

0.0537 

0.0411 

0.0235 


(-0.0011) 

(-0.0006) 

(-0.007 8) 

(-0.0065) 

(-0.0041) 

4p-»-f 

0.0123 

0.0065 

0. 0323 

0.0224 

0.0151 


(-0.0051) 

(-0.0005) 

(-0.0068) 

(-0.0046) 

(-0,0006) 

SpH^-f 


0.0055 

0. 007 2 

0. 0043 

. 0.0171 



(-0.0023) 

(-0.0009) 

(-0.0005) 

(-0.0041) 
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Table IV. 5 (. 

. , Contd. ) 


- 



Perturbation 

^ 2-1- 
Ru 





■NT 

- 

Sjp^ i- 





0.1202 






(-0.0315) 

3d.».s 

-0. 0034 

-0.0037 

-0.0165 

-0.0124 

-0.0073 


(0.0026) 

(0.0012) 

(o.ooii) 

(0.0003) 

(0.0007) 

4d->-s 

0.0021 

-0.0014 

-0.0116 

-0.00S7 

-0.0024 



(0.0006) 

(0.0029) 

(0.0022) 

(0.0004) 

5d-^ 


-O. 0001 



0.0005 






(-0.0004) 

3d-»-g 

0.0225 

0.0110 

0.0667 

0.0526 

0.0324 


(-0.0028) 

(-0.0012) 

(-0.0077) 

(-0.0065) 

(-0.0046) 

4d-*g 

(0.0059) 

0.0085 

0,0401 

0.0287 

0.0203 



(-0.0007) 

(-0.0081) 

(-0.0056) 

(-0.0005) 

5d-»g 


0.0028 



0.0160 



- 



(-0.0036) 

4f-^p 


-0.0045 



-0.0146 



(0.0003) 



(0.0021) 

4f^h 


0.0096 



0.0256 



(-0.0006) 



C-O. 0020) 

R^(Ang): 

0.1233 

0.0357 

0. 2977 

0.2273 

0.33 44 

Rg (Ang) ; 

(-0.0222) 

(-0.0137) 

(-0.0385) 

(-0.0293) 

(-0.0581) 

Tot (Ang) 

0.1011 

0. 07 20 

0.2592 

0.1980 

- 0. 27 63 
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Table IV. 5 (...Contd. ) 


Pe rturbat io n 

_ 2+ 

Ru 

^ 2+ 

Os 

•pi 2 + 

Pr 

Tm^'^ 

Mp 


-0.1076 

-0. 0430 

-0.2222 

-0.1795 

-0.1183 


(0.0847) 

(0.0299) 

(0.0190) 

(0.0138) 

(0.0175) 

SpH-p 

-0.1386 

-0.0579 

-0.5376 

-0.3961 

-0.1576 


(-0.0161) 

(+0. 0221) 

(0.2399) 

(0.1741) 

(0.0524) 

4ps.p 

0.0183 

-0.1299 

0.1903 

0.197 4 

-0.23 67 


(0.0613) 

(-0.0432) 

(-0.0769) 

(-0.0930) 

(0.0028) 

_5p-»-p 


-0.0275 

0.5130 

0.3279 

-0.1093 



(+0.1618) 

(-0.1684) 

(-0.0941) 

(0.1432) 

6p+p 





0.4390 






(-0.1943) 

3d^d 

-0.0336 

-0.0110 

-0.1462 

-0.1015 

-0.0470 


(-0. 0593) 

(-0.0189) 

(0.0736) 

(0.0542) 

(0. 03 48) 

4d-»-d 

0.0291 

-0,0374 

0.0648 

0.0573 

-0. 0580 



(-0.0438) 

(-0.0402) 

(-0.037 2) 

(-0,0036) 

Sd-^-d 


0.0281 



-0.0010 






(+0.0330) 

4f^-f 


-0.0283 



-0.0278 



(0.0146) 



(-0.0392) 

R^(Rad) 

-0.2324 

-0.3069 

-0.1879 

-0.0945 

-0. 3167 

R„(Rad) 

Cl 

(0.07 01) 

(0. 1175 ) 

(0.0470) 

(0.0228) 

(0.0466) 

Tot (Rad) 

-0.1623 

-0.1894 

-0.1409 

-0.0717 

-0.2697 

Met R 

-0. 061 2 

-0.1174 

+0.1183 

+0.1264 

+0.0066 
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net anti shielding effect of -C.0612 and -0.1174 is obtained in 

2 ~h S 2+6 

the cases of Ru (4d ) and Os "^(Sd ), respectively and these 
results are believed to be accurate to v/ithin 20%. 

The earlier calculations on rare earth ions have not 

included the exchange terms. It is clear from Table IV. 5 that 

although the individual exchange Contributions are important 

the net effect adds up to a small fraction of the total direct 

3+ 3-1- 

contribution. Thus, in the cases of Fr and Iti ‘ the R 
values become +0.1183 and +0.1264 respectively, while without 
exchange the corresponding values are +0.1098 and 0.1328 
respectively. These R values are in reasonably good agreement 
vrith the experimental estimate ' of R 't»0, 2-0. 4. 

6 - 4-1 

In the case of Np (5f ), the present calculations using 

nonrelativistic HFS v/ave functions give a net shielding of 

R = +0.0066. The only literature value available in this case 

149 

is the estimate of +0.32 obtained by Dunlap et al. from 
Mossbauer studies on (IO 2 ) and (NpO^) coordination 
compounds- These authors, hov/ever, assumed that the difference 
in the principal component of the field gradient tensor in 
the two cases is directly related to the valence electron 

. . ■ . ^ I ■ . . 

contribution due to 5f shell in Np .The relativistic 
effects for the actinides have been shown to be quite signi- 
ficant (''' 60%) in calculations and we expect an effect of - 
similar magnitude for R as well. 
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IV. 2. S Ionic Antishielding Factors in Metals and Their Relevance 
to Hossbauer and TDPAC Studies 

In recent years a mxaber of experimental studies on nuclear 
quadrupolar interactions in metals have been made using especially 
the Mossbauer and TDPAC methods. In the interpretation of 
quadrupole coupling , data on metals and dilute alloys it is more 
appropriate to include the screening effects of the conduction 
electrons in calciiLating values. Most of the calculations 
reported so far in the literature correspond to free ion cases, 
whereas, it is well ]:nov7n that in metals the core electrons are 
more similar to those in the neutral atoms with appropriate 
electronic configurations given by 3d^4s^, 4d^5s^ and 5d^6s^ 
for the three transition metal series. 

We have calculated using moment-perturbed wave functions 
for the neutral atom cores in third-, foiorth-, and fifth row 
elements in the periodic table. The transition metal elements 
have been ass’orned to be in 3d^4s^, 4d^5s^ and 5d^6s^ configura- 
tions in these three rows, respectively. VJe have also calculated 
Y contribution due to single electron (in the cases where d 

CD 

shell is occupied) by using the c coefficients (vide eqn. 11.32*) 
as ^ c (nl-^l). These results have been given in Table IV. 6. 

139 ■ ' 

Das and covrorkers have found that the use of appropriate 

neutral atom wave fxanct ions for Zn and Cd increases Y„ by 10?4 

2 + 2 + 

over the values based on Zn and Cd wave fimctions. In the 



Table IV, 6; Sternheimer antishielding factor corresponding to the closed shells and pe 
electron in the elements of third, fourth and fifth row of the periodic table 
The transition metal elements have been assumed in 3d’^4s^, 4d’^5s^ and Sd’^Ss^ 
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case of highly positive ions such as those in the intermediate 
regions in the three rows considered here, one expects a more 
significant difference. Indeed, the neutral atom values for 
Ti, V and Cr '.ro -13.4, -12.2 and -11.3 which are to be 
corapared with -S.8, -7.4 and -6.4 for Ti"^', and as 

obtained in Table III.l. Similarly, in the cases of Zr, Nb, 
and Mo, the values of -32.5, -29.7 and -27,6 are signifi- 
cantly larger in magnitude than the corresponding ionic values 
of -26, -22 and -19. 

The universal correlation obtained by Raghavan, Raufmann 

and Raghavan^ emphasizes the role of (l -Y^ ) in determining 

the extra-ionic field gradient contribution in metals and 

2 

alloys. From the available experimental data on e qQ on 

metals and dilute alloys obtained mainly from Mossbauer and 

140 

TDPAC studies, Raghvan et al. have extracted the extra- 
ionic field gradient, eq', by subtracting the anti shielding 
corrected lattice field gradient from the total field gradient 
according to 

eq’ = eq - (l - Y^ ) (IV. 5) 

where q^g^q^’s were calculated using the hxiint charge model. These 
authors have found that eq' is uiniversally proportional to 
the ionic counterpart (l “ Y^ ) and is opposite in sign. 

Such a proportionality indicates that the conduction electron 
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contribution to the field gradient in metals and dilute alloys 

cannot be adequately represented by (1-R) 

139 

coxvorhers have noted that the calculation of the field gradient 
for Cd metal based on the usual expression, 

q = <1 - Y„) + (1-R) (IV. 6) 

111 

leads to value of the nuclear quadrupole moment for Cd 

v/hich is three times larger than that derived from ionic 

17 8 

data. Patnaik, Thompson and Das have found in a detailed 
analysis of the field gradient contribution from the conduction 
electrons in Zn and Cd metals that the plane x-zave component 
is anti shielded by JgCl - ) instead of the usual factor 

of (1-r3. Now, such a prescription anphasizes all the more 
strongly the need of using more appropriate values of in 
the cases where the antishielding effect is sizeable. We, 
therefore, strongly feel that the presently calculated values 
are more appropriate to the cases of metals and alloys, and 
should be used to interpret the quadrupole coupling data in 
these systans. 



CHAPTER V 


Sternheimer Hexadecapole Antishielding Factors: 
nd-(n=3,4,5) Shell Ions, 4f-, and 5f-Shell Ions 
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CHAPTER V 

Sternheimer Hexadecapole Anti shielding Factors: 
nd-(n=3,4,5) Shell Ions, 4f-, and 5f-Shell Ions* 


V. 1 Introduction 

Nuclei vrith spin I_> 2 can possess a nonvanishing electric 

179 180 

hexadecapole moment, H. The nuclear calculations ' using 

2 

Nilsson model give K values of the order of 0.5-5b for the 

nuclei in rare earth (lanthanide) and transuranium (actinide) 

regions. Nuclear hexadecapole deformations have been studied 

1 81 

by Coulomb excitation methods for several Wolfram, lanthanide 
and actinide nuclei. However, the values of H cannot be 
always derived from such measurements. Danlcwort, Perch and 
Gebauer^^^ have measured H in the ground state of holmiim 
atom using atomic beam magnetic resonance technique. In solid 
state the interaction of H with the fourth derivative of 
electrostatic potential at the nucleus, i.e. the nuclear 
hexadecapole coupling e'^m^gH has been first measured by 
Wang“ for Sb (l = in SbBr^ via NQR. Mahler and 

■■ 'v ' 

The wo rlc presented in this chapter has been/is to be 
pxiblished in the following papers: 

■ . (l) K. D. Sen and P. T. Narasimhan, Phys. Rev. A 11, 1162 
(1975); (2) Phys. Rev. A (to be communicated); (3T~Pramana 

(to be communicated). 
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coworkers^^'^^ have detected \zltrasonically- induced hexadecapole 

115 150c 

transition of In in InAs. More recently, Dinesh and Sinith 

4 93 

have estimated the e raji for Hb in KbCl^ Xi_§: pulsed NQR 
studies. 

The effects of nuclear electric hexadecapole moments are 
generally maslced by the more dominant effects due to the nuclear 
quadrupolo moment (e^qQ/e^m^gH a- 10^ for ^^^Sb and ^^'^Sb in 
SbBr^). Sternheimer ' has shov/n for the first time that the 
hexadecapole moment induced in the closed d and f shells can 
be represented, analogous to the quadrupolar case, as 


H. 

ind 


- H 

CO 


(v.l) 


where ’^ca can be as large as 10^-10“^ for medium and heavy, atoms 

with closed d shells- For example, n,, for Cu^, Ag"^, Cs"^, In^"*", 

Ho^"^ and Hg^'^ has been calculated to be -1200, -8050, -670, 

15 Os c 

-3791, -660, and -63000 respectively. A few values '' of 

A 

e 'm. -K in solids are available and it is of interest to 
16 

calculate in the cases of several other experimentally 
iraportant ions. 

In addition to the direct hexadecapolar perturbation due 

to the nuclear hexadecapole moment (cf eq. V.l), the quadrupolar 

8c 

perturbation taken to second -order also has the transfoimia- 

tion properties resonbling the hexadecapole moment. In the 
complete interpretation of hexadecapole interaction therefore 
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it is essential to include the second-order quadrupolar effects 

The induced hexadecapole moment due to the second order quadru- 

3c 

polar effects can be written as 


ind 


H + H 
11 02 


(V.25 


arises from the electron density u^, where u^ is the first- 


order quadrupolar perturbation of u^ and is given by 



'^11 

= 

16 




CO 

with 

^11 

rr 

/ ^ 1 ' 
0 

(nl-vl) 

is 

the 

result 

bles and 

the 

sixfimation c 




2 4 


(V.3) 

(V.4) 


^ iri 

and (nf-»f)“, C values are 0.373 and 0,749 respectively. 

arises from the second-order perturbation U 2 of the 
Uq which is determined from the equation. 


(H - E ) u„ = - ( H, - E ) u^ + E Up^ 
o o 2 1 1 12 0 


(V.5) 


and is given by 


= 32 (nl-^1^-^12) ^02 


(V.6) 


where 


'^02 = ^ % ^2 
o 


(V.7) 


I 2 = ^2 ~ ±2* 
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Sternheirner^'^^ has shown that for Ho, the ratio 

a, 0. 34, which shoves that the second-order quadrtJ^Dolar 

effects are quite significant as compared to the direct hexa- 

decapolar effects due to H and therefore H. , shovild be 

end 

considered in interpreting the data on hexadocapole . interaction. 

Sell 

Sternheimer has recently calculated the valence shield- 
ing factor corresponding to the hexadecapolar polarization 

in the case of nolrai’um atom as +0.24. The shielding-corrected 

value of H for holmium atom as obtained from the atomic beam 

1 5 0c3 2 

experiments of Danlcv/ort et al. ' then becomes 1.06 b . We 
note here that and values are significantly larger than 
the corresponding values of snd R in the quadrupolar case. 

The hexadecapole interaction data on free atoms are now 
becoming available from atomic beam magnetic resonance experi- 
ment and in order to derive the nuclear hexadecapolar moment 
from such data it is essential to have the Icnovrledge of 
reliable values of Rr,. 

fl 

V. 2 Results and Discussion 

V, 2.A Sternheimer Hexadecapole Anti shielding Factors, for . 
nd-, nd^ and nd^^ Ions 

It appeared to us that nuclei with large ’ics and favourable 
natural abundance, quadrupole moment, magnetic resonance sensi- 
tivity should be more likely candidates for studying nuclear 
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Sternheimer^'^' has sliov/n that for Kb^ the ratio 

<v O. 34, which shows that the second -order quadri^jolar 

effects are quite significant as compared to the direct hexa- 

o 

docapolar effects due to H and therefore H. , slioi-ild be 

ind 

considered in interpreting the data on hexadecapole . interaction. 

Scfl 

Sternheimer has recently calculated the valence shield- 
ing factor (Ryjf) corresponding to the hexadecapolar polarization 
in the case of nolmium atom as 40.24. The shielding-corrected 

value of H for holmixm atom as obtained from the atomic beam 

150(3 2 

experiments of Dankvjort et al. ' then becomes 1.06 b . We 
note here that and values arc significantly larger than 
the corresponding values of ^oid R in the quadrupolar case. 

The hexadecapole interaction data on free atoms are now 
becoming available from atomic beam magnetic resonance experi- 
ment and in order to derive the nuclear hexadecapolar moment 
from such data it is essential , to have the knovrledge of 
reliable values of R... 

rl 

V. 2 Results and -Discussion 

V. 2.A Sternheimer Hexadecapole Antishielding Factors for 
nd-, nd^ and nd^® Ions 

It appeared to us that nuclei with large and favourable 

natural ab-ondance, quadrupole moment, magnetic resonance sensi- 
tivity should be more likely candidates for studying nuclear 
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hexadecapole interaction in ionic solids via NQR. For this 
purpose we have calculated n values for a number of positive 

CO 

ions using the moment-perturbed method of Stemheimer described 

in Chapter II. In Table V. 1, we have presented our values of 

—5 

individual radial contribution to n„, <r ^nd total 

for all the ions considered by us. 

In all tliese cases v/e have found that the nd-»d perturba- 
tions con-tribute most significantly and nf-»-f perturbations 
produce a much wealcer effect. In Fig. V. 1 X're have plotted 

vs Z for 3d, 4d, and 5d ions. It is found that in these 
three cases the value attains a maximum at Cu’*’, Ag"^, Au"^, 

respectively, v;here the nd shell is complete and most external. 

8t> 

VJith further increase in Z, again decreases. Sternheimer 

has earlier observed a similar trend in the case of 4-d ions. 

With increase in ionic charge within an isoelectronic series 
it is found that t values get drastically reduced. 

+ -f 

The a.ngular contributions to in the case of Gu , Ag , 

"f S 

and Au have been estimated frora 

( y / ) = — 

V ang'^ , ang HFS 5 

The values of y are obtained as 0.53, 0.79 and 1.15 for 

ca, ang 

' . 4. ■ ■ 4_ . ■ ■ ■ ■ 

Cu , Ag and Au respectively. The neglect of angular contri- 
bution therefore does not cause significant inaccuracy in n 

' CO . 

values reported in Table V.l. In the following discussion we 
shall assume rj — ti / radial. 

■ CO ' CO 



Table V. Is Individual radial contributions to the total nuclear hexadecapole anti 
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In Table V. 2 v/e have compiled the values of mass number, 
nuclear ^in, natural isotopic abundance, nuclear moment, 
nuclear magnetic sensitivity, nuclear quadrupole manents, and 
the presently calculated values for some favourable cases 

34 - 

up to the rare earth ions. It appears from Table V. 2 that Bi 
2 + 

and Re are the most suitable candidates for studying hexadeca- 
pole interaction. The ions In^"^, and Mn^"^ also appear to 

be prominent in this regard. 

93 

In the recent NQR studies on iTb in NbCl^, Dinesh and 
150c 

anith have estimated the hexadecapole coupling constant as 

372 kHz and 396 kHz at 2/m and m Nb sites, respectively. In the 
case of a perfectly ionic compound, niobium in NbCl^ would exist 

c V T O *9 

as a closed shell ion of Nb (3d 4s 4p ) and only 3dr>-d perturba- 
tion would contribute to the noo value. From Table V, l v/e obtain 
n., = -30.86 for Nb^’^ v/hich is a small value. In Table V.3 we 
have presented the Hoo values corresponding to Nb'*', Nb^'^, Nb^"^, 

Nb'^’^ and Nb^’^ ions, respectively. In actual case, however, 
one expects significant covalency effects in NbCl^ and it 
is therefore natural to anticipate contribution to due to 

the lov7- lying 4d orbitals in niobium. The contribution to 

' 44 - ■ 

due to tlie single 4d electron in the case of Nb 

(3d^*^4s^4p^4d^) is calculated to be -390.52 and this makes the 

total value as -421.77 which is an order of magnitude larger 
54- 

than for Nb case. Due to the significant increase xn 
value as a result of charge transfer covalency it is not 
surprising that the hexadecapole interaction in NbCl^ has been 
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150c 

detected since the values of the other parameters in 
Table V. 2 are quite favourable in this case. 

/. 2.B Stemheimer Hexadecapole Anti shielding Factors for Rare 
Earths and Actinides 

1 p o 

Extensive evidence has been obtained for the existence 
of hexadecapole deformation in the nuclei of rare earth and 
actinide region. We have therefore calculated values for 
all the tripositive rare earth (excepting Nd ) and tetra- 
positive actinide ions using nonrelativistic HFS vjave functions 
We note hare that the relativistic effects V7ill be extremely 
important in this region as has been already seen in the case 
of Y„ . In any case, the presently calculated values should 
prove to be useful in making quantitative estimates of relati- 
vistic effects x7hen the relativistic calculations become 
available. Our results of for rare earth and actinide 
series are presented in Table V. 4. As is evident fax)m these 
values, Hoo for tripositive lanthanides and tetrapositive 
actinides can be tal;en as approximately constant at -550 (+50) 
and -4500 (+500) respectively. Such a trend has already been 
observed in the case of for these ions where a value of -80 
and -100 has been arrived at by similar calculations. The 
approximate constancy of and in these regions can be 
qualitatively explained on the basis of almost similar sizes 
of the ions V7i thin each series- 
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Table V. 4: Hexadecapole antishielding factors for the 
lanthanide and actinide ions 


lon 

3 + 


Ce' 


Pr^ + 
Nd^ + 


Configuration 


4f 


4f2 


4f' 


616.61 

597.54 

- SCF Program 
does not converge 



4f '■ 

567.54 


4f" 

554.72 


4f^ 

545.31 

Gd^ 

4f’^ 

536.37 


4f3 

530.21 


4f^ 

524, 28 

Ho^’*' 

4f^° 

519.59 

Er^'*' 

4f 

516.01 


• 4fl2 

513.51 

^34- 

4f^3 

511.86 

Lu^-^ 

4f^'' ■ 

511.35 
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Table V. 4 

( . . . Contd . ) 

Ion 

Ctonf iguration 




5f^ 

4 - 4 - 

u • 

9 

5f 

4-1- 
Hp ^ 

5f^ 

4+ 

Pu 

5f ^ 


5f^ 


5f^ 


5f^ 


5f" 


5f^ 


5fl° 


5f^l 


5£^^ 

4-h 

LiW 

5f^^ 

Hypo the ti 

cal 5f^'^' 


- n^CFree) -n (crystal Ion) 

CO 

5002.27 

4753.13 7550.9 

453 4.68 

43 47.64 

4176.17 

4024.16 

3892.54 

3765.39 6904.94 

3659.01 

3553.23 

3459.47 

3377.36 

3296.44 

3194.12 5250.99 

3155.77 
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V. 2,C Effect of Crystalline Potential on HeKadecapole 
Anti shielding Factors 

We have estimated the effect of lattice potential on 
using the Watson model (see Chapt. II. 2. b) for Pa^"^, and 

Lv;* . These values corresponding to the core electrons 
have been given in Table V. 4. Similar to the case in 
calculations we have found that in the cases of positive ions 
in crystals as a result of expansion of electron density, 
values get enhanced with respect to the free ion values- The 
enhancement in values duo to crystal potential is found 
to be 'V/150%, Such a dramatic enhancement obtained in ^eo values 
for actinide ions in solids calls for a need to investigate 
the hexadecapole interaction in the cases of actinide nuclei 
in solids. 

V. 2.D Second Order Quadrupolar Effects, Ions and 

Ions in Crystal 

The additional induced hexadecapole moment (eq. V.3) 

44- 44- 44- 

has been calculated for Pa , Bk and Lw using free and 
crystalline wave functions. These values have been presented 
in Table V,5. It is evident from these values that the (5f.»-f) 
contribution to gets significantly enhanced by an order of 

magnitude as a result of the crystalline potential. We have 
not calculated the second -order contributions corresponding 
to (eq. V. 5 ) and it is difficult to comment about the 
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Table V. 5: The second order quadrupolar effect (eq. VI. 2) 

for Pa""*", and in free and crystalline 

states 


Pertur- 

bation 

Paf'^JSff) 

Free Crystal 
_ ion ^ 

Free 

_ iop 

{5_f^) 

Crystal 

Lw^"'' (5f^^,)_ 
Free Crystal 

ion 

Sd-^-d 

0.03 

0.03 

0.02 

0. 02 

0.02 

0.02 

4d-»-d 

1.08 

1.06 

0.77 

0.76 

0.58 

0.57 

Sd-Kl 

40.99 

53, 20 

28.01 

31.93 

21.00 

22.92 

4f^f 

0.18 

0.17 

0.11 

0.10 

0.07 

0.07 

5f-^f 

3.29 

221.99 

14.19 

455.92 

19.36 

285.09 

^11 

45.57 

281.45 

43.10 

488.73 

41.03 

3 08.67 


(Total) 
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Q 

total We have tested our program for the calculation 

of corresponding to (4dr»-d-*-d) polarization in Ag"*". Complete 

calculations of *^^nd rare earth ions are currently under 

progress. Our preliminary results, tiowevor, indicate that it 
is necessary to consider the second order quadrrpolar effects 
in analysing the experimentally measured hexadecapole inter- 


action. 



CHAPTER VI 


SternhGimor Shielding Factors for Atomic Quadrnpole Moments: 
np(n+l)s Configuration of C, Si, Ge, Sn and Pb 
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CHAPTER VI 

StGmlieimer Shielding Factors for Atomic Quadrupole Moments: 
np(n+l)s Configuration of C, Si, Ge, Sn and Pb 

VI. 1 Introduction 

The atomic quadrupole moment 0 is the first non-vanishing 

multipole moment for an atom. The initial measurements of the 

2 

atomic quadrupole moments of the ^2/2 states of A1 and 

In have shown that the experimental values are slightly lower 

1 2 

than the theoretical result of t< r > based on the central 

5 np 

184 1 35 

field theori’’ due to Buckingham. Sandars and Stev/art" 

have measured the atomic quadrupole moments of Ne, Ar, Kr and 

3 5 

Xe in ^2 states v/ith configuration np n's, where n is the 

principal quantim nixaber of the outermost p shell and n'=n+l. 

These authors have found that not only the experimental values 

are smaller in magnitude than the theoretical result of 
12 

— <r > but there is a change in sign of the values for Kr 

D 

1 86 

and Xe. Sandars" has suggested that the significant shield- 
ing effect is probably due to the n*s-^d type polarization of 

5 ' . 

the n‘s shell by the quadrupole moment of the np vacancy. 

ik 

The work presented in this chapter has been published in the 
following paper: . . 

K.D. Sen and P. T. Narasimhan, Phys. Rev. A (in press). 
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5 

Denoting the quadrijpole moment associated with the np 
vacancy and the n's-ti perturbation as 9 (np) and 0 (n' s-»d.) 
respectively, one can define a shielding factor Rg by. 


^"0 


0 (n ' 
0 (np) 


(VI. 1 ) 


The total quadrupole moment, ® , of the atom is now given by, 

®atom ” + 0(n’&»-d) = e(np) (l - R^) . (VI.2) 

13 

Sternheimer has formulated the perturbation tlieory for such 

a shielding effect and calculated values up to second 

order for Ne, A, Kr, and Xe in the np n‘s configuration. These 

theoretical values of Rg are in excellent agreement with the 

135 

experimental results of Sandars and Stewart which afford 
an important direct evidence on the existence and magnitude of 
Sternheimer shielding effects. 


Theoretically, it is interesting to calculate values 
for the excited state atoms with a single electron each in np 
and n's shells. In this chapter we present the R^ values for 
C, Si, Ge, Sn and Pb in the npn*s excited state configuration 
v/ith n=2 for C, n=:3 for Si, n=4 for Ge, n=5 for Sn and n =6 
for Pb. 


VI.2 Method of Calculation 

The follox>ring derivation is identical with the equations 
derived by Sternheimer in Ref. 13. 
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184 

In Buckingham's * theory of quadrupole moments, , the 

quadrupole moment, 0(np), associated with the np electron with 

magnetic quantum number m=+l is given by the expectation value 
2 

of P 2 r over the np wave function. On integration over the 
angular variables, it can be shown that 

00 

0 (np) = + ^ e / u^ r^ dr (VI. 3) 

. 5 ■' np 

o ^ 

where e is the magnitude of the electron charge and gives r 

times the radial wave function of the np electron with the 

normalization condition, 

00 

/ u^^ dr =1 (VI. 4) 

o ^ 

The repulsive electrostatic quadrupolar interaction between 
the np electron and the n's electron in Rydberg units is given by 

r 2 o ^2 

X /(oj)^ Sinead (VI. 5) 

where (r^, 6^) and (r 2 , 02 ^ coordinates of np and n's 

electrons respectively, P 2 is the Legendre polynomial for 1=2, 

and 0^ is the spherical harmonic ^ with 1=1, m=l obeying the 

normalization condition 

ir 

/ ( 0 ^) ^ sine d 0 = 1 

o 

Integrating over we get. 


(VI. 6) 
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Bj) = f V^ptr^) P2( (VI. 75 

vrhere minus the function in the square bracket of 

equation (VI. 5 ). gives the perturbation diae to the quadru- 

pole potential of the np electron vrhich polarizes the n's 
orbital to d states. The first order polarization (n's^) is 
given by; 


{ - + ^ + } ’^1 (n'si-d) = V _(r) v' (r) (VI. 8 ) 

dr"- r^ o o 1 np o 

where v' (n's-»d) and v’ (r) represent r times the radial part 

of the pertiirbed and the unperturbed v/ave functions respectively, 

V^ is the effective potential and is the energy eigenvalue 

Corresponding to the unperturbed n's shell. In eq. (VI, 8 ) r2 

has been replaced by r for convenience. The term V -E is 

o o 

6 7 c 

replaced by the Sternheirner ' local approximation. 


V -E 
o o 


^2 , 
T d V ' 

1 o 

v' ' 2' 
o dr 


(VI. 9) 


The quadrupole moment induced due to first order overlap 
density is given by. 


0 (n's^) = ^ v^ (n’s.^.d) r^dr } ("Vi.lO) 

o 

From eq. . (VT. 1 ) , the first-order shielding factor Rq 
given by. 


R 


0, 1 


=1 

5 


{ / V' (n's) v' (n’s^) 


2^ ,/ 2 

r dr }/<r > 


np 


(VI. 11) 
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We not© that eq. (VT.li) is effectively the some as the fijrst— order" 
contribution corresponding to the np^n's configuration, where 
v'(n's-»Ta) are obtained by replacing V by -V due to attrac- 
tion between the np^ vacancy and the n's electron. 


There are two types of second-order contribution to . 
The first arises from the square of the first order perturbed 
density and is given by 


f .2 2. 

J v^ r dr 


P — j. ® — . - - , 

0, 2a " ^175 2 

<r > 


(VI. 12} 


np 


The other arises from tlie second-order overlap density and 
is given by 


/ v'vi r^dr 
1 o 2 

•q _ j. Q 

0,2b " ^ 175 2 

<r > 

np 


(VI. 13) 


V7here v^ is r times the radial part of the second-order pertur- 
bation vrhich satisfies the following inhomogeneous equation. 


{ + ^ + V -E } vl(n*s-Ki) = V v’(n'&>d) (VI. 14) 

-,2,ioo2 nol 

dr r 

The total second-order contribution to R. is given by, 

o 

^0,2 ^ ^e,2a'^^e,2b ^ ^ ^ 

Finally, the total Rq is given by, - 

The perturbed Hartree-Focb- Slater (HPs) wave functions 
Corresponding to the excited state configxxration were gemrated 
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over the 521 point mesh using a modified Herman- Shi liman 
159 

program. It may be pointed out here that we have used an 

extended mesh of 521 points as compared to the usual 441 point 
mesh to satisfactorily describe the external radial v;ave 
functions for the excited states considered in this chapter. 
The integration of the difference equation corresponding to 
the inhomogeneous eq. (VI. 8) was started at the last mesh 
point r^ of the unperturbed function and v^(r^) was gradually 
varied in an iterative way till the solution became well- 
behaved near r=0. The appropriate boundary condition near 
r =0 is given by v^{r)ar as r -^-0 and v^(0)=0. In a similar 
manner, using v^' (n’s-Ki) and obtained previously, the 
second-order perturbed wave function v^Cn's-^) is obtained by 
solving eq. (Vi.14). Three to fovir iterations were sufficient 
to achieve tlie required }x)^andary condition. 

VI. 3 Results and Discussion 

In Table VI. 1 we have listed the energy eigenvalue 

corresponding-vijD n ' s orbital in Rydberg units, the expectation 
2 , 

value r over the np electron v/ave ftinction, the first-order 

shielding factor Rg the total second-order contribution 

Rq 2 , and the total atomic quadrijpole shielding factor Rq 

for all the cases considered by us- The total R value is 

9 

calculated as 1.041 for C, 1.390 for Si, 1.371 for Ge, 1.616 
for Sn and 1.703 for Pb. 
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Table 

VI. Is 

Values of the 
R ^ and other 

atomic quadrupole 
relevant quantities 

shielding 

factor 

Atom 

State 

-e„,^(ry) 

2 

<r > 

np 

Re 1 

^e, 2 

R 

0 

C 

2p3s 

0. 2686 

2.560 

1. 025 

0.016 

1.041 

Si 

3p4s 

0. 2268 

6.545 

1-345 

0,045 

1.390 

Go 

4p5s 

0. 2277 

7,067 

1.324 

0,047 

1.371 

Sn 

5p6s 

0. 2095 

9.049 

1.548 

0.068 

1 . 616 

Pb 

6p7 s 

0.2C35 

10.004 

1.627 

0.07 6 

1.703 
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In all the cases, the first-order perturbation is 
found to have a sign opposite to that of in the external 

region of r. T' ..e integral in the numerator of cq. (VI. ll) 
therefore becomes negative and the first-order shielding 
factor assumes a positive value, R 0 is always positive due 
to the v^ term in the numerator of eq. (VI. 12). The second- 
order perturbed wave function v^ has the same sign as v^ in 
the region of large r values. This makes the numerator of 
eq. (VI. 13) positive and therefore Rg is also positive. Thus, 
in contrast to the case of the vacancy-induced shielding 
effect, both the second -order terms make a positive contribu- 
tion to tlic total Rg. However, the total Rq ^ values are 
much smaller than tlio first-order term R 

0/1 

Except in the case of Ge, Rg ^ increases with increase 

in Z. The rather anomalous behaviour obtained as we go from 

Si(3p4s) to Ge(4p5s) is probably due to tiic fact that the 

binding energy of n's electron in both the cases is almost 

2 

similar whereas the values of <r > maintain the increasing 

np 

order of magnitude. Tire values of reported here are somewhat 

13 ' 

larger than those obtained earlier by Stemheimer with HPS 
wave functions on the np^(n+l')s states of rare gases and this 
may be due to the somewhat tighter binding of the np and (n+l 5 s 
electrons in the present case. 
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13 

It has he on indicated that the use of HPS wave func- 
tions to describe the unperturbed state leads to a slight 
overestimate in the calculated values of Rg. Ifev/ever, except 
for C(2p3s) xviiere is very close to unity, ' the present 
calculations predict that the total atomic quadrupole moment 
^atom atoms considered here should be smaller in 

magnitude and opposite in sign to that calcixLated from the 

1 PA 

central field theory due to Buclcingham. 
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CHAPTER VII 

Simmary and Future Work 

VII. A Summary 

The interpretation of experimental data on the nuclear 
quadrupole and hexadecapole interactions, and atomic quadrupole 
moments obtained from i-B-IR, 11 QR/ Mossbauer, TDPAC, and atomic 
beam experiments requires the knov/ledge of appropriate values 
of , R, n / ^ / which are more commonly knoxv-n as 

/ CO 

Sternheimer shielding-anti shielding factors. A survey of the 
experimental evidence in favour of Sternheimer factors in , 
this respect has been given in Chapter I. 

Most of the available theoretical values of Sternheimer 
factors pertain to the free ions and it is important to esti- 
mate the changes in the free ionic Sternheimer factors due to 
the crystal potential, charge- transfer, covalency, and relati- 
vistic effects. Further, a literature survey of the various 
theoretical calculations on Sternheimer factors revealed that 
for some ions of interest to the experimentalists even free^ 
ionic valuee of reasonable accuracy are not available. 

' In this thesis we have made an attempt to bridge the^^^^^ 
gaps mentioned above by adopting the reasonably accurate 
Sternheimer perturbation— nunerical method and HFS wave functions. 
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and calculated Stcrnhcimer factors for a variety of free atoms 
and ions, and ions in solids. We have used the Watson model 
to simulate the lattice potential around positive and negative 
ions in solids. 

In Chapter III we have presented the results of our cal- 
culations on 

(a) and values for the following free ionss 



Li^ 

3e2+ 


44 . 

C' , !■: 

,5 + 

• / 

S’ , 

Ma+, 


ai3+. 

. 4+ 
Si*\ 

p5- 

C1~, 

K+, 

Ca“'*’, 

Sc^+, 



Br", 


Sr^+, 

Yt^’*', 

Zr^'+, 

Mb' 

I", 

Cs 

Ba^+, 

La^-^; 

and 


At , 

Fr+ 

Ra^+, 

Ac^+. 




6+ 


S 

Cr^^; 

5 4 - _ 6+ 


(b) y^ and values for the following crystal ions: 

Li+, Be^+, 




1 

CM 

0 

F“ 

, Ka^, 

Hg2+ 

Al^^, Si'^"*', 


Si^“, 


s^“. 

Cl“ 



Sc^^ Ti^^ 



As^", 

Se^“, 

Br~, 


sr^t 

yt^+, Zr^^, 


Sn"^”, 

£b^". 

Te^", 


Cs"^, 

Ba^t 

,3 + 
ua ; 





At“, 

Fr'^, 


Ac^^. 



,6-f 


Cr 


Mo 


6 + 


6 + 


(c) mono and dinegative ions Y and a 


as a function of the 


radius of the Watson sphere, and expressed Y» and 

a in terms of polynomials in r7^ . 
q c . ion 



151 


^ a / ^ c values for the crystal ion of Ce^"^. 

q ^ 

For R and v/e have calculated all tlie direct and 
exchange contributions. 

(e) values for the free ions in the actinide series given 


by. 


Th^"' 

Pa^^ 

Aj- 

u-‘. 

4-i“ 4+ 

, Pu , 

, Qa ^ Blc , Cf , 

e4+^ 



4-h 4-r 

lio , Lw ^ 

U , li , U , arid Np 


(f) values for the crystal ions of Pa"*', and Lw'^'^. 


In Chapter IV we have presented. 


(a) 


Y ^55 5 

'co values ror 3d , 4d and 5d isoelectronic sequences of: 


Cr^ 

Hn^", Fe^+, 

Co^^ 

Wi^’*', Cu^% 

Mo^, 

Tc^^ Ru^+, 

Rh'^"^, 

Ag^+; 


Re^''', 

-r 4+ 
Ir , 

Pt^'*', AU^'*". 


(b) Polynomial fitting of the values for 3d^ and 4d^ sequences 

in teirms of P , vrhere P is the distance of the outermost 
m m 

maximum in tlie total radial electron density distribution 
of the ion. 


(c) R facbors for the following free ions, 

Cr, Mn"^, Fe^'*’, Co^"^, Ili^^ (all in 3d^ configuration); 

Ru^"^, Os^,^, Pr^ ', Np^'*'. 

All the direct and exchange terms in R have been computed. 

The field gradient due to valence electrons has been^^ ^ 

6 

correlated to P values in 3d series. 
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"'^oo values for the following atoms in their appropriate 

« 

electronic configurations in metals. 


!/■ 

Ca, 

Sc, 

Ti, V, Cr, I 

In, 

^4 

(D 

Co, 

Pi, 

Cu, 

Zn, 

Ga, 

Ge , 

As; 

Rb, 

Sr, 

Y, 

Zr, Hb, Mo, 

Tc, 

Ru, 

Rh, 

Pd, 

Ag, 

Cd, 

In, 

Sn, 

Sb; 

Cs, 

Ea, 

La, 

Hf, Ta, W, 

Re, 

Os, 


Ft, 

Au, 

Hg, 

Tl, 

Pb, 

Bi, 


In Chapter V vre have presented 

(a) values for the free ions of 


Ti^ + 

.1 2 T 

Mn 

, Cu'’', 

ry 2-f 

Zn , 

Ge^+ 

Rb"^, 

Zn^+, 

Nb^"^, Tc^'*', 

Ag"^, 

In^^ 

Sb^"^, 

Cs'^, 


Re^'*", 

Au'^, 

Hg^"^, Bi^"^, 


Th^'*', 

Am^+, 

Ce^^ 

Pr^ + 



, Eu^+, Gd^+, 

Tb^+, 

Dy ^ 


^ 3 + 
r Sr 


', Yb^”^ 

’, Lu^ 

- 1 - 44- 44- 

■, Th^ , Pa 

..4+ . 


_ 4+ 

* 4+ 

_ 4+ 

4+ 

_ 4 *i“ 

4+ 4+ 

u , 

NT 

Mp , 

Pu , 

An , 

Qn , 

BiC f 

Cf 

, E Md 

Aj. 

No ■ , 

and 

T 4+ 

Lw 







(b) ila, values for the crystal ions of and 

* 

(c) The second order quadrupolar contribution, to the 

total induced hexadecapole moment for the free ions and 

f 

44- 4-1- 

crystal ions of Pa' , Bk , and Lw , 

In Chapter VI we have presented our results of the 
shielding factor to the atomic quadrupole moment for the : ^ 
np(n+l)s excited states of C(n=2), Si(n=3), Ge(n=4) and Sn(n=5). 

For the positive ions our results of Steimheimer factors 
are in excellent agreement with those based on HF wave functions. 



153 


For the negative ions HF3 x\^ave fmctions are significantly more 
external. 

The free ion values for tetra-positive actinide ions 
appear to be constant at -100 (+5). A comparison with Feiock 

159 

and Johnson's" values based on relativistic HFS wave func- 
tions suggests that the relativistic effects would increase 
the magnitude of quadrupolar antishielding factors by a. 65%. 

The superimposed potential due to the Vlatson sphere 

reproduces the empirical trend that the Stemheimer factors 

in the case of positive ions in solids are larger than the 

free ionic values# vjhile for the negative ions there is a 

significant contraction. In the case of <^2 crystal ion of 

3 + 

Ce , hoviever, the free io'n values are larger in magnitude. 

The Yea values using the »Tatson model for positive ions are in 
excellent agreement with the experimental estimates for Li"*”, 

Na"^, Al^"*”, Sc^ ‘ . For the Cl case also the crystal ion calculatior 
of Y^ is in excellent agreement viith the experimental estimates. 

For Br~, and I~ a careful rechecking of the e^qaerimental 
estimates is required. 

The problQTi of estimating crystalline potential effects • 

has been attcaapted in tv/o different ways. Firstly, for the ; 

mono and dinegative ions, Y^ and a values have been expressed 

as a function of the ITatson sphere radius, r. . The calculated 

^ ion 

values of Y^ and a get significantly affected due to changes 
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in r. . If the valiis^ of r. „ can he estimated v/ithin a 
particular solid in some semiempirical manner, the correla- 
tions of "'^03 3ncl cs^ obtained in this thesis would be helpful 

in estimating and a values for ions in ionic solids in situ. 

o " ' " 

Secondly, the cin^stalline potential effects have been 
estimated for positive ions by correlating and in the 

case of experimentally important nd (n=3,4, 5) and 3d confi- 
gurations respectively to the position of the outermost peak 
in the total radial electron density distribution of the ions 
within a given isoelectronic series. ^4hen the experimental 
estimates of by means of single crystal X-ray studies become 
available these correlations would give the estimates of Yco 
and (1-R) for ions in solids. 

A large value of coupled vrith the favourable values 
of the parameters such as nuclear spin, nuclear quadrupole 
moment, ITIIR sensitivity, and natural abundance has been clio sen 
as a criterion to predict the systems with detectable hexa- 
decapolar interactions in solids. In this respect the ions of 
Bi^"^, Re^"^, Sb^"^, In^'*’, tripositive rare earths and 

tetrapositive actinides seem bo be high in the priority list. ^ 
The Co valency effects due to 4d participation in Mb-Gl bonding 
for KbClj. have been shovm to significantly enhance the t, value fcr 
niobium. A dramatic enhancement in values of *''150% for 
tetrapositive ions in solids suggests the need of examining 
these ions for hexadecapolar interactions in solids. 
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The calculation of the atomic quadrille shielding factor 
Rg for np(n+l)s excited states predicts that the measiired 
quadrupole moraent in such excited states of C, Si, Ge, Sn and 
Pb would be smaller in magnitude and opposite in sign to the 
theoretical estimate based on the central field theory. 


VII. 3 Future Work 

The future activity in the theoretical caicxilations of 
Sternheimer factors for atoms and ions is going to involve the 
use of better starting wave functions and a more complete 
perturbation procedure. The relativistic many-body perturba- 
tion theory (II3PT) calculations,, though expensive, have to be 
carried out at least in the cases of a few representative systems 
in allcali metals, VIII group elements, and rare earths and 
actinides. Tlris would ascertain the relativistic, consistency, 
and correlation effects on Sternheimer factors and help in 
evaluating the accuracy of various other coupled and uncoijpled 
calculations. 

The recently developed 'effective-operator' form of MBPT 

71 

by Lindgren and covrarkers seems to be extremely promising in 
cariy^ing out ab initio calculations of hyperfine interactions- 
owing to its local character, the HFS potential with the 
appropriate 'potential-corrections' can be used in this method 
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to reproduce tlie results based on HF potential as the starting 
point. rela.tivistic formulation of the effective-operator form of 
MBPT is highly desirable since this would enable the evaluation 
of relativistic, consistency and, correlation effects on Stem- 
heimer factors to be done simultaneously. 

In the regions where the two-electron effects on 

Sternheimer factors are shown to 3De small by MBPT calculations, 

it would be interesting to estimate the relativistic and crystal 

potential effects on shielding-antishielding factors by using 

KF wave functions. In the case of positive ions, however, we 

have shovjn that the IIPS results of shielding- anti shielding 

corresponding 

factors are in excellent agreement with the/HP calculations. 

It would be therefore more economical to perform relativistic 
KFS calculations of H and ri and compare the res\alts v/ith the 

CD 

presently calculated nonrelativistic HFS results and thus 
estimate the relativistic effects on these Sternheimer factors. 

Such a comparison has already been carried out for Yea in 
Chapter III. 

Very little is linovm about R„ and k 9 ^ even at the non- 
relativistic level of approximation. We are currently engaged 
in the calculation® of and 

crystal ions. 

With the first-order wave functions generated in the 
present thesis one can improve upon the values of Sternheimer 


for the rare earth free and 
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factors in an appjxjximate manner via the geometric -approxima- 
tion ' or more completely using the differential equation 
procedure of Ray et al. ^ Such improvements in the cases of 
more polarizable C and S ions would be welcome. Alterna- 
tively, one can use the presently calculated uncoupled perturbed 

wave functions as input to the more complete coupled calcula- 

I9b 20 

tions such as those of Kanelco ' ' and obtain the fially 
coTjpled resiilts of Sternheim.er shielding-anti shielding f actors. 

In the scamianpirical procedure adopted by Stemheimer 
to calculate R factors for high angular momentum states of 
alkalis some correlation effects are included indirectly into 
the potential. it would be interesting to see if the improve- 
ments can be achieved by using a similar method to calculate 
other Stemheimer factors. Hov;ever, as pointed out by Lindgren"^^^ 
one has to be extranely careful in such calculations since tvK? 
local potentials without ‘potenti al -corrections ' , though fitted 
to the same experimental parameter, yield significantly different 
results of R for nd excited states of alkali metal atoms. 

For several metals, the band structure calculations have 

recently become available using the renormalized atom approach. 

It is hoped that these results vjould be utilized in the calcula- 
the 

tion of/ conduction electron contribution to the field gradient. 

In such a case it ivould be more appropriate to use Stemheimer 
factors pertaining to the renormalized atoms. The calculations 
of for several renormalized atoms are currently being carried 
out by us. 
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WNLUll ,ViS21|,RSATCF(5211,l 
NFFMl|l|i£ (!7l|RLixCHi5211tl . 

?OFFCN*Vf|N(cI^B, Rsl!BftSivAL|#RUtRUEXCHfXl#X4# ^RSATCM# SNL 
iyLltFLSilLt IKaLCRe tFnUF i f IHaVlILctUCNn } #HiW#KaAiUrl 

OiFlNSXCN AC P (201 



.HUN132CCC1 

4 BSFpl»IlllJ| 

SCRTFfXi-SCPIiAl 










e 

1C 

n 

ICC 
1C2 
1C 3 


2CC2 
2CC3 
2 tie 
2C11 
SCCC 


FCBF#lCFi-i,«F7«5l 
FCBBilT III, 3E14.1tI4, 
F0flBITtF5,if3I4l 
F0BH4 J ( SF 6« c J 

F4.0,6H 

1 UH hVlLES* I4iliH 
FCPF/I1C«FS,2 ,F1C«31 
F0PHil7C!4,F5.1tF12.4l 

FOBPilTf f4|2F?-«j5I4) 
F0PPni2F4.Ci2F7«31 
FCBPATC6H 221- F4-0,8H 
1 ICH BE7/I2- F7.3) 
FCRPAlieH - 141 

FCIBP#T {6F ITEP,7XtlF 

licni ) 


FORPH"^ 

FCPPUT 




• IHZ 

14.7 

,1HZ 


E14.7I 

ZZZ- F4.0|6H Z» F4.0,11F NCORIS- I4f 
RCSFKS- 14/25H CCNTPCL CIRDS INCCRRECT. 1 

2fUiF<3.3H 

1 

IGF PC/IRCS- 141 
ZZ2- F4.0fl0H BETAl- F7.3t 

2.4?(5FC|LTiltTXf3GI-llCELI XKELl IlCllTl X 

E14,7,8XtlFZ, 13,141 
•13,14} 

.lX.lH, 13,141 
,13,141 


8CCC 

8CG1 


CCNlINtf. ^ 

R£«0 HEIOIRC C4RC. 

RliD.scce 

PRIM CCGC 

’^®RiRC®7|lJpiD!0SiRjniC 

FcSSSflli 

F^gmtniP^EifcjtnfcFS 

/) 

RERO CORTPCL CIRCS IPO 


,ILFI,IMIT,1TIP£ 


I 0»4|E14.7,4RITIC-*,|14.7i*ST-EN€tH OF 
iWITlI FDR kA1S0WI»*,l4,* !tiPE-*,I4I 

FCP ,I6,31H < USING F-F-S HIVE FtNCTIONS 1 
INPUT PCTENTIILS. CILCULITE TRIAL POTENTIAL 




12C3 


aCQ3 

€CC4 


NC4RCS-?0 

cSSIlRti 

R HSH 

1*1 

xdi»c«c 

DElljS«C.CC29. 

8g 15 

llphiliHJSlIslHpx 

lPlKit-U12C3,2€f,|4 

SHE !8: SEiif ;r=l;5Sii 

z|2***RU| 1 1) /2«C 
ZE3**"RUll I J ^ 2* G 

rLd^iS RKMC PCTEhTpt 
REIC l»JRt2(R>tl<*li437t4l 

yHp5SSpi|15JU!E/3.o. 

rgn ihft I , 


:c/3.oi 


TfcC|C^;ICR^Ie^,, 

IFfUAl.KE.l} 


GC TC 716 




RAClCW»PJIEHjS/C,52g 


716 CCKTIMI 

Z22-ICN41 
TI.CI22 » 222-*2ZZ 
OQ 224 I-ZiPESh 

R 


tlLFItlM«TtIT«RE 


224 

21C 

fiC6C 


EIC If tRM2(l}|lihiHL(XUEEUlf I«ltNCSPVS} 

REAC 8CflC|MCRIJlf J«l#KCSFVSi 
FCRPAlIZItZCIZ} 

Wfcli«C.C 
DC H2 I*liFCSFVS 

IFUESFIZ-H.C-iililt-ZZ2J-0.CCll 114,113.113 

ill Jgj5!j}j,®|hiiii|22Z,Z,FCCR|S,fiV4lES,KCSPVS 

1| iFC48SFt2E3-Z|2-Z+ZE3l-0.001l 18,18,1818 

illT FORfIt i27fI§tIr1iIg FCTEhTIILS INC Z*,F4.C,9F IN ERROR, 2F4.CI 
18 00 21 I»l,441 

Rtf! J»R03U J^l 
21 CeN|INOE 
GO 1C 211 
2C4 TiiCZ*2-»Z 

DC lie 1-1,437,4 
IK RCISi — RU2<|J*TltCZ 
Rti44|J«RL<4i7l 
Rli445)«Rll437l 

DOlll 1-1,437,4 
N«F-i 


J4RU3CI J-RIJ2III 


lice R£ll-»lJH|lic*Ro}J J + ll.0’»RUI|-i4l-RUC 


GO tC 111 



nci 


liil 

III? 

12U 


126C 


Rljinil«(ai.C*RU<n + 14*Q*BUUMl-3.0=*RU(|+Sl 1/32*0 

Rlh4|!-rllc55!j1!!+3oIoJSiil!J4l-3.0*RUl!4|!f/32lc 

coRjpyi 

RLli ?»BL3ni 

cgR“5Roi 

GO TC 21 1 
ZC2«l/2i3 

RLili-Ryllilizcz 

5«|I^^L|3S 

inSilf fzlllillcliziz 

DO 1213 I> 1 iR 

IFUisH-wf liUJ i214tl2l6 

i.iRn«R 
3Cl’« RESF 
;c»RiSH 
GO K 268 
CCMIMJI 
ica*c , 

?F ftcijll*261t262f263 
if! 1hC2J24Rt1lll 266,266,262 

I i£H’ «| 

i VtlS— TltC222/RUl 
GO to 26 

1 vni;BU|ll/f Cl» 
i Cly*iT|Ryi 

I 265,265,266 

' tpUisJ-R] 261,268,267 




2€1 COMINUI 

ns vfiU-Tiiclll/Rlii 


26S CGMlRUf 

Dfll* -ICCCCCC. 
NltEB-C 
NCM'Cli -3 
STIRT XTSPAIICft 
PRIB"; 9CS 


21 

4C2 

5CA 


" li» ii K I# S ** 

PRlNt 4lc 


r- 


MtER} 40Sf2a«28 


EC iC4 


C 

C 


Rll|l«BU|l ^ 

iBl3CII#RElRlllIItRUFRLHntRUlNL2CI»fRUFRL2<Iif BU 

/111 

Jic F8Bll*!^l4CHlIt3<FFE3fREIHl,REFhLltRUINI.2#RGFNL2tBU I 
491 PQRRill aSfF10*4f Ie14«7} 

GO TC 2C5 
28 DC 2S |-1|PESH 

Rlltil- C.O 

R|CCR|(J:j>C.C 

SCLvl“sCHRciclRG|R ECG4TICR FaR_liCH STATE IN TURN. ALSO CALCULATE 


3C DC 


SCHRCiCXNGER ECLATXCN FOR EACH S 
ctRE are vherce charge CENSXTY. 

9 H»ltRCSPVS 
-_(R 1 
NN-RftL2fPl/lGO 
LAP »NllL2(r>l/XO -10»*RR 

lap 

CALL SCH|C*l,EjLAP|RRtPRKtPISHtC#THRESH 
XF?P-NCCRES}30t3Ct32 

r>rk 9 1 * ti it 1 


31 




RSCO 
GO TC 34 






SCI 


SNlfSi|}-lHC(II 


37 


NKf(KfR}«KKK 
MClRCS-RCIPCS+a-HtUKK-llMOMe 
33 EE lf‘I»E 

CAlCllAtf ICTJil CH#PGE CE^SITY JlNC ATCRIC EXCHANGE POTENTIAL 
DO 37 I®ltPESH 

36 RSATCM|J*BSCCBECIHRSVILEII| 

Rl|XCHlil»-6.«nFA=i{C3.0*P(:)=iPSATOPlIll/315,a2734 }♦♦( 10/3. 
|i|ATE nCRiC CCLLCPfi ECTERTIAL 

C 


. 'C^C 
, JlP«C«C 
H»C.CC2S*C 
I»1 

XKlloQ.O 
00.38 K-liAC 


m tk % m, . 

m- 

SSifi 

A2-RSATCN(II/2.C 

82«RlAtCNfl}/(2«G4R{n} 

xllf}M|UP-»Al*H 
XUai-BlOP + BKH 
Ai»il-»A 


38 


ASLR>X] 


I 


a\ 


B|LR*>jni 
a!«A2 
B1»B2 
H«H+I 




Cl 





712 

714 


,1>6C 7C 717 

^IFIFllflcifpiCICM CC TC 715 
GG 1C 714 

3iJUI«ll4<:i47li3«IC^ 

CCNIULfi 

ocSl ^ V '^^' 

RIIMK 


llii-RiiRLZfl} 


C 

C 


42 

4i 

47 

46 

5 C 


RGfNLlCil"RlFNL2<l} 

Rl|NL 2 dr-Ptl|l 

RLFNL 2 U;>«X«(:} 

N|lER«R!tiF-ll 


51 

8 C 

61 

270 

271 


If s!f CRillpiChRCf^SATlSFlEC, CALCULATE NEXT TRIAL POTiNTlAL. 
PARl^CiT 

Ifunai, 

PiiRTi»l.. 

DO 


I.NE.IJ PARl-0*5 
.-PARI 


IF 

w 


P-IJINiSFJ _ , 


RlllG «<RL|P-BUIlIPniJ/IRU2P-XjlLIMlTII 
00 2<1 I«L!pniNE|H 
RLCIl*RU2P-iiT!cl1PLZP-XJCll I 


LIPS 


■PESF 



€2 CCMIM1 


C 

C 


< - 
le 


H 


IF DELTA IS ^CT >CACTCA|c C 


18A 

189 


FRA1T IFFPCVEPEM SChlRE 
CW-1 . 
tfiO 1 185 


CECPEASING FODB TiMESt BYPASS 


182 

83 

832 

84 

841 

842 

\h 

1845 


613 

Ifli 


r r» ■ i * f- r r* », ¥iu 

NOBPCB -NCBPCN-l 
IFlNCBFC^ieCieC 
ALFH*C«5 
PBA|] IPPPCYEffNI SCHEPE 
DC 845 I*2 <ICDi 

XNLPC|l«BLiMl«II*BLF8L2fII-FU2NL2(| ^♦RDFNL1^ 11 

f SBBlIl-RlFAyai-PLFAllI n-Fu!NL2l! I+RIJINLK II 
?IAilFI0|hP Ul/PHM2(in-0.CC01J 83f83f84 
CQRTlPui 
ALPH-C.B 
GC tC 848 

AlPH»iANLPn»/DENP(IHPUFRL2nH/SNtO{U 
IFfALPHl E41*845t842 
ALPH-C.C 
60 1C 846 

iFfC.S-ALFHl 843f845t845 

lPRSh*IPR58"' 1 
IFIPM) 1813il84Sfl873 



COPT! NOE 

ASlN«)l!fl)«IU3HXl(4HXll5} 



DC 814 I«3ilCCl 

A|{flJst!i-S(5li-2 HXIIl + 31 

IFtglTI Iei3,l814tl8t5 
CCMIhUE 
ICl»!C+i 

DC 815 f-ICltKESF 
KJ(|)«C.C 
R(j(U*RLFKL2 <I) 

Rcdj* RUF^lSClHXJdMSMCl II 

mmr2nt,uiu2Bio 

2878 SCll-FESH^ 

hrilViH^: 


€74 

1874 


175 

1875 

€7e 

87 


2871 

187C 

851 


854 

nil 


_ 2iFESH 

vlff-RclIl/P U I 

CCMlRCi 
IClt -C 

- cCt»i 

V^I|~TfCU2/aa. 
vai«Ru 
ScIflMUL 

xlfll«0«0 

nbi 7 b1ai eigen v#iLis p 

NbRES«5G 

gSlhaSlillSlWi 

H*C*GG2S4C 
ASI.N«C«C 


-X14ST 



c 

c 

c 


53 

54 

ss 


5€ 

sic 


i-i+i 

A2-Xi(l}«SI«L If fIH^2 
A1-A14A2*I‘ 

ASl»<-ISlP4/H- IA2/2.C HI' 

H*H4H 

Al«IA2/2,CMI- 
E|Cfl«illKMISil' 
NCIRCS*^CAPCS48*I«42 
GC 1C 21 

RfStiLTS TPAR5FEFPED FPCF 
OllPl 
CCF 

CCMIMJE 

PRIM l|fPC4POS|FCIRES 
PRIM 11 % flRCij 


iRTfPNRL RERORY TC OUTPUT TAPEfSI. 


71S 


581 

582 

51 
51 


588 


M 1U ijRCl >tI*ltRESH 
PRINT iH. I BSj»TCRIilf|»lfRESH 
PRINT H?l iflU3(!lt!-l.RESP 
FCRRAT|13i,10E12«5) 


1 


RESPl 

. rRPATIl]l«lGE12.5) 

NC«1 

PRINT 2CC3|2fN€0R£StNVAlESfl0NtlZtNC 
DO 583 I“l*4Al 

lF(TVCICN4HiNL2lin S83t583t58l 
DO 512 N-IiA41, 
r|IN|2|M*-TIiCICN 

ccn^IruI 

CGNTlNyi , 

CC 888 PIN «lt44C,S 
HA]i«>HN44 
NC«NC41 

CCNTlNUl 

NC«NC41 

DC 5 % H*ltNC5PYS 

NL2«NNL2fNl 

Ki(R«NRKR(R) 

XL»NLZ/10-2C4fNl2/lC0l 



SSI 


8CC8 


f J 

MELCCK-KKK/4C 
LP*)fl + I.C 

iriki uu, 

*.i5Ui:Siioi 

Aaf5i«>SHl (PfI-*lI/Pa + lH*Lf 
CAtl CRCSIP Hi 

PRlH?*2CC2,RLZ,XUEE(P|t>(h5Llf ItKfKf #I||5I. 

WRITE i1fE*C5» * 

FCRPATIll" il2.1Kfl2tlM tF4«Ctl3til4.1} 


IZfNC 


Kl-gPK-l, . , . 

DC 5S8 RlA*liKlt5 
NC»»iC4l ^ 

mIji- i'lfi+4 „ „ 

PRINT iiAVE FUh ET RSPCVINC IT 
5SE CCNTUUE 

^^"pRINT hA\£ FUA EK REPCVING IT 
DEni)i**CC2S 
i ■ 

4»1|FBICCK 

4 J>4444C 

fFI4-II 8C5C,8CS1,8C50 

ecsi N|||*41 

60 J5 8C52 

I ^^kRlIl Hi8C531 N6hfC|LTA3(t)'f44n 

3 FORNATllH fiafiFS.^I 



8C 


OElTAJi- 


Ifu-ii less, 8CS^, Seas 

8C84 WRITE 54,IC5jl IS?UP.ntI-lf44l 
8ci€ F0 Sp 1T(1H taS14*7J 

8C3S WRITE fSiicflEI t lSM,tPt I * 1 1* JJlt jJZf 2I#SNLIH# 443 1 1 





lie 

\i\ ki j “ 



J41 



LISTHTG OP THE MAIM COMPUTER PROGRAM 
TO CAIOJIATE AMD R FACTORS 

03 



non non nnnn non non 


D1><£KS1CH XRlfl3}tXL»>r>Uf4)flCi4}tMENTRY{ 13I.A0(;eiTU1}iN£tECTIl7liX 
/LMLNM?}t^CFe6(17lf6UKC(13t4l}>KEMTRYU3ltLENTRY( ll lfDei.KRT ISIrOtt 
/XRRn|}«€ARXlMa7.4}tRR(17t4|fCRBlTPf nt^IfNEXlTAf 11},|NcRGVC|?l 
DIREKSICN liFL{13f47).XR{13t4llfNUNfl6tll»A4)tFUNC(13t47}tA{4,5} 
COFRCR liP|i|XR,MUAtF(;k»A 
CCFfCA ipRG^ 

CIPENSICR SFDStl'<17f4),SP£SU«<17f4 J 
CAIL FLUNC32CCCJ 
RRloQ.O 

CALL FLCV( 32000) 
nAPi-4 

C READ NCAlJpf ^HE NUPEER CF THE ATCP/ICN TO 8E TREATED IN ONE RUN 
REhlAC 4 
READ SOCfPCAICP 
SCO FORNATfli) 

THIS DC-LCCP fEPFORNS THE CALCULATION OF ALL THE NOATOM ATOMS ONE BY 0 
DO 7 NAT0P»1|NCATCP 

READ AND PRINT THE NAPE CF THE ATCM/ION BEING TAKEN FOR THE CALCULATIO 
READ llTAPiieOCl) 

8Cei F0RPAmX|3?H J 

PRINT gooc 
SCCC FQRPATC//) 

PRINT 8001 
PRINT 9001 

9C01 FORMAT {125H 4«*4 4 4 

/«««« ) 

PRINT 9001 

READ Z. THE SCALE FACTOR USED IN DEFINING THE VARIABLE -XR- AS XR*Z*R 
H-F-S CALCULATICNS* 

.C« ^ 

ZZ-{0«l8534x38/Z)4A3 

READ NG0R8# THE TOTAL NUPEER CF THE ORBITALS IN THE ATOM UNCIR CCNSIOE 
READ inTAPEj8C04) NCCRE 

jtape*3 

NRITg I JIaIi. 80041 NCCRE 
PR2NT iCC4fNQCR8 
8C04 FORMAttlH ,131 

^ THIS DfHLCOP OCES the CAICULATICNS FOR THE VARIOUS ORBITALS ONE BY ONE 



c 

c 

c 

c 


kHEh OClhC CALCUL/ITICKS 
hCRBVA »3 

THE FCLLClilHG C4RC 
DO 5 NCRB«l|hCOPe 


SUFT 781 TO 789 AFTER gSC • REHOVfc 


BgAD ACRBITt THt 
RtlD HILECT, THj 
R|#D XlhlKf TFE 

rIbd ncreg* tf| 

THE UNPtRTURBEO 


CRBITBL 

HupeER m 


PICTRASCCPIC DE 

*HE electrons 


L VRILE^FCB. Tli|^|PBIT/»l 


I GMT ION 
N THE ORBITAL 


NOHBER CF THE 
hAVE FUNCTION. 


IONS USED IN THE TABULAR REPRESENTATI 


READ flTAFS.eOOS) ACRBlTlNCREIfNELECTINCRBltXLHUNfNCRBI.NOREGINORB 
U.iNERGVfNCRBI 

NRITE UTAPE.80C5) ACR6IT ( NCRB ) t NELECT ( NORB )t XLWUN CNORB )t 
INgRECtNCRBifENERGYCNCRei 

PRINT 8C05f ACRBITINCPBJ ,NELECT«NOfiE»tXLMUNCNORBJ, 

IHOREG INORB I. energy CNCREJ 
8CC5 FOBNATIIH t A2 .IH « ,12, IH » , F4. C t lit S 14.7 » 

HORE€*NCBEG INCPB > 



RYCNOR6I .THE LNFUBTLREfcC NAVI FUNCTION CF THE ORBITAL 

NENTRYtNCRGI.OELXRtNORGl.XRICNORG I 


NCRG-: 

ilTAPI 


fNCREG 

8C53} 


DO 

Rfnu t & iArc t e 
SCSI FORMTCIH •!3t2F9.4} 

NENTRY-HENTRYINCRGI 

READ fITAFE|i€|6} |MUN{NCR8tNCRG« J»1>NENTRY| 

^ |hUMNCPBtN0R6.4l iJ-1.40> 

(NUNfNCRBfNQRGt 4} .4»1.40} 


C 

C51 

C 

C 

C52 

C 

C5I 


IF INCRG-1 1 51.51.52 
GC TO S3 

MRIIg {4JAFEf80S6) 
PRINT 80S6j 
:ONTlNi 


8C58 FQRNATUH .SllA.T} 

NIIH the help CF XRUNCRGl 
REGION UNOIr CCNSIDERATICN 


DEIXRINCRGI-IJOELXPINCRGI 
XllNCRG.l)«XRl{NCRGI«Z 
DO 5 4-2.NENTRY 

XRiNCRG.4)«XRfNCRGt J«1)4DELXRINQRG} 
NGRSGoNGREGtNCRBVAl 

F|NCRG“Il 56.541.57 


{ HUN fNORe.NOR6»4 ).4«4»43 ) 
CHUN fNQRBtNORG.4 ). 4"4«43 I 


ANC OELXRINORG} HE NOH GENERATE THE XR MESH 
CNCRG-TH REGION} 



C 

C51 

£ 

C55 




Lt 5 pE , 80 S 6 I 

PRINT 8C5«, 

CCNTlftCf 
READ 2224iNCAL 
2 224 FORNATtpi 

- - gS TC (|226x2225t2225}«NCAL 

2225 r| 40 ^q^.KeRBvS 
sis FORNATtth 
10*0 

N0RE6-NCREGINCRB%AJ 
U-C.C 


pRCNORG , J ) yJ » l « 4 QI 

bPCNORGy J )» J - lt 4 Q } 

{ XR { NCFe » J |, J « 4 t 43 } 

tXRfNOFCtaHJ -^^ t ^ a } 



NCRG»ltNCPE(> 

NENTPTCNCRGl 

IFfNCRG-1} SSS4,SSS5fSS94 
SSS5 Jl«2 

GO TC ssse 

SSS4 Jl-1 

SSSfi DO ice4 J»41tNENTB¥ 

1CC4 FUNCINCRGtJJ «ihUh«NCRE¥AfNCB€tJl*^21/lXRfNCRGt4l’»*3 I 
CAP INTfGR {SXPVA,KCREG,PENTRY,LOyOEtXR| 

PRINT ICll.EXPVA 

ICll FORNATI/zIIf TNg EXFECTATICN VALUE C NUN=»1 1/IR#*3 1 l^HUN I FOR THE V& 

/lInce crbItal-.Sh.ti 

2226 SURRAQ»0*C 
SLRAN6«C,C 
SP|“!CT»C.C 
SPD1CT-0.C 
sr|ang«o«c 

SRCANG«Q«C 

srIrad-g.c 

|rcrao«c«c 

R"SLF«C«0 
Rl«0.0 
SU6RAD*0t 
SL6ARG*C*i, 

DO ir NCRB«ItN00BB 
NOREG-NCRlGfNCRe} 

E2ERC«ABS|ilt|RG¥<NCR6}} 

LNlNxXLhUNlNCRB ) 

: PRINT THE ORBITAL SPECIFICATICNS 


PRIN ^ 





onn no non 


-LESS INC THE^ ICE GIVEN IBfilTRlLLY. 

REIO ^*l«E>ni*NCPEI fNCRGN.JN 
MR|TEUTAFEfl) N|X|fl C NCRBI tNCPGNt J« 

PR|n 1 4t NEXillCNCRBItNCRGNf JN 

k F0PnA|(!2|2I3) 

NEI-NIXITIINGRBI 

IF THERE IS KC EXITITICN TC BE CCNSIOERfcCt SKIP TO THt NEXT ORBITAL 

1F<N|XI SClfTCCtSCl 

ICCPRiNiiCl 

Kl FaRNATC//HH NC EXCITATICN// 1 
n|XI 2AINCRBI“1 
RR(NCRB*1}>C«C 
GAMlNlNCPBi U-C.O 
GO TC 7T 

5C1 con«!nu| 

FLNCUjlI-C-C 
00 IC03 NCRG»I,hCREC 
N|NTRT«N|NTRYINCRGJ 
IFINCRG-ll SS9lfS992t999l 
9992 Jl-2 

GO TC 9993 
9991 Jl*l 


DO 1CC3 J»JlfNENTRY 
1CC3 FLNCINCRGiJl«liiUf»{NCREthCR6t J»**21/IXR<NCRGtJI**3l 
CALL INTEGB CEXPliUNfNCPI6,KlNTRyfL0f CELXRI 

1C12 FQRHATf//|fH THE EXPECTATICN VALUE I WUN=*C 1/IR**3 I l•^^UN » FOR THIS 
/RiI7AL?»ll4.TI 


AlI»l)«l*C 

! Aflfll-tuMilRBf 


1CG5 An;51-RUMRCRBtltl4l}/CXRfltl4i}**{LMUN4|i| 
CAlt 5IKULT <41 


1CC9 


CONSTANT A*t£I4«7//» 


READ XLRPUt THE L VALUES INC C* THE C-fLHUNtLHPUl VALUES FOR ALL THf 
EXITAT:CKS IN the QRBlTIL 

wllT|t5TS?i5lf loSiffHNiRfifNlxiViJSiJulNlxllflilNixlllJfiEi^^^^ 

Vr!nT 6« {0R8ITP(NER6tNEXlT}»XLNPUINEXlTlvClNEXlTltNEXlT»lt 



uuuu 


'fSP* 


KATf4fA2i2F8*Sn 

THIS CC-ICCP fEPFCRHS THE CALCUL«TIGNS FOR ALL THE EXITITIONS IN THfc 0 
CHE BY ORE# 

00 777 Hi3(IT?IaH|X 

(NOReivLMUHfCReiTPfNCRetHEXITtfLNPU 


ww iff 

PRiHT^IafiAclBlrlhui^c !•!.«' 
BB JJIxIt J-0-0 

|RESUBfHCRetBEXlT}«0.< 


222 


S 

221 


8 

22Q 

16 


;7/16F 'THf EXClfil 


'li:APif222l ACRI 


Cff f #2 ffSH IL*f 1 2f 5r I TO 


2flH»// 


rSU! 

/) 

DLfcUhP»(LliPL>*(LHPU+U-LHlJH*UliUH + lll 
IFJDLtiOkPI Sf8t9 
PR|hT 221 

F0R»'AT«29H TUS IS # ANGULAR iXITATICN //» 

EXTEST-C.C 
GO TC 16 
PRINT 22 C 

FaiBATl28H THIS IS A RACIAL EXITATION //I 
EXTEST*l»C 

CALL BURMA V(LNPUtLMLN»EXPl)UNfNCRGHtJBfNCBEGtHENTRYtNQR8tC£LXRtZZf 
lAllRCr 




1C4C 00 11 

T 

11 FUNCI 

CALL ^ 
PRINT 9CI 


|04Cf 1041 

ill 


RY*B|NtPY|NCI 

CVEBLAP INTEGRAL IS 


*tE14.7l 


9C11 FaBHATI __ 

DO 12 J«ltNENTBY 

12 MPLfNCRGtJ|«>hPUtNCRGtJ)*-CLAP4WUN{N0RetN0BGt4} 

LET US NCM CHECK THE ORTHGGCNALITY BETMEEN HPU AND WUN. 


DO 9CQ8 NCRG«lfNCRfG 
NiNTRI-NiNTBTiRSBGI 
DO 9CQi I»ith 

,cc, 

prIht 90 
9tCf FOiNATI 


rCl 


NTRY 

Al7kcREGtNEfitl!Y,L< 


1 1 « MON fNQRetNGRGyi I 


RGtl} 


THE OVERLAP INTEGRAL IS NON »tE14*7} 



nevTi rtnrtn 


NOh PRINT THE fUPTUREED hRVE FUNCTION MHICH IS OBTHOGONIL TO THfc iNPUR 
ONE 

ISI^ FoISItIz/SCF THE PUPTUREEC W/Vi FUNCTION FOB THfc EXCITATION IS//I 
ICIC FOBNATCIH fSEH«7} 

NOli EVAilATE THE INTEGRAL t NPUA < F>*A2 I AHUN I BEGUIRfcD TO CALCULATE GANA- 


2 227 
2 228 


GO TC {2228f2227i2228},NCAL 
GAPAlNINCBBiNEXni-C.O 


2222 

2223 

6861 

€U6 


is I| N?il»ltNOB|G 
NiNTBY»NENTR1lfNClG) 

0 io ']L ^ kJI ^ JL Pi 1 ^ IN If ^ 

FUNCfNQRGialakiUNiNCReiNCBGt J}AMPU(NOB 6 t JlAfXBt^ 

CALL INT|2r IGANAfNCRiG,N|NTPY,LO,CELXftr 

GAPAINCNCPB.NEXnf-CCNEXITIAGANA 

GNl»GNi-»GAPAlNfNCReiNEXITI 

IFINCAi-ll 2223,222*t2223 

RRfNCRB|NEXI1}«0«0 

GO TC 1111 

IFtNCRB-NCRIYA} £666f6661 i66E£ 

"giSHf 

nIIS Glvis’’TCTAL NG CF CIFFEPENT CONTRIBUTIONS 


Jl*IXRtNQBG,Jl*«2J 

I 


:|fi Glvis’’fcTAL 


CALCULATED 


FOBNAtilXjaSfFS.S) 

00 A72 NT1p|-1sNG|n 

'^i sKll i! ni-“ 

KflR|6V«NCRE€fNCREVA} 
NQRSoNCRB 

N0R|G«NINUNCltNC2} 
NORloNOReYA 
N0B6-NGRB 

NOREGY-NOEG 

oS^ill NSli|»ltNCR|i 


fNCBli 

PEVAl 

CltNC2} 


lSfS|Rl5-l|p|lrl86!i81t887 

N|NTRY*N|NTRY-2 
GO to ifil 
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8^6 




613 

8 ^C 


II 


89 «« 8 S 5 i 8 gi 


867 

JRl-5 
DC 884 dl 
IFIJf3-3l 
KC*2 
FUTl-O.C 
FI».t2-0,0 

IF(JR-l1 885 |876,88i 
DO in NCRG-itNCPGR 
Ki8fBVCf»CBGl»l5EMRV(r.CP6l 
Kl81fil<8CPCP 1*JR42 
IFUR-NSN1PV411 e72t873t873 
XF(NCRGR-I^CPEG} 890,872»872 

V<Ptf8CRGR,RERTRV-illi»|WFCjChCRCP-»lt4l+WPUIRCRGR+lf 5IM0.5 . ^ . .. 

WL|iil8CRB#RCPGRfNENTPY+l)*l KUMN0R8fNCRGR4lf4 |+HUNtNOR8*NC]RGR-Hf & I ) 

^IiJlj 5 iRBV 4 *RCRGR,KEMRY+ll»thURCNCRBV 4 tKCRGR + l« 4 »+WDNCNOR 0 VAt 
lN0RGR4lt5l!’*iC.5i 

XRINCRGR#8SbiRY4U» lXR(NCBGB4l,4HXRINCR€R'»l,5 n*CC.5> 

' ii “ ‘ 


812 


418 

814 


tcSt'*f» 5 lBV'» 21 »liPU(WCP 6 F 4 l. 5 l 
R6t8CPGRiNEK7PY42i«WURIR0R8tNORGR4l,S} 

PCRGP|REMBY42)«RCMNCBBV4#NCP6R4 1,5> 
PlBY42|«XBCPCBGP4l,5» 


SEiiJDg 

MiiKfRORBVA 
XRIIiCRGRtB, 

DO 814 ^CPG«liNCPGP 
NJ«R|RTRYl«CPG» 

» f* i E* # 81 ^ C C' 70 '^7 € 

FURC|80RGf4l*«ONINCR5f NCRGf4l>»«PUIN0RG# JI^IXRfNORGf I 


GC 1C 
FCRCIR 


l}sM(JM(RCP5tPCRGt4MHPU{NORGt4} 
'*'^i|NlitPCRGP,K|NTRYfLO#DH.XR> 



994 
9 93 

816 

an 


;s-wcBGP 

il-iilsG-PCRGS+l 

,4R-RiN?RY4l| 994tgg4t9g3 
R *3 

TC 817 
R*4 


i 

kluii 


N«|iGRGS«l 
fORGl»K 


CP 


NOI ^ 
RR»jR 

Q 818 


GR 6 G- 6 CRGS 41 

R8«1«8CRGL 





SS2 

«€3 

8 S 2 

881 

6 84 

610 

UC 7 


1CC6 

CC6 

S 8 


OS 6 

095 

891 


876 LEIi1RYCftN|-F6h767*»ihi 

LptR^UI-fEMBYiKCfiGSI-JPP-H+^Q 

00 861 fiN«ltKCRGL 

N»P+1 

DilXBRCPRI-EglXBiNI 
NJ-LfKTBYC«B I 
ifinh-lY 683 |8a2f883 
J-JRP-I-KC 
GO 7C 892 
J«C 

DO 881 

^J| jQGB *^1"' J 

^ ' FLRCiRP.JJl-WtNfPCBStRtJMhPUCN.J }/(XPihf JI**(L4 JIJ 
CALL INTEGR (FlM2|PCPGC.L|NTRY,K0tDlLXPBl 
IF«4B-K|MPY-»1| 884,991,884 

Fl«CliPLfPCPCP,JPMhLhCACB5.PCP6P,JRJl/l <XRfN0RGR.JR|4*u*lilJ 
F2«( fcPL(«CP€F,JR-»2l=*»ilJR««CP5,RGBGB,4R + 21l/ltXRIN0R6B, JR*2I*^IL + 1 

FIM2-F1N72- CfFl4F2HDELXB|hCFGRI>iCl«/2. n 
GLKCCNCPGP,JPI«CFIMl4FIRt2<CXBIN0RGR,4RH*i2*L-il ! H 
GURCII«CPGP4!,n«GUNCCRCB€R,NfRTBY-6J 
G0RCfRCBGB-il,2l*GljRCCRCRGR,NpTRY-4l 
GliPC pCRGP+1 ,3l*GtPCCPCPGP,HlhTRY-2l 


17 


;LPCCRCR6P41,4}«GLRCCRCPGP,(i|£RTRY) 

;g»itirue 

iliRCfRqRfG,A|NTRY+lJ«GlJHC|RCP|G,NENTRY) 

GU8CfR0REG,plNTPU2}«GiJAC{ACPlG,NENTRY} 


NgP 61 «NGREG 4 l 

DO 1C08 NCPG>NGRG1,RGREGV 

NiRTPY-PfPTPYlWCPGf 

DO 1CC6 J»l,8fSRTPY 

GL|CfRCRG,4J»6LltIT 

PRINI 998 

F0R8ATt//41l' THE GEhERAL CUACRUPOLE POLARIZABILITY ISV/I 

FyACIl,U«C«C 

00 891 RCBG*1,R£P|6V 

NEhlPY-PEPTPYtHQPG) 

IFINCRG-ll 996,997,996 




JO fc 
41-1 


995 


DO 891 4-41,PEN|PY 

FyAC{R0P6,4]«&liRC(KGPG,4}«{l«URtNaRBVA,li0RGf Ji4HUNINQR6tNQRG,4ll/ 
l{XRlRCR6,4f*«{L4U| 


OACi:iRTiGR_ |R,NCRE6Y,PIHTRY,L0tDEUR» 





IFINllK 
RDiLI<«S 
rIsir-sr 

PRIHI 416 



2^1 

1CC2 


ICC I 
ICGC 

nil 


SRCStRCRCREfN 
SRgSLR<RCRBfM 
fRCREtREXin 
_^R|«CP 6 ihiX|TI 


HR! 


?(|T)«SRCSURINOBBtNfX|T 
x!n»SRfSyRfNOR0tNfx!T HRl 


tRCSUMfRfSyR.RRtNQRBtNgXITl 
IT HSRESURlKCRBtNEXlT J 


1 


RRfRCRBtN - ^ - - - - - ^ - 

FoSS!lil!<t*RtlRE»<fI3f>*L •♦tI3t*C0EF-*,F8#5#*R FOR 
l*R-DIRECT»*iil4,1i5R-ixCH-*,|14.7t*R-TOT#L**tiH.7 
PRINT 1C02 

FGRRilTi//2eF THE RAIIC GAMiK F IS//I 

DC ICCO NCBG«lfRCRECV 
NiMR 7 »RERtFnNCPGi 
DC ICCl J«lfREMRY 

GtNCIRCBGi4)«GyhC(NCRGf Jl/GARA 
CCRTJME 

PRINT 444fRR{hCRE,hEXIT}tCAPA«GAMAIN(NCRe»N£XlT} 
NCRIG-PCPEGtNCREl 
NCRBl- 11 

lFfNCRB.Nf.ll } GC TC 324 

liRlIlE i 3f32l) {NCR|G,NCReit 
PCNOi 321> (NCR|G,NCReif 


7HlS**#E14«7t 


322 


PRINT 

CC 

nInt 

Wfilt 

prin! 

kRl 

PtHCH 

PRINT 

PiNCH 

NR 

PLNCF 

PRINT 


I21i (NCREGt 
122 NCBG«1#NCR£G 
lY«PGNTRYfNGR6) 


NCRBlt 


NTRY(NQRG},NlQRG«ltNOREGn 

NTRYfN0R61fNaRG«ltNQREGI» 

NTR¥fNQRG}fNORG»l«NQREG}) 




(Nl{N{NCR61»NCRGtI}tl»ltN|NTRY} 
(WUNfNCRBltNCRGtllfl^ltN 



XRfN 


NENTRYJ 
I 
1 

RYI 
^Ry| 
RTRY I 


611 

1€21 


FORPAH//l2)i2lHSTEBNHEIPEfl R-F ACTCR»«E 14«7//4X39tlTH| |XP£CTATICJN 
/ALEE uMCB’»*2Hl<UM»tfI4.7//29X14HGARA-lNF:NlT¥»,F14.7//l 
RRl-RRl+RRINCRetNEXITI 
PRINT 677t GRltRRl 

of r-factcr»*, 

|iRlNG-lDRANG4RR iNCRlf NiXITr 
SiGANG^SOGANG^GAPAlNCNCRBtREXlTJ 



1C22 

U13 

1C14 

C 

ICIS 

icie 

1C18 

ICll 

IClS 

IC2C 

1C23 

1C24 


R{:Ar<€»S8pK§4SRCSU8{t»Cf«BtKe](!T» 
^R|Ali€«SRfAA€4SR£SLI'lhCBBtNI)liTI 
6C_TC_717 

URRAC 4 Rfl|fcCRe,HEXni 
,U€RAC4GARAlR{hCR8tf«EXITi 
RCPIC 4 |RCSUMRCRetNiX!T» 
RERAC 4 SRESUMh£RetRlxn} 


liRRAO* 

St6RAC» 

SRCRAC* 

SRERAO* 

SUR^ 

iwp 

SRSt 
PRINT 


?-SljRARG-i 
T-SUGAR 64 
7 - SR* 
T-SBC 

sccc 


UPRAC 
yGRAC 
AR64SPERAC 
RC-iSRCRAC 


PRINT 8001 
PRInT !C13 
PRINT 1C13 
FORNATCSSh 
PRINT 9CCC 
|IN7 1C14 


4 4*4«44 4«4444*««4 49l444i44>ii4i|i«*3|3**9|4r4>*«3l43|443t3|i*>|9i>3|t4>4i4>t4E)ti I 
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FQRNAT ( Ix • 1 1 3ii44 44 4 444444 4444444444444444 4444444444444444444i4«44i44s 
X444444444444 4444 44 4 44 4 4 44 44 4 44 44 4444 4444 44 44444444444444 44444 I 


PRINT 1015 
FGRPATf lXtU3H4 
14 

PRINT 1018 
FQRNATdXf 
I • 

PRINT 1015 


4 4 4 1 

113H4LNPR-CBEITAL* P6RTUR6EC CBBITAt * 6ANNA-INFINITV 
F-CIRECT 4 F-iXCNANGI 4 R-TOTAL ♦! 


PRINT 1C14 
DO ICIO NC 

*f* 8 SSK!?§; 

lt2Xf IH 4 . 2 I 
NsNT 1015 
PRSnT 1014 
PRINT qccc 

■’gsj.ili' 

FCRPATla 


RE-IfNCCRE 

CNCRi} 



• NEX 



tF14.7«2XtlH4f2XtFI4.7f2Xt iMt 


NORB,N£Xnit 
Ixf lN4,2XtF14«7 


tSliORlCylRCRACf SRlRAOtSURRAO 

IXtFU« 7 « 2 XtlH*» 

1H4* ■ : ■ 




UX.FH.i 
IMiFUi 


^ - = if2XtlH=*,2X,F14.7t2X,lF*f2XtFl4«1t2XtlH*,ZX,Fl4,7.2X, 1F*| 

7. CC^UFUf 
223C CCM^MJI 
S7CP 
|«C 

ItERClTlNE FlRM/IVCLliPUtLMURtEXPMUNtHtJtRCRIEfMENTRVtNORef CELXRt22t 
1A2ERC} 


SUeRClTlRE IRTEGP CFURCIf hCREG# JENTRV,KCfCELXR» 
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C 

C 


NTEGRANDS FtiNC G 
“ USED* THE iNT“ 


tUlS SieRQGlIAf I^AiljXTfiS THE IRTIGRALS FUNCI CF THE 3N 
VARIuLrIsh PCIA-I. 7HE PETHCC CF FlNuTE CIFFEPENCES IS 
IS PERFORPEC tFFCUGH ICJCEFT INTERVALS ANC TFE FCURTH DIFfIrINCES INCL 

DIPENSICN JEFTPV(13),CELXRC13I 

DIpInSICN NFlfl3tA7)|XRI13,47l«WUN(16»12ff4A}tFUNC(iit4?}tA{4f5} 

CQPNCN liPLfXRtliUNfFLNCtA 

FLNCI«C-C 

DO 3| NCRG»l,NOREG 
FliNC!l»C.C 
FLNC!2-C.C 
FtlvCl3-C*C 
NgNTRV;J|NTPVfNCRGi 
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18 

32 

If 

2 C 


22 


33 

IS 

24 


INCRG-iri7 
tKCj 


17,35 




6,17 


XI 
IF 

N »1 

L*| 

Q« 2 » 

KRsQ 

GO TC IE 

K *3 

L*| 

Sl-l 

ft?CJlrFyNCfN^R 6 |K| 4 FUNCINGflG,K+l» 
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! FINEST 

FCNCfl-FUNCl 14FUAC CNCRG,4HFUNCiN0RGt J4l } 
F« 4 rAEN 7 RY 43 ) 2 C, 2 Ci 22 


A , 

J* J4l 

IF I4-hcr» ii« T'*3i 

F0NGl2«FUNCfNCRG,L43!4CC-2*l>»FUNClN0RG,L42|4( !•- 2«*C l♦FUNCf NORG, t+ 

/IHQfFONCINCPGjLI 
IFINENTRY-SI 33,30,33 
:F«RM |f,||,4l 
!fIn|nTRV-6I 26,26,38 

^ol€l2«FUNC124FONC(NCRGt443}«FUNCtNORG,J42HFUNCfNORGf441}4FyiiCIKO 

/R6,4) 



2« 


14 

!fi 


3C 


31 


j«4+ 1 

IF 


4“^ 1 

_ |J-^E^TFK■l41 24f24,26 
FtftCI3»FU^CUCfiG,6HfC-4,>>iFllftC«NCRG,5)-K6*-4.4QI*FUNClN0RG,41*(6. 
/♦C-4t l*FtJ^C«^CRGi3H ll#-4#’»C)#FuNCCNCRGt2HG*FUNCCNCRGf II 
IFIREMR'if-ei 34t3CfS4 

J ® 2 ' ' 

FUKCI3«FURC!34FURCIhCR€.J45)-3*4FUNC(HQRG.a+4l-»2«*FyNC(N0RG»J-»3l42 
/.♦FtRCCIiCPGf d42l-3.4FUhC(hCPCtJ4l|4F«jNCINCBGt4) 

4" j+ 1 

IFf4-«EMPV4SI 28,28,30 
FtjhCn-FURCil*CD£LXP^NCRGI)4ti,/2-l 
FtKCl2-FURC|24(0iLXRUCPGI|4|-l./24. I 


FliliC!3«FiRC!3*CDilXF<hCSGIMCl|,/1440, I 
FyRC!*FyNCl4FUNCll4FUNC224FURC!3 
CCHtlRyg 
RE1LRR 
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C 

C 
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|wc 

SifiRCyTlNE 


siFon (M 


SlRlLt^MCUS la^tlCF SCLVER 
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DIFERSXCN kFl.{13t47}|XRC13t41}tMUNfl6f 12,44lf FUNCC13t47 }fA{4»5) 
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